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1.  INTRODUCTION 


Artificial  graphitic  materials  have  been  used  for  more  than  the 
past  decade  in  reentry  vehicle  nosetips.  The  exacting  requirements  im- 
posed in  their  use  necessitate  accurate  stress  analysis  techniques.  An 
integral  part  of  every  stress  analysis  is  the  stress-strain  relationship 
or  material  model . 

The  stress  analysis  problems  inherent  to  reentry  vehicle  nosetip 
design  were  discussed  by  Jones  [1-1]  in  1967  along  with  numerous  specif- 
ic problems  by  other  authors  in  the  same  conference  proceedings  volume. 
Since  that  time,  periodic  reviews  of  nosetip  stress  analysis  technology 
have  been  made.  The  most  recent  review  by  Jones  and  Koenig  [1-2]  is 
addressed  to  the  material  modeling  characteristics  necessary  for  graphite 
and  carbon-carbon.  Two  of  the  significant  deficiencies  of  current  mate- 
rial modeling  that  they  point  out  are  (1)  biaxial  softening  and  (2)  dif- 
ferent moduli  under  tensile  loading  than  under  compressive  loading. 

These  characteristics  are  described  along  with  other  characteristics  of 
graphite  and  carbon-carbon  in  the  following  paragraphs. 

1.1  BIAXIAL  SOFTENING 

Biaxial  softening  Is  characterized  by  the  development  of  slightly 
larger  strains  in  biaxial  tension  titan  in  uniaxial  tension,  as  shown  in 
Figure  1-1.  This  behavior  of  generally  decreasing  Poisson's  ratios  is 
in  contradiction  to  what  might  be  anticipated  on  the  basis  of  conven- 
tional Poisson  effects  (where  v increases).  This  phenomenon  was  apparent- 
ly first  observed  by  Jortner  [1-3  thru  1-6]  for  graphite  and  is  attri- 
buted to  plastic  volume  changes  resulting  from  internal  tearing  or  micro- 
cracking.  Jones  and  Nelson  [1-7]  developed  a material  model  for  descrip- 
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FIGURE  1-1  BIAXIAL  SOFTENING  OF  GRAPHITE 

tlon  of  the  deformation  behavior  of  ATJ-S  graphite  under  biaxial  tension. 
Their  model  is  used  in  the  SAAS  III  program  to  obtain  predicted  strains 
for  Jortner's  biaxial  test  specimen  shown  in  Fig,  1-2.  The  predicted 
strains  are  shown  along  with  Jortner's  experimentally  observed  strains 
in  Figure  1-3  for  room  temperature  behavior  at  a constant  principal 
stress  of  3550  psi.  The  Jones  and  Nelson  strain  predictions  are  within 
3%  of  the  equal  biaxial  tension  strains  and  are  identical  to  the  two 
uniaxial  tension  cases. 

Actually,  the  Jones-Nelson  model  is  more  a general  model  for  non- 
linear behavior  of  orthotropic  materials  than  just  a biaxial  softening 
model.  Thus,  the  Jones-Nelson  model  should  be  considered  for  use  in 
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figure  1-2  HOLLOW  GRAPHITE  BIAXIAL  SPECIMEN 


FIGURE  1-3  BIAXIAL  STRAIN  RESPONSE  OF  A HOLLOW  ATJ-S  GRAPHITE  SPECIMEN 
AT  ROOM  TEMPERATURE  (70° F)  AND  3550  psi  PRINCIPAL  STRESS 
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modeling  other  materials.  Specifically,  carbon-carbon  will  be  shown  to 
be  a nonlinear  orthotropic  material,  and  the  basic  Jones-Nelson  model 
will  be  proposed  for  analysis  of  carbon-carbon. 

1.2  DIFFERENT  MODULI  IN  TENSION  AND  COMPRESSION 

Many  composite  materials  behave  differently  under  tensile  and  com- 
pressive loads.  Both  the  elastic  moduli  (stiffnesses)  and  the  strengths 
in  principal  material  property  directions  of  these  orthotropic  materials 
are  different  for  tensile  loading  than  for  compressive  loading.  This 
characteristic  behavior  is  shown  schematically  in  the  stress-strain 
curve  of  Figure  1-4.  This  phenomenon  is  but  one  of  several  differences 
that  make  composite  materials  more  difficult  to  analyze  (and  hence  de- 
sign) than  the  more  common  structural  materials  such  as  aluminum. 

Both  fiber-reinforced  and  granular  composite  materials  have  differ- 
ent moduli  in  tension  and  compression  as  displayed  in  Table  1-1.  Uni- 
directional glass  fibers  In  an  epoxy  matrix  have  compression  moduli  20% 
lower  than  the  tension  moduli  [1-8].  For  some  unidirectional  boron/epoxy 
fiber-reinforced  laminae,  the  compression  moduli  are  about  15-20%  larger 
than  the  tension  moduli  [1-9].  In  contrast,  some  unidirectional  graphite/ 
epoxy  fiber-reinforced  laminae  have  tension  moduli  up  to  40%  greater  than 
the  compression  moduli  [1-9].  Other  fiber-reinforced  composites  such  as 
carbon-carbon  have  tension  moduli  from  two  to  five  times  the  compression 
moduli  [1-10],  Thus,  no  clear  pattern  of  larger  tension  than  compression 
moduli  or  vice  versa  exists  for  fiber-reinforced  composite  materials.  A 
plausible  physical  explanation  for  this  puzzling  circumstance  has  yet  to 
be  made. 

For  granular  composite  materials,  the  picture  is  no  clearer.  ZTA 
graphite  has  tension  moduli  as  much  as  20%  lower  than  the  compression 

4 


,5 


TABLE  1-1 


TENSION  AND  COMPRESSION  MODULI  RELATIONSHIPS 
FOR  SEVERAL  COMMON  COMPOSITE  MATERIALS 


MATERIAL 

FIBROUS 

OR 

GRANULAR 

REPRESENTATIVE 

MODULI 

RELATIONSHIP 

GLASS/EPOXY 

FIBROUS 

Eta1*2Ec 

BORON/EPOXY 

FIBROUS 

Ec  - 1.2Et 

GRAPHITE/EPOXY 

FIBROUS 

Et"  ]*4Ec 

CARBON/CARBON 

FIBROUS 

Et  " 2-5Ec 

ZTA  GRAPHITE 

GRANULAR 

Ec"1*2Et 

ATO-S  GRAPHITE 

GRANULAR 

hm  1‘2Ec 

FIGURE  1-4  STRESS-STRAIN  CURVE  FOR  A MATERIAL  WITH 

DIFFERENT  MODULI  IN  TENSION  AND  COMPRESSION 
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moduli  [1-11].  On  the  other  hand,  ATJ-S  graphite  has  tension  moduli 
as  much  as  20%  more  than  the  compression  moduli  [1-12]. 

Many  other  materials  have  different  tension  and  compression  moduli. 
Which  modulus  is  higher  may  depend  on  the  fiber  or  granule  stiffness 
relative  to  the  matrix  stiffness.  Such  a relationship  would  influence 
whether  the  fibers  or  granules  tend  to  contact  and  hence  stiffen  the 
composite.  A general  physical  explanation  of  the  reasons  for  different 
behavior  In  tension  and  compression  Is  not  yet  available.  Investigation 
of  the  micromechanical  behavioral  aspects  of  composite  materials  may 
lead  to  a rational  explanation  of  this  phenomenon.  Until  such  an  expla- 
nation is  available,  the  apparent  behavior  can  be  used  in  analyzing  the 
stress-strain  behavior  of  materials.  That  is,  even  without  knowing  why 
the  materials  behave  as  they  do,  we  can  model  their  apparent  behavior. 

Actual  stress-strain  behavior  is  probably  not  as  simple  as  shown  in 
Figure  1-4.  Instead,  a nonlinear  transition  region  may  exist,  between 
the  tension  and  compression  linear  portions  of  the  stress-strain  [1-13]. 
The  measurement  of  strains  near  zero  stress  1,s  difficult  to  perform 
accurately,  but  the  stress-strain  behavior  might  be  as  shown  In  Figure 
1-5  wherein  replacement  of  the  actual  behavior  by  a bilinear  model  is 
offered  as  a simplification  of  the  obviously  nonlinear  behavior.  For 
most  materials,  the  mechanical  property  data  are  insufficient  to  justify 
use  of  a more  complex  material  model.  However,  one  possible  disadvantage 
of  the  bilinear  stress-strain  curve  approximation  is  that  a discontinuity 
in  slope  (modulus)  occurs  at  the  origin  of  the  stress-strain  curve. 

Given  that  the  uniaxial  stress-strain  behavior  is  approximated  by  a 
bilinear  representation,  the  definition  remains  of  the  actual  multiaxlal 
stress-strain,  or  constitutive,  relations  that  are  required  in  structural 
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FIGURE  1-5  COMPARISON  OF  BILINEAR  MODEL  WITH  ACTUAL  BEHAVIOR 

analysis.  Over  the  past  ten  years,  Ambartsumyan  and  his  co-workers 
[1-14  thru  1-17],  in  the  process  of  obtaining  solutions  for  stress  analy 
sis  of  shells  and  bodies  of  revolution,  defined  a set  of  stress-strain 
relations  that  will  be  referred  to  herein  as  the  Ambartsumyan  material 
model.  Jones  [1-18]  applied  the  model  to  the  problem  of  buckling  under 
biaxial  loading  of  circular  cylindrical  shells  made  of  an  isotropic  mate 
rial.  However,  in  application  of  the  Ambartsumyan  material  model  to 
orthotropic  materials,  certain  deficiencies,  such  as  a nonsymmetric 
compliance  matrix  In  the  stress-strain  relations  [1-19],  are  apparent. 

Jones  [1-20]  also  applied  modified  bilinear  stress-strain  relations 
to  buckling  of  shells  with  multiple  layers  of  orthotropic  materials  hav- 
ing different  moduli  in  tension  and  compression.  His  modifications  con- 
sist of  weighting  tension  and  compression  compliances  according  to  the 
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proportions  of  the  principal  stresses  in  order  to  obtain  a symmetric 
compliance  matrix.  Isabekian  and  Khachatryan  [1-21]  made  the  Ambartsumyan 
material  model  have  a symmetric  compliance  matrix  by  enforcing  certain 
relations  between  the  material  properties.  Both  Jones’  and  Isabekian 
and  Khachatryan' s relations  are  used  in  the  modified  Jones-Nelson  model, 
but  only  Jones'  weighted  compliance  matrix  material  model  is  used  in  the 
present  report. 

When  the  different  moduli  in  tension  and  compression  characteristic 
is  combined  with  the  biaxial  softening  characteristic,  the  Jones-Nelson 
material  model  leads  to  predicted  versus  experimental  strains  shown  In 
Figure  1-6.  There,  the  predicted  strains  In  the  mixed  tension  and  com- 
pression quadrant  are  within  3 to  9%  of  the  measured  strains  at  room 
temperature  for  the  data  shown.  Similar,  but  less  accurate,  results 
(from  9-12%  error)  are  shown  for  2000°F  in  Figure  1-7.  The  experimental 
data  In  Figure  1-7  are  much  less  accurate  than  the  data  in  Figure  1-6 
because  of  testing  difficulties  at  elevated  temperatures.  Thus,  the 
Jones-Nelson  graphite  material  model  is  validated  by  favorable  compari- 
son with  a well-defined  set  of  biaxial  experimental  data. 

1.3  CHARACTERISTICS  OF  GRAPHITE 

Graphites  used  in  reentry  vehicle  nosetips  are  macroscopical ly 
homogeneous,  transversely  Isotropic,  and  generally  fall  In  a brittle 
manner.  The  typical  stress-strain  curve  shown  in  Figure  1-8  is  non- 
linear to  failure.  A typical  initial  modulus  versus  temperature  rela- 
tionship is  also  shown  in  Figure  1-8.  There,  the  modulus  actually  in- 
creases from  Its  room  temperature  value  until  a temperature  of  about 
3500°F  is  reached  and  subsequently  decreases  to  nearly  zero  as  graphite 
approaches  sublimation.  In  addition,  at  all  temperatures,  the  axial 
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FIGURE  1-7  BIAXIAL  RESPONSE  OF  A HOLLOW  ATJ-S  GRAPHITE  SPECIMEN 
AT  2000° F AND  3550  psi  PRINCIPAL  STRESS 
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FIGURE  1-8  GRAPHITE  BEHAVIOR 

modulus  is  lower  than  the  modulus  in  the  circumferential  and  radial 
directions. 

Graphite,  as  mentioned  previously,  exhibits  the  biaxial  softening 
phenomenon  and  has  different  moduli  and  stress-strain  curves  in  tension 
than  in  compression.  These  characteristics  were  successfully  modeled 
by  the  Principal  Investigator  in  Air  Force  Contract  F3361 5-73-C-5124. 
However,  the  temperature-dependent  characteristic  has  not  yet  been 
coupled  to  the  other  characteristics  nor  has  an  actual  nosetip  been 
analyzed. 
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1.4  CHARACTERISTICS  OF  CARBON-CARBON 

Carbon-carbon  materials  used  in  reentry  vehicle  nosetips  are  macro- 
scopically  inhomogeneous  because  of  large  fibers  in  the  axial  direction 
of  the  nosetip.  These  materials  can  be  characterized  as  orthotropic 
if  the  fibers  are  in  orthogonal  directions,  but  are  anisotropic  if  fibers 
at  other  than  90°  angles  are  inserted.  Carbon-carbon  fails  in  a pro- 
gressive manner  as  illustrated  in  Figure  1-9.  There,  the  material  is 
stressed  in  the  direction  of  axial  fibers  which  apparently  slip  relative 
to  the  matrix  material  as  stress  is  applied.  The  Initial  modulus  versus 
temperature  relationship  is  also  shown  for  the  circumferential  and  radial 
directions  in  Figure  1-9.  The  three  curves  shown  are  interpretations  of 
the  same  experimental  data  by  different  people.  Thus,  considerable  dis- 
agreement exists  as  to  the  actual  modulus  versus  temperature  relation- 
ship. The  axial  modulus  for  this  particular  carbon-carbon  material  is 
about  twice  the  circumferential  modulus.  Such  a relation  (quite  dif- 
ferent from  graphite)  is  not  unexpected  when  the  large  fibers  in  the 
axial  direction  are  considered. 

Carbon-carbon,  like  graphite,  exhibits  different  moduli  in  tension 
than  in  compression;  however,  the  differences  are  strikingly  greater 
for  carbon-carbon  than  for  graphite.  No  evidence  currently  exists  that 
carbon-carbons  exhibit  the  biaxial  softening  phenomenon.  Some  of  the 
specific  mechanical  properties  of  carbon-carbon  are  given  by  Legg, 

Starrett,  Sanders,  and  Pears  [1-22]. 

The  most  evident  difference  of  carbon-carbon  from  graphite  is  its 
three-dimensional  woven  character  as  opposed  to  the  fine-grained  struc- 
ture of  graphite.  The  fibers  in  carbon-carbon  are  placed  in  three  mutu- 
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FIGURE  1-9  CARBON-CARBON  BEHAVIOR 
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ally  perpendicular  directions.  Thus,  carbon-carbon  is  a highly  ortho- 
tropic  material  in  the  r-6  plane  of  a nosetip  (as  opposed  to  the  Isotropy 
of  graphite  in  this  plane).  Geiler  [1-23]  used  a linear  elastic  model 
in  the  ASAAS  program  due  to  Crose  [1-24]  to  account  for  the  circumfer- 
entially varying  orthotropy.  Geiler  obtained  apparently  good  results. 
However,  the  ASAAS  program  would  be  very  difficult  to  adapt  to  nonlinear 
analysis  because  of  the  already  highly  coupled,  time  consuming  internal 
workings  of  the  computer  program. 

A related  characteristic  of  carbon-carbon  is  that  the  size  of  the 
fibers  is  not  negligible  in  comparison  to  the  size  of  the  billets  in 
which  it  is  manufactured  or  in  comparison  to  the  critical  dimensions  of 
the  nosetip  into  which  the  billets  are  machined.  (Note  the  Z-direction 
fibers  are  the  white  streaks  in  the  axial  direction  of  the  nosetip  in 
Figure  1-10,)  Another  way  of  saying  the  same  thing  is  that  the  fiber 
spacing  is  of  the  same  order  of  magnitude  as  the  distance  over  which  the 


FIGURE  1-10  CARBON-CARBON  PLUG  NOSETIP 
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stresses  change  rapidly.  Thus,  we  would  anticipate  possible  difficulties 
in  applying  a macromechanical  or  continuum  mechanics  model  to  carbon- 
carbon  materials.  Not  enough  work  has  been  done,  however,  to  resolve 
or  even  clarify  the  macromechanics  versus  micromechanics  issue. 

1.5  STATEMENT  OF  THE  PROBLEM 

Two  principal  efforts  are  involved  In  this  research:  one  on  graph- 
ite material  modeling  and  the  other  on  carbon-carbon  material  modeling. 

The  graphite  modeling  is  essentially  a continuation  of  efforts  begun 
under  Air  Force  Contract  F33615-73-C-5124.  In  that  contract,  graphite 
multiaxial  stress-strain  behavior  was  successfully  modeled  under  both 
biaxial  tension  and  mixed  tension  and  compression  load  at  room  tempera- 
ture and  at  2000°F.  The  remaining  tasks  include  (1)  incorporating  a 
temperature-dependent  character  in  the  material  model  (previously  men- 
tioned results  are  for  a constant  temperature),  (2)  validating  the  model 
by  comparison  with  further  experimental  data,  and  (3)  exercising  the 
model  in  a thermostructural  analysis  of  an  actual  reentry  vehicle  nose- 
tip.  Upon  completion  of  these  tasks,  the  graphite  material  model  should 
be  ready  for  routine  use  in  Air  Force  reentry  vehicle  nosetip  analysis. 

The  carbon-carbon  modeling  is  a new  effort  with  the  objective  of 
applying  the  basic  concepts  of  the  successful  graphite  model  to  the  analy- 
sis of  carbon-carbon  stress-strain  behavior.  Generally,  carbon-carbon 
stress-strain  curves  are  more  jagged  than  those  of  graphite.  Thus,  some 
modifications  to  the  graphite  model  are  anticipated.  The  first  step  In 
carbon-carbon  modeling  is  to  describe  and  evaluate  the  stress-strain 
curve  characteristics.  Next,  a revised  model  will  be  formulated  based 
on  these  characteristics  and  on  discussions  with  researchers  who  have 
been  dealing  with  carbon-carbon  for  some  time  (at  Southern  Research 

15 


li Hi ' ff ii iMMu » m ■ • *#'-* ■ ■*»»►»•  ,\  * Iff* u,  -i*  11  > 

„■  ....  ■■■  


3 


1 r 


I 


Institute,  McDonnell -Douglas,  AFML,  SAMSO,  Prototype  Development  Asso- 
ciates, arid  Weiler  Research,  Inc.).  The  model  will  then  be  correlated 
with  available  experimental  data  in  a validation  stage.  If  the  model 
does,  or  can  be  refined  enough  to,  give  good  correlation  with  experi- 
mental data,  then  the  model  will  be  evaluated  for  implementation  in  AFML 
and  SAMSO  nosetip  thermostructural  analysis  computer  programs. 

1.6  STATEMENT  OF  RESEARCH 

The  present  research  is  divided  in  two  najor  phases,  graphite  and 
carbon-carbon,  each  of  which  are  further  divided  as  follows: 

Phase  G - Graphite 

G-I  - Model  Formulation 
G-II  - Correlation 
G-III  - Nosetip  Demonstration 
G-IV  - Implementation 
G-V  - Reporting 
Phase  C - Carbon-Carbon 

C-I  - Data  Evaluation 
C-II  - Model  Formulation 
C-I 1 1 - Correlation 
C-IV  ••  Implementation 
C-V  - Reporting 

These  phases  are  described  in  the  following  paragraphs. 

1.6.1  PHASE  G - GRAPHITE 

Phase  G-I  - Model  Formulation 

The  graphite  material  model  has  but  one  essential  characteristic  to 
be  Incorporated  prior  to  use  in  actual  nosetip  analysis.  That  charac- 
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1:'  teristic  is  temperature-dependent  material  behavior.  A new  scheme  must 

r;  be  devised  to  express  the  material  model  as  a function  of  temperature 

based  on  data  at  a finite  number  of  temperatures.  Basically,  the  objec- 

j tive  is  a material  property  versus  temperature  interpolation  scheme. 

[■  However,  this  scheme  is  complicated  by  the  presence  of  many  more  mate- 

! 

[;■  rial  property  characterization  constants  in  the  present  model  than  in 

[' 

i.  ; previous  models  for  which  such  interpolation  is  well-known. 

I;'  Phase  G- 1 1 - Correlation 

The  graphite  material  model  will  continue  to  be  correlated  with 
■'i;  available  experimental  data.  All  known  biaxial  data  generated  on  con- 

ji.  stant  temperature  tube  specimens  by  Jortner  of  McDonnell-Douglas  has  been 

: successfully  correlated.  The  next  logical  step  is  to  attempt  correla- 

I ' tion  with  data  on  specimens  with  a nonconstant  temperature,  i.e.,  a 

t •- 

| temperature  gradient.  The  data  for  the  wedge-shaped  disc  shown  In  Figure 

j 1-11  generated  in  the  Temperature/Stress  Test  developed  by  Southern 
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Research  Institute  [1-25]  is  most  appropriate.  Limited  correlation 
studies  will  be  performed  with  that  data.  Calculated  temperature  pro- 
filges  (verified  by  measurement)  through  the  disc  will  be  used  in  con- 
junction with  the  temperature-dependent  material  model  to  predict 
disc  diameter  changes.  These  predictions  will  be  compared  with  Southern 
Research  Institute  measurements  to  further  validate  the  material  model. 

Phase  G-III  - Nosetip  Demonstration 

The  final  stage  in  the  development  of  the  material  model  Is  to  exer- 
cise the  model  in  thermostructural  analysis  of  a reentry  vehicle  nose- 
tip.  Although  "experimental"  data  are  limited  (obviously  only  flight 
tests  with  limited  Instrumentation  can  be  considered  experiments  in  the 
present  context),  the  comparison  of  predictions  of  the  present  model  with 
previously  used  models  is  essential.  Without  such  comparisons,  the  real 
worth  of  the  present  graphite  material  model  for  reentry  vehicle  nosetip 
stress  analysis  cannot  be  established. 

The  finite  element  data  cards  for  the  nosetip  and  its  thermal  (and 
mechanical)  loading  will  be  chosen  by  AFML  and  supplied  to  Southern  Meth- 
odist University  In  the  SAAS  III  format.  These  data  cards  shall  have 
been  previously  verified  to  work  on  SAAS  III  by  another  contractor  for 
an  all-elastic  analysis  or  an  elastic-plastic  bilinear  analysis.  Accord- 
ingly, the  nosetip  demonstration  will  involve  a single  variable,  the 
material  model.  All  other  variables  will  be  provided  to  SMU  in  ready- 
to-run  form. 

Phase  G-IV  - Implementation 

The  graphite  material  model  developed  and  validated  in  previous 
phases  and  in  Air  Force  Contract  F33615-73-C-5124  will  be  incorporated 
in  a version  of  the  SAAS  III  computer  program  [1-26].  That  program  is 
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the  basic  finite  element  computer  program  operational  at  Southern 
Methodist  University.  Duplicate  computer  decks,  listings,  test  cases, 
and  output  will  be  provided  to  AFML  should  an  operational  deck  be  de- 
sired. In  view  of  Southern  Methodist  University's  mission  of  graduate 
education  and  research,  limited  manpower,  and  small  computer,  the  graph- 
ite material  model  will  not  be  implemented  in  any  program  other  than 
SAAS  III. 

Phase  G-V  - Reporting 

The  graphite  material  modeling  efforts,  correlation  studies,  and 
nosetip  analysis  are  described  in  the  present  Technical  Report. 

1.6.2  PHASE  C - CARBON-CARBON 

Phase  C-I  - Data  Evaluation 

Carbon-carbon  materials  have  many  different  manufacturing  processes 
and  hence  many  different  characteristics,  as  alluded  to  earlier.  The 
objective  of  the  data  evaluation  phase  is  to  examine  the  available  mate- 
rial property  data  and  Isolate  the  significant  characteristics  that  must 
be  modeled  to  accurately  predict  thermostructural  response.  This  objec- 
tive will  be  met  by  review  of  published  data,  review  of  published  mate- 
rial modeling  efforts,  and  consultation  with  AFML,  SAMSO,  Southern  Re- 
search Institute,  McDonnell-Douglas,  Prototype  Development  Associates, 
and  Weller  Research,  Inc.  The  latter  consultation  should  take  place  both 
In  this  phase  and  In  subsequent  phases.  Because  of  this  consultation 
and  the  expectation  that  important  material  properties  will  likely  be 
found  to  not  have  been  measured,  data  evaluation  will  be  regarded  as  a 
phase  continuing  throughout  the  remainder  of  the  program. 

Phase  C-II  - Model  Formulation 

The  basic  material  model  used  for  graphite  and  described  In  Reference 
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1-7  will  be  fit  to  the  thermostructural  characteristics  of  carbon-carbon. 
Should  a more  complicated  relationship  be  necessary,  studies  will  be  ini- 
tiated to  determine  an  appropriate  relationship. 

The  behavior  of  carbon-carbon  will  be  studied  in  this  initial  effort 
from  the  standpoint  of  the  axi symmetric  macromechanics  or  continuum  me- 
chanics analysis  in  SAAS  III  [1-26].  In  addition  to  the  basic  SAAS  III 
analysis  with  the  new  material  model,  an  essentially  one-element  model 
with  the  new  material  model  will  be  used  in  modeling  feasibility  studies 
such  as  were  carried  out  in  Air  Force  Contract  F33615-73-C-5124.  Neither 
the  asymmetric  linear  elastic,  analysis  of  ASAAS  [1-24]  nor  micromechan- 
ical analyses  will  be  attempted.  Evaluation  of  the  need  for  and  benefit 
from  more  geometrically  sophisticated  models  will  be  made  In  conjunction 
with  the  data  evaluation  phase.  That  Is,  we  must  be  certain  the  quality 
and  sophistication  of  the  data  merit  the  Increased  expense  of,  e.g., 
asymmetric  micromechanical  analyses. 

Carbon-carbon  will  be  modeled  in  uniaxial  on-axis  and  off-axis 
stress  states  for  which  material  property  data  exist.  Moreover,  limited 
biaxial  stress  state  data  will  also  be  modeled  to  the  level  of  sophisti- 
cation possible  within  the  scope  of  two-dimensional  or  axi symmetric  macro 
mechanics  theory. 

Phase  C-III  - Correlation 

The  material  model  developed  in  the  previous  phase  will  be  used  to 
obtain  stress-strain  predictions  for  situations  In  which  carefully  ob- 
tained experimental  results  are  available.  These  experimental  data  will 
be  selected  in  cooperation  with  AFML.  These  experiments  will  Include 
uniaxial  on-axis  and  off-axis  tests,  the  most  logical  starting  point,  for 
any  material  modeling  effort.  In  this  manner,  the  material  model  will 
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either  be  validated  by  comparison  with  experimental  data  or  will  be  in- 
validated and  Improvements  will  be  made. 

Phase  C-IV  - Implementation 

The  carbon-carbon  material  model  developed  and  validated  in  previous 
phases  will  be  incorporated  in  a version  of  the  SAAS  III  computer  pro- 
gram [1-26].  That  program  is  the  basic  finite  element  computer  program 
operational  at  Southern  Methodist  University.  This  research  program  is 
not  likely  to  advance  the  state  of  the  art  of  carbon-carbon  modeling  to 
the  point  where  the  model  Is  judged  completely  ready  for  widespread  Air 
Force  use  in  reentry  vehicle  nosetip  stress  analysis.  Instead,  the 
present  effort  is  best  described  as  a bold  step  toward  that  goal  with 
some  hope  of  a reasonable  model  being  obtained  within  the  next  year. 
Southern  Methodist  University  expects  that  further  work  on  carbon-car- 
bon modeling  will  be  necessary  before  the  real  question  of  Implementation 
of  the  carbon-carbon  model  in  Air  Force  computer  programs  arises. 

The  carbon-carbon  material  modeling  studies,  including  data  evalua- 
tion and  characterization,  model  formulation,  and  correlation  activities, 
are  described  in  the  present  report. 

1.7  SCOPE  OF  REPORT 

The  actual  accomplishments  during  the  contract  are  presented  in  the 
following  sections.  First,  the  Jones-Nelson  nonl Inear  material  model  is 
reviewed  in  Section  2.  Also,  that  model  is  extended  to  temperature-de- 
pendent material  behavior  and  to  treatment  of  extrapolated  stress-strain 
curves  in  what  is  called  the  Uones-Nel soi  -Morgan  model.  Then  in  Section 
3,  the  graphite  modeling  efforts  are  described,  first,  the  correlation 
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studies  for  the  Southern  Research  Institute  thermal  stress  disk  test  are 
discussed.  Then,  the  AFFDI.  50  MW  nosetip  correlation  studies  are  de- 
scribed. The  carbon-carbon  modeling  efforts  are  discussed  in  Section  4. 
The  general  modeling  is  first  described  and  then  a characteristic  of 
carbon-carbon  in  bending  tests  is  treated.  The  current  contract  efforts 
are  summarized  in  Section  5. 
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2,  JONES-NELSON-MORGAN  NONLINEAR  MATERIAL  MODEL 

2.1  INTRODUCTION 

First,  the  basic  Jones-Nelson  nonlinear  material  model  is  briefly 
reviewed  in  Section  2.2,  and  some  new  aspects  of  its  behavior  are  de- 
scribed. In  particular,  procedures  are  developed  for  temperature  Inter- 
polation of  temperature-dependent  material  behavior.  Then,  methods  of 
extending  the  range  of  applicability  of  the  material  model  are  discussed 
in  Section  2.3.  The  major  accomplishment  in  that  section  is  the  develop- 
ment of  an  extended  stress-strain  curve  version  of  the  Jones-Nelson  model 
which  is  called  the  Jones-Nelson-Morgan  nonlinear  material  model.  Final- 
ly, the  JNMDATA  computer  program  is  described  in  Section  2.4.  This  pro- 
gram Is  used  to  convert  the  various  measured  stress-strain  curve  data 
directly  to  the  parameters  of  the  Jones-Nelson  and  Jones-Nelson-Morgan 
material  models  suitable  for  use  in  the  SAAS  HIM  finite  element  stress 
analysis  computer  program  which  is  a modification  of  the  SAAS  III  pro- 
gram [2-1].  The  JNMDATA  program  is  a very  useful  aid  in  the  modeling  of 
a material  because  the  results  are  obtained  automatically  and  are  pre- 
sented visually  for  rapid  evaluation  of  the  model. 

2.2  JONES-NELSON  NONLINEAR  MATERIAL  MODEL 

The  Jones-Nelsori  nonlinear  material  model  was  developed  under  USAF 
Contract  F3361 5-73-C-51 24  and  reported  in  AFML-TR-74-259  [2-2].  Several 
other  related  and  more  accessible  publications  are  condensed  from  Ref. 
2-2,  namely  Refs.  2-3  thru  2-6.  That  work  will  be  summarized  in  Section 
2.2.1  for  the  sake  of  convenience  in  reading  this  report.  The  necessary 
further  details  will  be  referenced  where  required.  Then,  a new  discus- 
sion of  how  to  Implement  the  model  is  presented  in  Section  2.2.2.  There, 
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the  recent  experience  In  application  of  the  model  is  reflected.  Finally, 
the  model  is  extended  in  Section  2.2.3  to  interpolation  of  temperature- 
dependent  material  behavior  at  temperatures  between  available  data. 

2,2.1  BASIC  APPROACH 

The  basic  problem  is  the  stress  analysis  of  nonlinear  elastic  bodies 
whose  stress-strain  behavior  is  described  with,  for  the  example  of  an 
orthotropic  axisymmetric  body  under  axisymmetric  load,  the  equations: 
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where  the  directions  denoted  with  the  subscripts  r,  z,  and  0 are  princi- 
pal material  directions.  The  material  properties  in  the  compliance  ma- 
trix of  Eq.  (2.1)  are 

Er  = Young's  modulus  in  the  r direction 

Ez  = Young's  modulus  in  the  z direction 

Eq  = Young's  modulus  in  the  0 direction 

\>rz  - -ez/er  for  the  loading  or  = a (all  other  stresses  zero) 

vro  = “eQ /Er  ^or  i°ac*in9  ar  ~ a (a^  other  stresses  zero) 

vze  = "eG^ez  <r°r  tfie  loading  oz  3 a (all  other  stresses  zero) 

Grz  = Shear  modulus  in  the  rz  plane 
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The  reciprocal  relations  of  orthotropic  elasticity 


vrz^r  ~ vzr^Ez 


Vr(/Er  v0r^E0 


vz0^Ez  = voz/Ee  '2t2^ 


can  be  used  to  express  alternative  definitions  for  the  Poisson's  ratios 
in  terms  of  the  seven  independent  material  properties  in  Eq.  (2.1).  The 
material  properties  in  Eq.  (2.1)  are  a function  of  stress  level  because 
the  material  is  nonlinearly  elastic.  However,  we  do  not  examine  the 
unloading  behavior  nor  any  subsequent  reloading  behavior  of  the  body. 

The  basic  stress  analysis  problem  could  be  more  complicated  than  is 
represented  with  Eq.  (2.1).  For  example,  the  material  could  have  prin- 
cipal material  directions  at  some  angle  to  the  r-z-e  coordinate  system. 

Or,  the  material  could  have  different  stress-strain  behavior  in  tension 
than  in  compression.  The  stress-strain  relations  for  multimodulus  mate- 
rials are  derived  by  Jones  and  Nelson  [2-2,  2-5]  and  are  applied  to  stress 
analysis  of  graphitic  materials  in  Ref.  2-6.  The  foregoing  considerations 
are  obviously  more  complicated  than  what  is  represented  with  Eq.  (2.1). 
However,  those  equations  will  suffice  for  our  discussion  of  the  Jones- 
Nelson  material  model. 

The  basic  premise  of  the  Jones-Nel son  nonlinear  material  model  Is 
that  the  mechanical  properties  of  a material,  e.g.,  the  material  prop- 
erties in  Eq.  (2.1),  are  expressed  in  terms  of  the  strain  energy  of  the 
body  with  the  approximate  equation 

Ci 

Mechanical  Property ^ = A^[l  - (U/UQ  ) ] (2.3) 

where  the  are  the  elastic  values  of  the  material  property,  the  and 
C.  are  related  to  the  Initial  curvature  and  rate  of  change  of  curvature, 
respectively,  of  the  stress-strain  curve  [2-2,  2-4]  (slightly  different 
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interpretations  exist  when  the  mechanical  property  is  a Poisson's  ratio), 
and  U is  the  strain  energy  density  of  an  equivalent  elastic  system  at 
each  stage  of  nonlinear  deformation: 

U = (o  e + oc  + o g + t Y )/2  (2.4) 

v r r z z 00  rz'rz  k ' 

The  strain  energy  density  U is  normalized  by  U0  in  Eq.  (2.3)  so  that  B.. 
and  are  dimensionless.  Typical  stress-strain  curves,  the  correspond- 
ing mechanical  property  versus  strain  energy  curves,  and  the  associated 
mechanical  property  equations  for  a Young's  modulus  and  a Poisson's  ratio 
are  shown  in  Fig.  2-1. 

The  nonlinear  stress-strain  model  is  actually  much  more  complicated 
than  Eq.  (2.3).  When  mixed  tensile  and  compressive  stresses  are  excited, 
the  strain  energy  used  in  Eq.  (2.3)  could  be  a weighted  combination  of 
the  strain  energy  of  compression  and  that  of  tension.  Moreover,  all 
coefficients  have  different  values  In  tension  than  in  compression.  The 
choice  of  which  properties,  tension  or  compression,  should  be  used  is 
made  in  the  Ambartsumyan  superposition  manner  [2-7]  after  rotating  the 
stress-strain  relations  to  principal  stress  directions  as  described  in 
Ref.  2-2. 

The  stress-strain  relations,  Eq.  (2.1),  and  the  mechanical  property 
versus  energy  equations,  Eq.  (2.3),  are  a set  of  indeterminate  relations 
which  are  solved  with  the  iteration  procedure  shown  in  Fig.  2-2.  That 
is,  the  stresses  and  strains  depend  on  the  mechanical  properties  (through 
Eq.  2.1)  which,  in  turn,  depend  on  the  stresses  and  strain  [through  Eqs. 
(2.3  and  (2.4)].  The  determination  of  mechanical  properties  and,  con- 
sequently, the  stress-strain  relationships  is  based  on  both  the  propor- 
tions of  the  principal  stresses  and  on  the  magnitude  of  an  energy  func- 
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(a)  DIRECT  MODULUS 


AXIAL  STRESS 


(b)  POISSON'S  RATIO 


FIGURE  2-1  REPRESENTATION  OF  STRESS-STRAIN  RELATIONS 
FOR  DIRECT  MODULI  AND  POISSON'S  RATIOS 
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| EXPRESS  MATERIAL  PROPERTIES  IN  TERMS  OF  TOTAL  ENERGY  U FROM  UNIAXIAL  DATA 

FORM  COMPLIANCE  MATRIX  IN  PMD  WITH 
INITIAL  LINEAR  TENSION  VALUES  OF  MATERIAL  PROPERTIES 

I 


CALCULATE  STRESSES,  STRAINS,  AND  STRAIN  ENERGY 
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1 

CALCULATE  NEW  MATERIAL  PROPERTIES 

h 

| FORM  ALL-TENSION  AND  ALL-COMPRESSION  COMPLIANCE  MATRICES  IN  PMD  | 

ROTATE  COMPLIANCE  MATRICES  TO  PSD 

FORM  MULTIMODULUS  COMPLIANCE  MATRIX  IN  PSD  | 




| CALCULATE  NEW  STRESSES,  STRAINS,  AND  STRAIN  ENERGY  | 

EVALUATE  REL  AU  ■ (U^  ] 

i i ’ 

I F t REL  AU | < SPECIFIED  VALUE,  STOP 

I | REL  AU | > SPECIFIED  VALUE 

PMO  - PRINCIPAL  MATERIAL  DIRECTIONS 
PSD  ■ PRINCIPAL  STRESS  DIRECTIONS 


FIGURE  2-2  ITERATION  PROCEDURE  FOR  NONLINEAR  MULTIMODULUS  MATERIALS 
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tion.  Each  step  in  the  iteration  procedure  is  described  in  Refs,  2-2 
and  2-5. 

Two  different  energy  functions  - total  strain  energy  and  weighted 

strain  energy  - can  be  used  in  the  Jones-Nelson  material  model.  The 

total  strain  energy  is  defined  in  Eq.  (2.4).  On  the  other  hand,  in  the 

weighted  strain  energy,  the  total  strain  energy  is  separated  into  two 

components:  (1)  the  contribution  from  the  tensile  principal  stresses 

and  (2)  the  contribution  from  the  compressive  principal  stresses.  Then, 

the  effective  energy  level,  u,  in  terms  of  the  tension  and  compression 

w 

components  of  the  total  strain  energy  is 

Uw  = (Ut2  + yc2)/U  (2>5) 

This  energy  U,,  is  used  to  determine  both  the  tension  and  compression 
material  properties. 

2.2.2  IMPLEMENTATION  OF  THE  MATERIAL  MODEL 

Much  of  the  work  in  Implementing  the  Jones-Nelson  nonlinear  mate- 
rial model  is  in  calculating  appropriate  values  of  A,  B,  and  C in  the 
governing  equation,  Eq.  (2.3).  To  reach  the  point  where  these  calcula- 
tions can  be  made,  we  must  first  determine  the  secant  values  of  the  me- 
chanical properties  and  the  corresponding  values  of  strain  energy  from 
uniaxial  stross-strain  curves  of  the  material  under  investigation.  For 
example,  the  values  of  the  secant  shear  modulus  and  corresponding  values 
of  U can  be  determined  from  a shear  stress  - shear  strain  curve  such  as 
shown  in  Fig.  2-3.  There,  for  two  shear  stress  levels,  the  correspond- 
ing shear  strains  are  found  from  the  experimentally  determined  shear 

stress  - shear  strain  curve.  Then,  G1?  and  U are  calculated  from 

sec 
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FIGURE  2-3  NONLINEAR  SHEAR  STRESS  - SHEAR  STRAIN  CURVE 


(2.6) 


U “ T 


12Y12/2 


(2.7) 


Similarly,  the  remaining  secant  mechanical  properties  and  their  corres- 
ponding strain  energies  are  found  from  the  appropriate  stress-strain 
curves  In  principal  material  directions. 

Next,  all  secant  mechanical  properties  are  plotted  versus  the  strain 
energy  as  in  Fig.  2-4.  Now,  we  can  begin  to  calculate  or  otherwise  find 
the  values  of  A,  B,  and  C In  Eq.  (2.3).  First,  the  constant  A Is  the 
Initial  (elastic)  value  of  the  mechanical  property.  That  is,  It  is  the 
Initial  slope  of  the  stress-strain  curve  In  Fig.  2-3  or  the  intercept  of 
the  mechanical  property  versus  energy  curve  in  Fig.  2-4.  The  value  of  A 
Is  higher  than  you  would  expect  from  the  normal  procedure  of  placing  a 
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FIGURE  2-4  REPRESENTATIVE  MECHANICAL  PROPERTY  VERSUS  U CURVE 


straightedge  on  a stress-strain  curve  such  as  Fig.  2-3.  We  draw  this 
conclusion  on  the  basis  of  many  observed  fits  of  the  material  model  to 
stress-strain  curves  with  simultaneous  fits  to  a corresponding  mechani- 
cal property  versus  energy  curve.  That  Is,  the  curve  on  a plot  of  me- 
chanical energy  versus  energy  "heads"  for  a much  higher  intercept  at 
U-0  than  you  would  expect  from  looking  at  the  stress-strain  curve  alone. 
This  conclusion  Is  not  surprising  when  viewed  In  the  context  that  the 
mechanical  property  is  the  slope  (first  derivative)  of  the  stress-strain 
curve  and  hence  Is  a more  sensitive  indicator  of  the  behavior  when 
plotted  against  energy  than  when  visually  determined  from  a necessarily 
somewhat  inaccurate  stress- strain  curve. 

The  constants  B and  C are  determined  from  data  at  two  points  on  the 
material  property  versus  U curve  as  shown  In  Fig.  2-4.  The  values  of  the 
mechanical  property  and  the  values  of  U at  these  two  data  points  are  sub- 
stituted in  Eq.  (2.3).  Two  equations  in  the  two  unknowns,  B and  C,  re- 
sult. These  two  equations  are  solved  simultaneously  for  B and  C to  get 
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(2.8) 


C = 


log 


a-(m.p.)2“ 
_a-Tm.p.)3  _ 


a-(m.p.)2/u2 

B ■ k 


-c 


(2.9) 


where  (M.P.)2  and  (M.P.)3  are  the  values  of  the  mechanical  property  at 
the  two  data  points  chosen  from  the  mechanical  property  versus  U curve. 
The  values  U2  and  U3  in  Eqs.  (2.8)  and  (2.9)  are  the  strain  energies  at 
the  two  chosen  data  points.  The  proper  choice  of  data  points  from  the 
mechanical  property  versus  U curve  is  an  Important  part  In  the  applica- 
tion of  the  material  model  to  a specific  material  and  will  be  discussed 
next. 

The  use  of  Eq.  (2.3)  to  approximate  each  mechanical  property 
corresponds  to  the  mathematical  procedure  of  three-point  Interpolation. 

At  the  three  data  points  used  in  determining  A,  B,  and  C,  the  approximate 
values  of  the  mechanical  property  are  the  same  as  the  actual  values.  In 
the  region  between  these  three  points,  Eq.  (2.3)  Is  a reasonable  approxi- 
mation of  the  mechanical  property  versus  strain  energy  data.  Thus,  an 
Interval  of  strain  energy  for  which  Eq.  (2.3)  is  a valid  representation 
of  the  mechanical  property  Is  defined  by  the  positions  of  points  1 and 
3 in  K1g.  2-4.  This  Interval  of  strain  energy  is  bounded  on  the  left  by 
U=0,  the  strain  energy  at  point  1 , and  on  the  right  by  the  strain  energy 
at  point  3.  The  approximate  mechanical  property  versus  U curve  must  pass 
through  point  2 so  the  shape  of  this  curve  between  points  1 and  3 is  de- 
termined from  the  position  of  point  2.  Obviously,  more  information  than 
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just  data  to  calculate  B and  C is  gained  about  the  mechanical  property 
approximation  from  the  data  points  1,  2,  and  3 in  Fig.  2-4.  The  func- 
tions of  these  points  can  be  summarized  as  follows: 

(a)  Point  1 is  used  to  determine  the  constant  A in  Eq.  (2.3)  and 
is  the  left  bound  of  the  interval  of  strain  energy  for  which 
the  approximate  mechanical  property  versus  U curve  is  valid. 

(b)  Point  3 is  used  in  determining  B and  C in  Eq.  (2.3)  and  is 
the  right  bound  of  the  interval  of  strain-energy  for  which  the 
approximate  mechanical  property  versus  U curve  is  valid. 

(c)  Point  2 is  also  used  in  determining  B and  C in  Eq.  (2.3),  and 
the  shape  of  the  approximate  mechanical  property  versus  U curve 
Is  defined  by  the  position  of  point  2 in  the  interval  between 
points  1 and  3. 

The  choice  of  data  points  used  to  determine  B and  C is  quite  impor- 
tant. If  B and  C are  calculated  from  data  at  points  2 and  3 on  the  actual 
mechanical  property  versus  strain  energy  curve  (solid  line)  in  Fig.  2-5, 
the  dashed-dotted  curve  labeled  2-3  is  the  result.  This  curve  is  a rea- 
sonable representation  of  the  actual  data  in  the  Interval  0 < U < U3 
where  U3  is  the  value  of  the  strain  energy  at  point  3.  For  strain  ener- 
gies larger  than  U3,  the  "2-3"  curve  does  not  and  is  not  supposed  to  re- 
present the  mechanical  property  versus  strain  energy  behavior  accurately. 
The  mechanical  property  equation  is  valid  over  a larger  Interval  of 
strain  energy  if  a point  further  out  than  point  3 on  the  actual  mechani- 
cal property  versus  U curve  is  used  in  the  calculation  of  B and  C.  For 
example,  when  point  4 is  used  in  determining  B and  C,  Eq.  (2.3)  Is  valid 
In  the  interval  0 < U < where  U4  is  the  value  of  the  strain  energy 
at  point  4.  Both  the  "2-4"  curve  (B  and  C determined  from  data  at  points 
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FIGURE  2-5  REPRESENTATIVE  MECHANICAL  PROPERTY  VERSUS  U BEHAVIOR 
AND  POSSIBLE  APPROXIMATIONS 

2 and  4)  and  the  "3-4"  curve  (B  and  C determined  from  data  at  points  3 
and  4)  In  Fig.  2-5  are  reasonable  approximations  of  the  actual  mechanical 
property  versus  strain  energy  data  for  0 < U < U4.  However,  in  the  in- 
terval 0 < U < U3  neither  the  "2-4"  nor  the  "3-4"  curve  are  as  good  a 
representation  of  the  actual  data  as  the  "2-3"  curve.  Thus,  the  me- 
chanical property  equation  Is  valid  over  a large  Interval  of  strain  ener- 
gy when  point  4 Is  used  in  finding  B and  C,  but  the  accuracy  of  the  ap- 
proximation In  smaller  subintervals  of  the  overall  Interval  is  sacrificed. 

Although  the  "2-4"  and  "3-4"  curves  are  valid  over  the  same  Interval 
of  strain  energy,  the  two  curves  have  different  shapes  In  the  Interval 
because  different  combinations  of  data  points  are  used  In  determining  B 
and  C.  The  "2-4"  curve  must  pass  through  point  2 whereas  the  "3-4"  curve 
must  pass  through  point  3.  (Both  curves  must  pass  through  points  1 and 
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4.)  Because  the  mechanical  property  at  point  3 is  less  than  the  mechani- 
cal property  at  point  2 and  because  U3  is  greater  than  U2 , the  "3-4" 
curve  is  steeper  than  the  "2-4"  curve  for  U < U2  and  flatter  for  U > U^. 

In  the  interval  U2  < U < the  two  curves  have  basically  the  same  shape. 

As  a result,  the  "2-4"  curve  is  a better  approximation  in  the  interval 

0 < U < Ug  and  the  "3-4"  curve  is  a better  approximation  in  the  interval 
U3  < U < U^.  Obviously,  In  the  interval  U2  < U < U^,  both  approximations 
are  about  equally  good;  the  "2-4"  approximation  is  better  for  the  strain 

energies  nearer  U2,  and  the  "3-4"  approximation  is  bettet  for  the  strain 

energies  nearer  U^.  The  interval  U3  < U < in  Fig.  2-5  is  larger  than 
the  interval  0 < U ' U2  so  the  "3-4"  approximation  is  better  than  the 
"2-4"  approximation  over  a large  portion  of  the  overall  interval  of  valid- 
ity 0 < U < U^.  Thus,  the  interior  data  point  used  to  determine  B and 
C should  be  chosen  so  that  the  approximate  mechanical  property  versus  U 
curve  takes  on  the  shape  desired  by  the  user  of  the  material  model  and 
represents  the  actual  data  accurately  over  the  desired  subinterval  of 
the  overall  range  of  validity  of  the  approximation. 

With  the  criteria  discussed  above  as  a guide,  the  actual  determina- 
tion of  the  constants  A,  B,  and  C seems  quite  easy.  However,  measured 
stress-strain  data  and  corresponding  mechanical  property  - strain  energy 
data  are  somewhat  random  by  nature  and  do  not  usually  plot  as  smooth 
curves.  When  A,  B,  and  C are  determined  without  regard  for  the  random 
nature  of  the  data,  several  pitfalls  arise  which  are  not  obvious  at  first. 
These  problems  can  occur  if  a single  set  of  constants  (A,  B,  and  C)  for 
a single  mechanical  property  is  calculated  by  hand  or  if  many  sets  of 
constants  for  many  mechanical  properties  are  calculated  with  the  aid  of 
a computer. 
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The  procedure  for  determining  A,  B,  and  C can  be  quite  tedious  when 
performed  by  hand  calculation.  The  measured  stress-strain  data  must 
first  be  converted  to  mechanical  property  - strain  energy  data.  Then, 
a value  for  the  constant  A must  be  determined.  By  definition,  A is  the 
initial  elastic  value  of  the  mechanical  property,  but  obviously  this  ini- 
tial value  cannot  be  determined  from  the  data  at  the  origin  of  the  stress- 
strain  curve  (o/e  = 0).  Hence,  the  initial  slope  of  the  measured  stress- 
strain  curve,  i.e.,  the  tangent  modulus  at  the  origin  of  the  stress- 
strain  curve,  is  often  used  as  the  value  of  A.  An  alternative  approach 
for  finding  A Is  to  calculate  the  value  of  the  mechanical  property  at 
the  data  point  corresponding  to  the  lowest  measurable  stress-strain  level 
and  then  to  arbitrarily  use  this  value  of  the  mechanical  property  as  the 
value  of  A.  Thus,  in  this  approach,  a point  on  the  mechanical  property 
versus  U curve  for  which  the  strain  energy  is  small  Is  translated  to  the 
mechanical  property  axis.  One  consequence  of  determining  A in  this  manner 
Is  that  the  value  of  A is  less  than  the  value  obtained  by  using  the  ini- 
tial slope  of  the  stress-strain  curve.  Another  consequence  is  that  data 
at  this  point  cannot  be  used  in  the  calculation  of  the  constants  B and 
C.  After  A is  found,  the  actual  mechanical  property  - strain  energy 
data  is  plotted,  and  points  to  be  used  in  determining  B and  C .ire  chosen. 
Once  B and  C.  are  calculated,  the  approximate  mechanical  property  versus 
U is  plotted  to  ensure  that  a reasonable  representation  of  the  actual 
data  is  obtained. 

For  materials  with  more  than  one  stress-strain  nonlinearity,  the  use 
of  hand  calculations  in  determining  the  constants  for  each  mechanical 
property  Is  inefficient  becan.e  of  the  large  amount  of  time  involved. 

As  a result,  the  procedure  for  determining  A,  B,  and  C should  be  auto- 
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mated  so  that  the  calculations  are  performed  by  a computer  and  so  that 
the  curves  are  plotted  by  an  associated  mechanical  plotter.  In  the 
JNMDATA  computer  program  written  for  this  purpose,  the  measured  stresses 
and  strains  are  input  data.  The  corresponding  mechanical  properties  and 
strain  energies  are  calculated  in  the  program,  and  the  value  of  the  me- 
chanical property  at  the  first  input  data  point  is  used  as  the  value  of 
A.  The  user  of  the  program  is  able  to  specify  which  data  points,  other 
than  the  first  one,  are  used  in  determining  the  constants  B and  C.  After 
the  constants  B and  C are  calculated,  the  program  is  designed  so  that 
the  actual  mechanical  property  - strain  energy  data  and  the  approximate 
mechanical  property  versus  strain  energy  curve  are  plotted  on  the  same 
page.  Also,  the  measured  stress-strain  data  and  the  stress-strain  curve 
implied  from  Eq,  (2.3)  and  calculated  in  the  program  are  plotted  In  a 
similar  fashion.  Checks  of  the  input  data  are  performed  In  the  program 
as  a precautionary  measure  to  avoid  the  pitfalls  which  arise  In  choosing 
data  points  to  determine  B and  C.  These  pitfalls,  due  mainly  to  the 
random  nature  of  the  measured  stress-strain  data,  are  mentioned  briefly 
in  a preceding  paragraph  and  are  discussed  in  detail  in  the  following 
paragraphs. 

One  of  the  pitfalls  which  arises  if  B and  C are  determined  from  an 
arbitrarily  chosen  combination  of  data  points  is  that  the  value  of  the 
constant  C can  be  indeterminant.  If  B and  C are  chosen  from  data  points 
2 and  3 in  Figure  2-6a,  the  mechanical  property  at  point  2,  (M.P.)2,  is 
greater  than  A,  and  the  mechanical  property  at  point  3,  (M.P.)^,  is  less 
than  A.  As  a result,  the  argument  of  the  logarithm  in  Eq.  (2.8)  Is  nega- 
tive so  C cannot  be  determined.  If  a mechanical  property  versus  strain 
energy  curve  passed  through  points  2 and  3 in  Fig.  2-6a,  the  correspond- 
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ing  stress-strain  curve,  also  shown  in  Fig.  2-6a,  would  be  initially  con- 
cave upward  and  then  would  become  concave  downward.  This  type  of  stress- 
strain  behavior  is  highly  unusual  and  will  probably  never  be  encountered. 

However,  measured  stress-strain  data  often  has  the  characteristic  that 
at  least  one  data  point  is  out  of  line  with  the  rest  of  the  data.  If 
this  data  point  is  one  of  the  points  used  to  determine  B and  C,  the  sit- 
uation shown  in  Fig.  2-6a  where  (M.P.)j,  is  greater  than  A and  (M. P . is  ; 

less  than  A is  encountered.  When  this  situation  arises,  the  process  of 
finding  the  values  of  B and  C,  whether  performed  by  hand  or  with  the  aid 
of  a 'Computer,  should  be  stopped  and  new  data  points  chosen.  j 

The  concave  downward  then  concave  upward  stress-strain  behavior  in  ] 

Fig.  2-6b,  like  the  stress-strain  behavior  in  Fig.  2-6a,  will  probably  j 

never  be  encountered.  However,  a point  at  which  the  mechanical  property  | 

is  less  than  A,  such  as  point  2 in  Fig.  2-6b,  and  a point  at  which  the  j 

mechanical  property  Is  greater  than  A,  such  as  point  3 in  Fig.  2-6b,  j 

could  easily  be  chosen  as  the  two  data  points  to  determine  B and  C,  If  \ 

two  such  points  are  used,  the  argument  of  the  logarithm  in  Eq.  (2.8)  is 
again  negative.  The  constant  C is  again  indeterminant  so  the  process  of 
calculating  B and  C should  be  stopped,  and  new  data  points  should  be  cho- 
sen. j 

} 

The  initially  linear,  then  nonlinear,  and  finally  linear  stress-  j 

strain  behavior  in  Fig.  2-6c  is  also  highly  unusual  but  is  shown  to  11- 

i 

i 

lustrate  a problem  which  occurs  when  B and  C are  determined  from  another  1 

combination  of  data  points.  The  mechanical  property  at  point  2 in  Fig.  j 

2-6c  is  not  equal  to  A,  but  the  mechanical  property  at  point  3 is  equal  ,j 

'] 

to  A.  When  these  data  points  are  used  in  finding  the  values  of  B and  C,  j 

the  denominator  of  the  argument  of  the  logarithm  in  Eq.  (2.3)  is  zero  so 
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C cannot  be  calculated.  Thus,  this  combination  of  points  must  be  avoided 
in  choosing  data  points  to  use  in  determining  B and  C. 

A different  type  of  problem  arises  if  the  mechanical  property  at 
point  2 is  equal  to  the  mechanical  property  at  point  3 but  is  not  equal 
to  A as  in  Fig.  2-6d.  Because  of  random  measured  data,  this  situation 
can  occur  when  points  2 and  3 are  chosen  close  to  each  other.  For  this 
combination  of  data  points,  the  constant  C is  zero  so  the  approximate 
mechanical  property  is  not  dependent  on  the  strain  energy.  The  value  of 
B is  A - (M.P.)?.  The  approximate  mechanical  property  is  a constant  but 
is  neither  equal  to  A nor  to  the  mechanical  property  at  point  2 or  point 
3.  That  is,  the  stress-strain  behavior  in  Fig.  2-6d  is  approximated  by 
a straight  line  with  slope  less  than  A but  greater  than  the  value  of  the 
mechanical  property  at  point  2 or  point  3.  Thus,  the  condition  that  the 
approximate  mechanical  property  curve  must  pass  through  the  three  data 
points  used  to  determine  A,  B,  and  C is  violated.  The  unusual  stress- 
strain  curve  with  two  linear  portions  of  different  slope  in  Fig.  2~6d 
cannot  be  represented  exactly.  Only  an  approximation  of  the  behavior,  in 
which  the  two  linear  portions  are  modeled  as  being  nonlinear,  can  be  ob- 
tained when  Eq.  (2.3)  is  used  to  approximate  the  mechanical  property. 

To  obtain  this  mechanical  property  approximation,  point  2 in  Fig.  2-6d 
must  be  located  on  the  nonlinear  portion  of  the  stress-strain  curve  so 
that  values  for  B and  C can  be  found. 

Pitfalls  similar  to  those  discussed  above  can  arise  even  when  common 
types  of  stress-strain  behavior  are  modeled  if  the  points  used  to  deter- 
mine B and  C are  chosen  indiscriminately.  The  stress-strain  behavior  of 
many  materials  is  represented  by  the  initially  linear  and  then  nonlinear 
stress-strain  curve  in  Fig.  2 - 7 a . However,  Eq.  (2.3)  cannot  be  used  to 
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FIGURE  2-7  MECHANICAL  PROPERTY  VERSUS  U BEHAVIORS  WHICH 
CAUSE  DIFFICULTIES  IN  DETERMINING  B AND  C 
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approximate  the  corresponding  mechanical  property  behavior  if  points  2 
and  3 in  Fig.  2-7a  are  used  in  determining  B and  C.  The  mechanical  prop- 
erty at  point  2 is  equal  to  A but  not  equal  to  the  mechanical  property 
at  point  3.  The  constant  C is  indetermi nant  hecause  the  argument  of 
the  logarithm  in  Eq.  (2.8)  is  zero.  Hence,  the  procedure  for  calculating 
B and  C must  be  stopped.  This  problem  is  avoided  by  choosing  point  2‘  on 
the  nonlinear  portion  of  the  stress-strain  curve  to  be  used  to  determine 
B and  C instead  of  point  2.  The  implied  stress-strain  curve  resulting 
from  this  approximation  is  nonlinear  at  all  levels  of  stress  and  strain 
and  is  a reasonable  approximation  of  the  stress-strain  curve  in  Fig. 

2-7a,  but  the  measured  linear  then  nonlinear  stress-strain  behavior  can- 
not be  modeled  exactly  with  Eq.  (2.3). 

Even  for  linear  stress-strain  behavior  shown  in  Fig.  2-7 b f B and  C 
must  be  determined  carefully.  Obviously,  the  value  of  the  mechanical 
property  at  points  2 and  3 in  Fig.  2-7b  is  equal  to  A,  but  If  B and  C 
are  determined  with  Eq.  (2.8)  and  Eq.  (2.9),  the  argument  of  the  loga- 
rithm in  Eq.  (2.8)  is  indeterminant  (0/0)  so  C cannot  be  calculated. 
However,  linear  stress-strain  behavior  is  implied  from  Eq.  (2.3)  if  B 
is  zero  and  C is  finite.  Thus,  the  problem  of  C being  indeterminant 
Is  solved  by  skipping  the  calculations  In  Eqs,  (2.8)  and  (2.9)  and  assign 
ing  B the  value  of  zero  and  C some  arbitrary  but  finite  value. 

With  B and  C determined  from  a proper  choice  of  data  points  on  the 
mechanical  property  versus  strain  energy  curve,  the  mechanical  property 
equation  is  a vaMd  approximation  of  the  actual  mechanical  property  data 
over  a specfied  range  of  strain  energy.  However,  under  many  loading 
conditions,  a value  of  the  mechanical  property  is  desired  for  a strain 
energy  larger  than  the  maximum  value  of  U for  which  the  actual  mechani- 
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FIGURE  2-8  UNIAXIAL.  STRESS-STRAIN  BEHAVIOR  AND  CORRESPONDING 
MATERIAL  PROPERTY  VERSUS  U BEHAVIOR 

cal  property  versus  U curve  Is  defined.  For  example,  the  mechanical 
property  versus  U curve  in  Fig.  2-8  is  defined  for  strain  energies  in 
the  range  0 < U < U This  curve  is  obtained  from  the  uniaxial  stress 
strain  curve  also  in  Fig.  2-8  which  is  not  defined  for  stresses  greater 
than  and  strains  greater  than  Hence,  „ is  not  defined  for 

strain  energies  greater  than  Um„u  B e„.au/ 2.  A problem  arises  if  a 

ITlo.X  ITlaX  maX 

value  of  the  mechanical  property  Is  needed  for  a strain  energy  greater 
than  U . Thus,  the  actual  mechanical  property  versus  strain  energy 

max 

curve  must  be  extrapolated  in  some  way.  Two  extrapolation  procedures 
for  the  material  model  are  described  iri  Section  2.3. 

2.2.3  TEMPERATURE  INTERPOLATION  OF  DEFORMATION  BEHAVIOR 

Many  materials  have  temperature-dependent  deformation  behavior. 
Typically,  that  deformation  behavior  is  measured  at  discrete  values  of 
temperature  over  some  pertinent  range  of  temperatures.  We  propose  to 
use  such  data  in  finite  element  stress  analysis  of  bodies  subjected  to 
temperature  gradients.  Then,  the  basic  problem  is:  given  two  stress- 
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strain  curves  at  two  temperatures  with  associated  mechanical  property 


versus  strain  energy  curves  for  which  we  known  A,  B,  C,  and  U0,  find 
the  correct  property  for  an  intermediate  temperature  at  a specified  ener- 
gy level  along  with  the  associated  point  on  the  stress-strain  curve  for 
the  intermediate  temperature.  Repeated  answering  of  this  question  for 
: any  energy  levels  would  lead  us  to  a mechanical  property  versus  energy 
curve  and  a stress-strain  curve  for  the  intermediate  temperature.  These 
curves  must  lie  appropriately  between  the  respective  two  original  curves 
for  the  temperatures  between  which  we  interpolated  the  deformation  be- 
havior. 

At  least  three  approaches  are  possible:  parameter  interpolation, 

property  Interpolation,  and  stress-strain  curve  interpolation.  Each  of 
these  approaches  is  defined  and  discussed  In  this  section.  Moreover, 
the  consequences  of  using  each  approach  are  discussed  relative  to  what 
is  perceived  to  be  the  desired  result.  The  word  perceived  Is  used 
because  the  actual  variation  of  mechanical  properties  with  temperature 
is  often  not  known  except,  as  stated  before,  at  discrete  temperatures. 
That  1$,  the  actual,  smooth  variation  with  temperature  of  properties  or 
stress-strain  curves  is  not  known  so  we  have  little  with  which  to  compare 
our  results. 

The  basic  example  for  all  three  approaches  will  Involve  the  inter- 
polation of  ATJ-S  graphite  properties  at  1403°F  when  we  are  given  the 
properties  at  70°F  and  2000°F.  We  examine  only  the  deformation  behavior 
in  the  Isotropic  plane,  i.e.,  o0  versus  e0  and  E0  versus  U (or  or  versus 

e arid  E_  versus  U).  The  stress-strain  curves  are  shown  for  70°F  and 
r r 

2000°F  in  Fig.  2-9  and  are  taken  from  an  SoRI  report  [2-8].  The  boxes 
in  Fig.  2-9  are  the  SoRI  data  points  and  the  curves  are  those  obtained 
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FIGURE  2-10  ATJ-S  GRAPHITE  MECHANICAL  PROPERTY  (E„  - 
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VERSUS  STRAIN  ENERGY  FOR  70°F  AND  2000°F 
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with  the  JNMDATA  program  which  will  be  described  in  Section  2.4.  The 
associated  mechanical  property  versus  strain  energy  curvts  are  displayed 
in  Fig.  2-10  for  which  the  mechanical  property  constants  in  Eq.  (2.3) 
are  given  in  Table  2-1 . 

TABLE  2-1  ATJ-S  GRAPHITE  MECHANICAL  PROPERTY  CONSTANTS 

FOR  E„  = E„  VERSUS  ENERGY 
r u 


T 

A, psi 

B 

C 

U0,psi 

70°F 

2,000,000. 

.182 

.337 

1. 

2000° F 

2,000,000. 

.0651 

.583 

1. 

2.2.3. 1 Parameter  Interpolation 

In  parameter  interpolation,  the  values  of  the  parameters  (mechani- 
cal property  constants)  A,  B,  C,  and  UQ  for  the  intermediate  temperature 
T are  found  by  linear  Interpolation  of  the  corresponding  parameters  at 
the  bounding  temperatures  T^  and  T^.  For  example,  if  T B (T1  + T2 )/2 , 
then 

A - (A1  + A2)/2  B ■ (B]  + B2)/Z 

C ■ (C,  ♦ C2)/2  U0  . (U0i  * 1^1/2 

We  might  expect  the  resulting  mechanical  property  versus  energy  and 
stress-strain  curves  to  lie  midway  between  the  respective  interpolated 
curves.  However,  they  do  not.  In  fact,  at  a temperature  of  1403°F  for 
which  A = 2,000,000  psi,  B » .101,  C = .507,  and  UQ  ■ 1 psi  (from  linear 
interpolation),  at  stresses  above  7,000  psi  in  Fig.  2-11  the  resulting 
stress-strain  curve  lies  below  the  70°F  curve  and  not  fairly  close  to 
the  2000°F  curve  as  we  would  hope.  Moreover,  at  U = 40  psi,  the  result- 
ing E„  Is  less  than  E„  for  70°F.  Part  of  this  undesirable  result  may 
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stem  from  the  simultaneously  constant  A's,  decreasing  B's,  and  increasing 
C's  (but  not  from  the  constant  u ).  At  any  rate,  parameter  interpolation 
is  obviously  undependable  since  the  resulting  interpolated  curves  may  not 
even  lie  between  the  curves  being  Interpolated! 

2. 2. 3.2  Property  Interpolation 

In  property  interpolation,  the  mechanical  property  for  the  inter- 
mediate temperature  T is  found  by  linear  interpolation  of  the  mechanical 
properties  at  the  bounding  temperatures  and  Tg  for  a specified  ener- 
gy level.  That  is,  if  for  example  T lies  midway  between  and  T-,  then 

Mechanical  PropertyT  = ij-  (Mechanical  Property^  + Mechanical  Property^  ) 

(2.10) 

The  resulting  mechanical  property  versus  energy  curve  lies  precisely  mid- 
way between  the  two  curves  being  Interpolated,  by  definition.  However, 
the  resulting  stress-strain  curve  may  be  closer  to  the  higher  of  the  two 
curves  being  interpolated.  When  T a 1403°F,  the  Interpolated  stress- 
strain  curve  in  Fig.  2-13  is  "about  the  right  distance"  from  the  two 
curves  from  which  It  is  indirectly  Interpolated.  Also,  the  Interpolated 
mechanical  property  versus  energy  curve  is  precisely  where  it  must  be 
when  we  assume  that  the  mechanical  properties  are  a linear  function  of 
temperature  between  the  discrete  temperature  values  at  which  the  mechani- 
cal behavior  is  measured.  We  delay  an  appraisal  of  the  value  of  this 
method  until  after  the  next  sub-section  on  stress-strain  curve  Interpo- 
lation. 
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FIGURE  2-14  MECHANICAL  PROPERTY  VERSUS  ENERUV  CURVE  AT  1403°F 
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2. 3. 3. 3 STRESS-STRAIN  CURVE  INTERPOLATION 


In  stress-strain  curve  interpolation,  the  coordinates  of  points  on 
the  interpolated  stress-strain  curve  are  found  from  the  information 
known  for  the  stress-strain  curves  at  the  bounding  temperatures.  Actual- 
ly, we  are  assuming  that  the  stress-strain  curves  vary  linearly  with 
temperature  between  the  curves  at  input  discrete  temperature  levels. 

For  example,  the  stress-strain  curve  for  a temperature  midway  between 
the  two  temperatures  at  which  the  stress-strain  curves  are  prescribed 
lies,  by  definition,  precisely  midway  between  the  two  prescribed  stress- 
strain  curves.  However,  since  there  is  no  direct  equation  for  the 
stress-strain  curves,  we  must  construct  the  Interpolated  stress-strain 
curve  by  calculation  of  stresses  and  strains  for  a succession  of  energy 
levels. 

The  calculation  procedure  is  straightforward,  but  several  steps  are 
necessary  for  each  energy  level: 

(1)  Calculate  E1  and  E2  at  T1  and  T^,  the  bounding  temperatures 
(the  1 and  2 subscripts  will  be  used  hereafter  in  association 
with  T'|  and  T2,  respectively). 

(2)  Calculate  the  corresponding  strains  and  c2  by  realizing  that 

a1  " E1g1 

a2  = E2e2 
but 

U - . \ 

u “ Jr  02f" 2 “ 7 ^2e2^e2 
so 

- (2U/E, )1/2 
e2  " (2U/E2)1/2 
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(3)  Calculate  the  corresponding  stresses  and  Og  from 

o.,  - E -j e-j  = (2UE-,  )1/2 
ci 2 - E 2r2  * (2UE2)1/2 

(4)  Interpolate  (a-j  .c.j ) and  (og.r^)  along  a straight  line  joining 
them  to  get  (oy,eT): 


°T 


a^)  + o 


1 


AT 


CT  “ T2  - ^ *e2  " E1  ^ + C1 


(5)  Calculate  the  mechanical  property  at  T: 


Et  ■ — 

T tT 


The  foregoing  five  steps  are  repeated  for  increasing  energy  levels  until 
the  stress-strain  and  mechanical  property  versus  energy  curves  are  suf- 
ficiently well  defined. 

With  this  approach,  the  Oy  - ey  curve  lies  midway  between  the  two 
stress-strain  curves  if  T « (T^  + Tg)/2.  Also,  the  interpolated  mechani- 
cal property  versus  energy  lies  closer  to  the  lower  specified  curve  than 
to  the  upper  curve.  However,  for  the  ATJ-S  graphite  example  at  1403°F, 
the  resulting  stress-strain  curve  is  indistinguishable  from  that  obtained 
with  property  interpolation  in  Fig.  2-13  (slight  numerical  differences 
exist,  but  they  are  t visually  detectable  at  the  scale  of  Fig.  2-13). 
Moreover,  the  resulting  mechanical  property  versus  energy  curve  is  indis- 
tinguishable from  that  obtained  with  property  interpolation  in  Fig.  2-14. 


> 


2, 2. 3. 4 Summary 

Three  different  approaches  to  interpolation  of  deformation  behavior 
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are  discussed.  One,  parameter  interpolation,  is  found  to  be  inaccurate 
and  unreliable.  The  other  two,  property  interpolation  and  stress-strain 
curve  interpolation,  are  found  to  be  equivalent.  The  two  methods  are 
both  applied  to  the  case  where  the  temperature  is  midway  between  the  two 
input  temperatures.  Although  the  interpolated  mechanical  property  versus 
energy  curve  is  precisely  midway  between  the  two  input  curves  for  proper- 
ty interpolation,  the  interpolated  stress-strain  curve  might  theoretical - 
ly  be  a little  too  high.  On  the  other  hand,  the  interpolated  stress- 
strain  curve  is  precisely  midway  between  the  two  input  curves  for  stress- 
strain  curve  interpolation,  but  the  interpolated  mechanical  property  ver- 
sus energy  curve  might  theoretically  be  a little  too  low.  The  two  methods 
are  also  applied  to  an  actual  case  for  ATJ-S  graphite  where  the  tempera- 
ture is  1403°F  and  stress-strain  data  are  known  at  70°F  and  2000°fr.  In 
this  practical  case,  the  differences  between  the  two  Interpolation  ap- 
proaches are  negligible.  That  is,  property  interpolation  leads  to  an 
interpolated  stress-strain  curve  which  is  also  linearly  interpolated 
between  the  respective  input  stress-strain  curves.  And,  stress-strain 
curve  Interpolation  leads  to  an  interpolated  mechanical  property  versus 
energy  curve  which  is  also  linearly  interpolated  between  the  respective 
input  mechanical  property  versus  energy  curves.  Thus,  the  two  approaches 
are  equivalent  for  practical  purposes.  The  property  interpolation  ap- 
proach is  used  in  the  SAAS  HIM  program  and  the  steps  involved  In  proper- 
ty interpolation  are  simpler  than  those  for  stress-strain  curve  interpola- 
tion. 
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2.3  EXTRAPOLATION  PROCEDURES  FOR  MATERIAL  MODELS 
2.3.1  INTRODUCTION 

The  strain  energy  can  exceed  the  defined  range  of  the  mechanical 
property  versus  U curve  for  two  reasons.  First,  the  nonlinear  model  is 
applied  to  materials  subjected  to  multiaxial  stress  states  in  which  the 
strain  energy  is  higher  than  in  the  uniaxial  stress  states  where  the 
properties  are  measured.  All  the  stresses  and  strains  of  a multiaxial 
stress  state  contribute  to  the  value  of  1),  and,  thus,  the  multiaxial 
strain  energy  is  often  larger  than  the  maximum  U attainable  from  a 
uniaxial  stress  state.  The  second  reason  for  the  existence  of  strain 
energies  outside  the  defined  range  of  the  mechanical  property  versus  U 
curve  is  that  orthotropic  materials  have  drastically  different  load  capa- 
cities and  hence  drastically  different  strain  energy  capacities  in  dif- 
ferent directions.  The  strain  energy  capacities  for  loading  in  the  2- 
direction  and  for  shear  loading  are  generally  much  lower  than  the  strain 
energy  capacity  for  loading  in  the  1 -direction.  For  example,  the  repre- 
sentative longitudinal  stress  - longitudinal  strain  curve  in  Fig.  2-15 
Is  associated  with  strain  energies  as  high  as  250  psi.  However,  the 
maximum  value  of  U which  can  be  used  to  define  the  mechanical  property 

E£  from  the  transverse  stress  - transverse  strain  curve  in  Fig.  2-8 
sec 

is  only  40  psi.  Thus,  strain  energies  corresponding  to  some  uniaxial 
loading  conditions,  such  as  loading  in  the  1 -direction  or  off-axis  load- 
ing, can  be  much  larger  than  the  Um_w  for  which  E9  is  defined. 

max  Sec 

The  mechanical  property  versus  U curve  can  be  extended  past  its  de- 
fined range  by  the  two  extrapolation  procedures  described  in  this  section. 
The  procedure  described  In  Section  2.3.2  consists  of  using  Eq.  (2.3)  for 
all  values  of  strain  energy  including  those  which  exceed  the  maximum 
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FIGURE  2-15  NORMAL  STRESS  - NORMAL  STRAIN  BEHAVIOR 
OF  AN  ORTHOTROPIC  MATERIAL 


strain  energy  for  which  the  actual  mechanical  property  data  are  defined. 
When  the  extrapolation  procedure  in  Section  2.3.2  fails,  the  extrapola- 
tion procedures  described  in  Section  2.3.3  are  used.  There,  extensions 
of  the  defined  stress-strain  curve  are  involved  as  opposed  to  the  direct 
extension  of  the  mechanical  property  versus  strain  energy  curve  In  Sec- 
tion 2.3.2.  The  stress-strain  curve  is  extended  by  following  the  trend 
of  the  defined  portion  of  the  curve  until  the  slope  attains  some  preas- 
signed value,  i.e.,  until  ^ • For  strains  larger  than  e.*,  the 

strain  at  which  pp  = . the  stress-strain  curve  is  extended  as  a 

straight  line  with  slope  f J to  represent  the  behavior  of  ductile  fi- 
ber-reinforced composite  materials.  Equation  (2.3)  is  used  to  represent 
the  mechanical  property  - U behavior  for  strain  energies  corresponding 
to  points  on  the  stress-strain  curve  to  the  left  of  e*.  Another  expres- 
sion which  will  be  developed  in  Section  2.3.3  is  used  to  define  the  me- 
chanical property  for  strains  larger  than  c*. 
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2,3,2  EXTENDED  MECHANICAL  PROPERTY  VERSUS  STRAIN  ENERGY  CURVE  APPROACH 


The  simplest  means  of  extrapolating  the  mechanical  property  versus 
strain  energy  curve  is  to  use  the  mechanical  property  equation  for  all 
values  of  U including,  because  of  extrapolation,  those  beyond  the  defined 
range  of  strain  energy.  A representative  mechanical  property  \ v-sus  U 
curve  is  shown  in  Fig.  2-16.  This  curve  is  defined  for  strain  energies 
less  than  or  equal  to  Umax.  Also  shown  in  Fig.  2-16  is  an  extended  ap- 
proximate mechanical  property  versus  U curve  corresponding  to  Eq.  (2.3). 
For  this  curve,  B and  C arc  determined  from  the  data  point  with  strain 
energy  equal  to  Umax  and  another  point  with  strain  energy  close  to  Umax. 
This  approximate  curve  is  an  accurate  interpolation  of  the  actual  data 
for  strain  energies  between  zero  and  Umax  as  discussed  earlier.  In  addi- 
tion, this  curve  has  the  same  shape  as  the  actual  mechanical  property 
versus  U curve  in  the  neighborhood  of  Umax  and  is  a reasonable  extension 
of  the  actual  data  for  strain  energies  greater  than  II  . Thus,  if  data 
at  the  point  corresponding  to  LI  = Umax  and  data  at  another  point  near  the 


MECHANICAL 

PROPERTY 


FIGURE  2-16  ACTUAL  AND  EXTRAPOLATED  MECHANICAL  PROPERTY  VERSUS  U BEHAVIOR 
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limit  of  the  defined  mechanical  property  versus  U curve  are  used  to  deter- 
mine the  constants  B and  C,  the  use  of  Eq.  (2.3)  is  potentially  a valid 
means  of  extrapolating  the  actual  mechanical  property  - U data. 

However,  the  use  of  Eq.  (2.3)  as  an  extrapolation  is  restricted, 

At  some  large  value  of  U,  the  approximate  mechanical  property  curve  in 
Fig.  2-16  crosses  the  U-axis,  i.e.,  the  extrapolated  value  of  the  mechan- 
ical property  becomes  negative.  The  strain  energy  at  which  the  mechani- 
cal property  is  zero  is  designated  as  0’  in  Fig.  2-16.  If  a mechanical 
property  is  less  than  or  equal  to  zero,  a thermodynamic  constraint  on  the 
mechanical  properties  of  a material  is  violated.  This  constraint,  imposed 
on  the  properties  to  avoid  the  creation  of  energy,  is  that  the  work  done 
by  the  stresses  applied  to  a material  must  be  positive.  Lampriere  [2-9] 
interprets  this  constraint  as  meaning  that  both  the  stiffness  and  com- 
pliance matrices  of  an  orthotropic  material  must  be  positive  definite. 
These  matrices  are  positive  definite  only  if  the  mechanical  properties 
ly  E2,  and  G-|2  are  positive.  Thus,  the  mechanical  properties  of  a mate- 
rial cannot  be  defined  with  Eq.  (2.3)  for  strain  energies  greater  than 
or  equal  to  U.  The  value  of  U can  be  determined  from  Eq.  (2.3)  and  de- 
pends only  on  the  constants  B and  C: 


1 

U.(I)C  (2.11) 

For  a constant  value  of  B,  the  value  of  "U  increases  when  C decreases. 

The  stress-strain  behavior  implied  from  Eq.  (2.3)  must  be  a reason- 
able extension  of  the  actual  stress-strain  behavior  of  the  material.  For 
the  implied  stress-strain  curves  in  Fig.  2-17,  the  stresses  reach  a maxi- 
mum value  and  then  decrease  with  Increasing  strain.  As  the  strain  in- 
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creases  without  bound,  the  stress  approaches  zero  asymptotically.  Also, 
the  strain  energies  determined  from  the  implied  stress-strain  curves  in 
Fig.  2-17  approach  but  never  reach  the  value  O’,  i.e., 


Lim  U 
a -►  o 

e -*  « 


(2.12) 


Hence,  O'  is  the  maximum  strain  energy  which  can  be  associated  with  both 
the  mechanical  property  equation  and  its  implied  stress-strain  curve. 

The  stress-strain  curve  implied  from  the  mechanical  property  equa- 
tion is  not  necessarily  representative  of  the  actual  stress-strain  be- 
havior of  fiber-reinforced  materials  in  the  range  0 < U < 0".  The  schema- 
tic implied  stress-strain  curve  shown  In  Fig.  2-18  reaches  a maximum  at 
the  point  (a,  ?.).  The  portion  of  the  curve  to  the  left  of  (d,  if)  Is 
representative  of  actual  stress-strain  behavior.  The  implied  stress- 
strain  behavior  to  the  right  of  (a,  c)  where  the  stress  decreases  as  the 
strain  Increases  has  not  been  observed  for  fiber-reinforced  composite 


FIGURE  2-18  REPRESENTATIVE  IMPLIED  STRESS-STRAIN  BEHAVIOR 
CORRESPONDING  TO  JONES-NELSON  EQUATION 
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materials.  Thus,  the  extended  mechanical  property  versus  strain  energy 
curve  approach  cannot  be  used  as  an  extrapolation  for  strain  energies  as 
large  as  U,  but  must  be  restricted  to  strain  energies  less  than  0 where 

0 ■ 1/2  or  (U  < U)  (2.13) 

2,3.3  EXTENDED  STRESS-STRAIN  CURVE  APPROACH 

An  extension  of  the  mechanical  property  versus  strain  energy  curve 
by  the  use  of  the  mechanical  property  equation  is  indirectly  an  exten- 
sion of  the  actual  stress-strain  curve;  i.e.,  the  measured  stress-strain 
curve  is  extended  past  its  defined  range  with  an  implied  stress-strain 
curve  similar  to  the  one  in  Fig.  2-18.  However,  for  strain  energies 
greater  than  0,  the  stress-strain  behavior  must  be  extrapolated  with  some 
procedure  other  than  the  extended  mechanical  property  versus  strain  ener- 
gy curve  approach.  A reasonable  extrapolation  is  to  extend  the  measured 
stress-strain  curve  directly,  then  convert  the  extrapolated  stress-strain 
data  into  extrapolated  mechanical  property  versus  strain  energy  data, 
and  finally  develop  an  expression  for  the  extrapolated  mechanical  property 
data  In  terms  of  the  strain  energy.  Several  extensions  of  the  stress- 
strain  data  are  discussed  in  this  section.  Linear  extensions  of  the 
stress-strain  behavior  with  zero  slope  and  with  nonzero  slope  are  dis- 
cussed in  Section  2. 3. 3,1  and  2. 3. 3. 2,  respectively. 

2. 3. 3.1  Linear  Stress-Strain  Curve  Extensions  wi Lh  Zero  Slope 

One  linear  extension  of  the  defined  stress-strain  curve  involves 
estimating  the  stress  level  at  which  the  nonlinear  curve  has  zero  slope. 
The  stress-strain  curve  In  Fig.  2-1 9a  is  defined  up  to  point  P.  The 
trend  of  the  defined  stress-strain  curve  is  followed,  and  the  curve  is 
arbitrarily  extended  until  the  slope  becomes  zero  at  point  Q,  The  : tress 
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(b)  MECHANICAL.  PROPERTY  VERSUS  STRAIN  ENERGY  CURVE  EXTRAPOLATION 

FIGURE  2-19  LINEAR  STRESS-STRAIN  CURVE  EXTRAPOLATION  WITH  ZERO  SLOPE 
BY  ARBITRARY  EXTENSION  OF  STRESS-STRAIN  DATA 


and  strain  associated  with  point  Q are  a*  and  e*,  respectively.  For 
strain  energies  less  than  e*  (U  < U*  = 1/2  o*e*),  the  mechanical  proper- 
ty is  expressed  as  a function  of  the  strain  energy  with  the  Jones-Nelson 
mechanical  property  expression  in  Eq.  (2.3).  The  constants  B and  C are 
determined  from  data  at  points  P and  Q.  For  strains  greater  than  e*» 
the  stress-strain  curve  is  considered  to  be  linear  with  zero  slope,  and 
the  mechanical  property  and  strain  energy  can  be  expressed  as: 

Mechanical  Property  * — e > e*  (2.14) 


U » a*e/2  e > e*  (2.15) 

Equations  (2.14)  and  (2.15)  are  combined  to  obtain  an  expression  for  the 
mechanical  property  as  a function  of  II: 


Mechanical  Property 


u > u* 


(2.16) 


Thus,  when  the  stress-strain  curve  is  extended  horizontally  from  the 
point  of  zero  slope,  the  expression  for  the  mechanical  property  for  all 
strain  energies  Is: 


Mechanical  Property 


The  mechanical  property  versus  strain  energy  curve  in  Fig.  2-1 9b 
corresponding  to  Eq.  (2.17)  approaches  the  U-axis  asymptotically  so  the 
mechanical  property  is  always  positive.  Hence,  the  horizontal  stress- 
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chanical  property  curve  which  crosses  the  U-axis  as  a strain  energy  of  U. 

The  major  problem  in  extending  the  stress-strain  curve  horizontally 
in  the  foregoing  manner  is  the  determination  of  the  point  Q in  Fig.  2-1 9a . 
The  stress  a*  associated  with  point  Q is  determined  by  using  engineering 
judgment  to  estimate  the  stress  level  at  which  the  extended  stress-strain 
curve  becomes  horizontal.  Obviously,  the  exact  stress  level  at  which  the 
curve  becomes  horizontal  is  ambiguous  and  depends  on  how  the  curve  is 
extended.  Many  different  choices  of  a*  could  be  made.  After  a*  is 
chosen,  a complicated  trial  and  error  procedure  based  on  mathematical 
constraints  on  the  mechanical  property  expressions  in  Eq.  (2.17)  is  used 
to  determine  e*.  The  objective  of  this  procedure  is  to  determine  the 
strain  (e*)  at  which  the  slope  of  the  implied  stress-strain  curve  is  zero. 
A value  of  e*  is  chosen,  and  then  the  constants  B and  C are  determined 
from  data  at  point  P and  the  point  (cs*,  e*).  Associated  with  each  mechan- 
ical property  equation  is  an  implied  stress-strain  curve  which  Is  obtained 
by  substitution  of  the  definitions  of  the  mechanical  property  and  strain 
energy  in  Eq.  (2.3)  as 


2-  * A 1 


(2.18) 


For  a specific  value  of  strain,  the  corresponding  stress  is  the  root  of 
the  nonlinear  equation: 

| - a[i  - B(S“)C]r-  0 (2.19) 

The  slope  of  the  implied  stress-strain  curve,  obtained  by  differentiation 
of  Eq.  (2.18)  and  given  by 
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(2.20) 


* —< 

A (211  )C-B(C+1  ) (o*(.*)C 
l—  2 -J. 

(2Uo)C  + ABC  (o*)C"1(t:*)'C+1" 

is  then  evaluated  at  (a*,  t*).  Then,  the  strain  t*  is  adjusted  and  the 
constants  B and  C are  recalculated  until  the  slope  in  Eq.  (2.20)  is  zero. 
The  numerical  technique  of  Interval  halving  described  in  the  Appendix  is 
used  in  adjusting  t:*  so  that  the  search  for  the  point  of  zero  slope  is 
orderly  and  converges  rapidly. 

The  procedure  for  finding  point  Q Is  ambiguous  and  quite  complica- 
ted. In  addition,  this  horizontal  extension  of  the  stress-strain  curve 
may  lead  to  a poor  representation  of  the  measured  stress-strain  data. 

As  a result,  this  method  of  extending  the  stress-strain  curve  horizon- 
tally Is  used  only  when  very  little  stress-strain  data  are  available  and 
when  the  trend  of  the  stress-strain  data  Is  such  that  any  of  the  other 
extensions  to  be  discussed  subsequently  is  not  valid. 

The  mechanical  property  expression  in  Eq.  (2.17)  Is  also  valid  for 
another  linear  stress-strain  curve  extension  with  zero  slope  which  Is 
less  complicated  than  the  one  just  described.  Moreover,  this  extension 
is  more  representative  of  the  measured  stress-strain  data  for  fiber- 
reinforced  composite  materials.  This  extrapolation  Involves  adding  a 
linear  extension  with  zero  slope  to  the  implied  stress-strain  curve  which 
is  the  best  representation  of  the  defined  stress-strain  data  as  shown  in 
Fig.  2-20a.  The  constants  A,  B,  and  C in  Eq.  (2.17)  are  chosen  so  that 
the  best  approximation  of  the  mechanical  property  versus  strain  energy 
data  is  obtained.  Then,  the  implied  stress-strain  curve  corresponding 
to  this  approximation  Is  used  not  only  as  a representation  of  the  stress- 
strain  data  but  also  as  an  extension  of  the  stress-strain  data  from  the 
last  data  point  to  the  point  (u*,  e*)  in  Fig.  2-20a.  The  point  (a*,  c*) 
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(a)  STRESS-STRAIN  CURVE  EXTRAPOLATION 


MECHANICAL 

PROPERTY 


(b)  MECHANICAL  PROPERTY  VERSUS  STRAIN  ENERGY  CURVE  EXTRAPOLATION 

FIGURE  2-20  LINEAR  STRESS-STRAIN  CURVE  EXTRAPOLATION  WITH  ZERO  SLOPE 
WITH  BEST  FIT  EXTENSION  OF  STRESS-STRAIN  DATA 
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is  the  point  on  the  extended  stress-strain  curve  at  which  a change  in  the 
extrapolated  stress-strain  behavior  occurs.  For  this  extrapolation, 

(a*,  c*)  in  Fig.  2~20a  corresponds  to  the  point  (a,e)  in  Fig.  2-18;  i.e., 
the  slope  of  the  implied  stress-strain  curve  in  Fig.  2-20a  is  zero  at 
(a*,  e*).  The  implied  stress-strain  curve  is  not  a reasonable  extension 
of  fiber-reinforced  composite  material  behavior  for  strains  greater  than 
e*;  so  the  stress-strain  curve  is  extended  as  the  horizontal  line  o = o* 
in  Fig.  2-20a. 

This  method  of  extending  the  stress-strain  curve  with  a straight 
line  of  zero  slope  Is  reasonable  and  simple  because  the  constants  A,  B, 
and  C in  tq.  (2.17)  ure  determined  by  fitting  only  the  defined  data  and 
because  the  point  (a*,  c*)  is  easily  determined.  The  point  (o*.  e*)  is 
the  point  at  wnlch  in  Eq,  (2.20)  is  zero  and  is  determined  with  an 
interval  halving  technique  which  is  described  in  the  Appendix. 

The  mechanical  property  versus  strain  energy  curve  in  Mg.  2-20b 
corresponding  vo  the  zero  slope  linear  extension  of  the  implied  stress- 
strain  curve  is  similar  to  the  mechanical  property  versus  strain  energy 
curve  corresponding  to  tne  horizontal  extension  in  which  a * is  estimated 
from  the  trend  of  the  data.  Both  mechanical  property  curves  approach  the 
U-axis  asymptotically  at  large  strain  energies.  However,  the  curve  in 
Fig.  2-20b  is  a good  representation  of  the  defined  data,  whereas  the 
curve  in  Fig.  2-19  is  not. 

The  linear  stress-strain  curve  extensions  with  zero  slope  discussed 
in  this  section  are  reasonable,  but  their  use  is  limited.  The  mechanical 
property  approaches  zero  for  large  strain  energies  although  it  Is  never 
zero  nor  negative.  This  behavior  leads  to  problems  when  the  mechanical 
property  expression  is  used  in  stress  analysis  problems.  The  Iteration 
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procedure  designed  for  use  in  the  solution  of  stress  analysis  problems 
diverges  at  large  strain  energies  when  the  zero  slope  linear  extension 
of  the  stress-strain  data  is  used.  The  limit  of  validity  for  the  hori- 
zontal extensions  depends  on  the  particular  loading  condition  under  in- 
vestigation. Thus,  an  extrapolation  which  can  be  used  regardless  of  the 
value  of  the  strain  energy  is  not  obtained  with  the  linear  extension  of 
zero  slope  approaches. 

2.3.3.?  Linear  Stress-Strain  Curve  Extensi ons  with  Nonzero  Slope 

The  stress-strain  data  can  be  extrapolated  by  using  a linear  exten- 
sion with  nonzero  slope.  This  extension  can  be  used  for  all  strain  ener- 
gies and  can  be  accomplished  in  several  ways.  One  approach  is  to  approx- 
imate the  stress-strain  data  (mechanical  property  - U data)  with  the  Jones- 
Nelson  mechanical  property  expression  in  Eq.  (2.3)  and  to  arbitrarily 
specify  the  slope  that  the  linear  extension  should  have.  The  Implied 
stress-strain  curve  as  shown  in  Fig.  2-21a  Is  initially  used  as  the  ex- 
tension of  the  measured  stress-strain  data.  When  the  slope  of  the  Im- 
plied curve  reaches  the  prescribed  value,  the  stress-strain  behavior  is 
represented  by  a straight  line  with  the  prescribed  slope.  This  linear 
stress-strain  curve  extension  is  expressed  as: 

0 *(aO  E + °o  E>E*  (2-2’> 

f d V* 

where  ( J is  the  prescribed  slope  and  u is  the  stress  at  which  the 

linear  extension  Intersects  the  stress  axis.  The  constant  cr  Is  easily 

o 

determined  once  the  point  (a*,  e*)  in  Fig.  2-21a  Is  found.  The  point 
(a*,  (•:*)  at  which  the  linear  extension  begins  is  calculated  from  the 
slope  expression  in  Eq.  (2.20). 
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(a)  STRESS-STRAIN  CURVE  EXTRAPOLATION 
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*.  (b)  MECHANICAL  PROPERTY  VERSUS  STRAIN  ENERGY  CURVE  EXTRAPOLATION 
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f Y,,  FIGURE  2-21  LINEAR  STRESS-STRAIN  CURVE  EXTRAPOLATION  WITH  NONZERO  SLOPE 

f:  M 

F | WITH  BEST  FIT  OF  STRESS-STRAIN  DATA 


The  constants  A,  B,  and  C in  the  mechanical  property  expression  are 
chosen  so  that  a good  representation  of  the  measured  stress-strain  data 
is  obtained.  This  representation  of  the  data  is  not  altered  by  the  addi- 


tion of  the  linear  extension.  In  addition,  the  secant  mechanical  proper- 


ty in  Mg.  2-21  b approaches  the  slope 
large  strain  energies;  i.e.. 


of  the  linear  extension  for 
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(2.22) 


The  slope  is  nonzero  so  the  extrapolation  procedure  of  using  the 

linear  extension  with  nonzero  slope  can  be  used  for  all  strain  energies 
greater  than  U*.  Actually,  the  linear  extension  with  zero  slope  de- 
scribed in  the  preceding  section  is  analogous  to  this  linear  extension 

with(df)  = °* 

Another  approach  very  similar  to  the  one  just  described  is  to  use 
ihe  Jones-Nelson  mechanical  property  expression  to  fit  the  stress-strain 
behavior  only  ov*_r  the  defined  range.  At  the  last  data  point  in  Fig. 
2-22a,  the  stress-strain  behavior  is  extended  as  a straight  line  with 
slope  equal  to  the  slope  of  the  implied  stress-strain  curve  at  the  last 
data  point.  For  continuity  between  the  implied  stress-strain  curve  and 
the  linear  extension  of  the  stress-strain  behavior,  the  implied  stress- 
strain  curve  must  pass  through  the  last  measured  data  point.  Thus,  this 
data  point  is  used  in  calculating  the  constants  B and  C in  Eq.  (2.3). 

For  this  extrapolation  approach,  the  strain  energies  Umax  and  U*  in  Fig. 
2- 22b  coincide.  This  approach  of  using  a linear  extension  from  the  last 
defined  data  point  is  somewhat  questionable  in  that  the  material  may  be 
represented  as  stiffer  than  it  actuary  is  at  high  strain  energies. 
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(a)  STRESS-STRAIN  CURVE  EXTRAPOLATION 


(b)  MECHANICAL  PROPERTY  VERSUS  STRAIN  ENERGY  CURVE.  EXTRAPOLATION 


FIGURE  2-2?  LINEAR  STRESS-STRAIN  CURVE  EXTRAPOLATION  WITH 
NONZERO  SLOPE  EQUAL  TO  SLOPE  AT  LAST  DATA  POINT 
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The  mathematical  technique  of  extrapolation  is  of  course  inexact. 

The  best  of  the  linear  extensions  with  nonzero  slope  cannot  be  determined 
because  comparison  with  measured  stress-strain  data  is  obviously  not 
available  or  else  the  extrapolation  would  not  be  necessary.  However,  a 
choice  between  the  two  approaches  is  best  made  on  an  individual  material 
basis.  Both  approaches  should  be  used  in  the  material  model  to  predict 
material  behavior  under  multiaxial  or  uniaxial  loading  conditions  for 
which  extrapolation  is  needed  and  for  which  measured  data  are  available. 
The  approach  which  leads  to  the  best  correlation  between  predicted  and 
measured  behavior  should  be  used  in  the  model  for  that  particular  mate- 
rial. The  various  material  models  in  Section  2.2  and  2.3  are  collectively 
referred  to  hereafter  as  the  Jones-Nel son-Morgan  nonlinear  material  model. 

2.4  THE  JNMDATA  COMPUTER  PROGRAM 

The  basic  objective  of  the  JNMDATA  computer  program  is  to  automate 
the  use  of  the  Jones-Nel son-Morgan  nonlinear  material  model.  In  partic- 
ular, the  Jones-Nel son-Morgan  material  model  parameters  are  calculated 
from  input  measured  stress-strain  curve  data  points,  the  associated  me- 
chanical property  versus  strain  energy  curves  are  plotted,  and  the  im- 
plied stress-strain  curves  are  plotted.  In  both  types  of  plots,  the 
actual  input  data  are  also  plotted  so  that  a vivid  comparison  of  the 
material  model  approximation  to  the  data  is  available. 

The  actual  stress-strain  curve  data  points  in  Fig.  2-23a  are  con- 
verted in  the  JNMDATA  program  to  mechanical  property  versus  strain  energy 
data  points  and  are  plotted  in  the  form  of  Fig,  2-23b.  However,  these 
two  types  of  data  are  tedious  to  plot  by  hand.  Thus,  the  plotting  fea- 
ture of  the  program  for  these  two  types  of  information  is  quite  attrac- 
tive. Moreover,  different  Jones-Nel son-Morgan  model  fits  to  the  data  can 
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be  obtained  as  in  Fig.  2-5  depending  on  the  specific  data  points  used  tn 
calculate  the  Jones-Nel son-Morgan  parameters.  The  Jones-Nel son-Morgan 
model  parameters  are  calculated  in  the  program  after  the  data  points  to 
be  used  are  specified  in  the  input  data.  Thus,  the  many  model  fits  to 
the  data  with  many  combinations  of  points  can  be  calculated  and  plotted 
as,  for  example,  in  Fig.  2-24.  The  best  fit  of  all  the  possible  fits  is 
the  one  for  which  both  the  stress-strain  data  and  the  mechanical  property 
versus  strain  energy  data  are  simultaneously  best  approximated  with  the 
respective  curves. 

Generally,  the  problem  is  to  obtain  .he  best  fit  to  the  material 
property  data  for  a single  material.  That  material  is  usually  more  com- 
plicated than  an  isotropic  material  and  has  different  stress-strain  be- 
havior at  different  temperatures.  Thus,  the  process  of  obtaining  these 
best  fits  for  all  properties  at  all  temperatures  involves  a lot  of  work 
even  with  the  JNMDATA  program  (and  would  be  almost  unthinkable  without 
such  an  automating  tool).  Many  successive  computer  runs  are  necessary 
In  order  to  achieve  the  best  fits.  The  two  specified  data  points  for 
each  property  at  each  temperature  as  In  Figs.  2-4  and  2-5  are  changed  in 
each  computer  run  along  with  the  specified  initial  slope  of  the  stress- 
strain  curve  until  both  the  mechanical  property  versus  strain  energy 
curves  and  the  associated  implied  stress-strain  curves  fit  the  data  by 
passing  through  them  as  in  Fig.  2-24. 

The  basic  steps  in  the  JNMDATA  computer  program  are  shown  schemati- 
cally in  the  simple  flow  chart  of  Fig.  2-25  and  are  described  as: 

(1)  Read  input  stress-strain  curve  data  as  individual  pairs  of 
stress  and  corresponding  strain  values  along  with  the  identi- 
fication of  the  points  to  be  used  in  the  model. 
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(a)  MECHANICAL  PROPERTY 
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(b)  STRESS-STRAIN 
RESULTS 


FIGURE  2-24  JONES-NELSON  NONLINEAR  MATERIAL  MODEL  FOR 
DATA  OF  FIGURE  2-23 
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FIGURE  2-25  JNMDATA  COMPUTER  PROGRAM  FLOW  CHART 


(2)  Calculate  the  Jones-Nel son-Morgan  model  parameters  from  the 
data  at  the  specified  points. 

(3)  Punch  material  property  cards  for  the  Jones-Nel son-Morgan 
model  in  the  SAAS  HIM  format  for  direct  use. 

(4)  Plot  all  the  mechanical  property  versus  strain  energy  curves 
for  each  mechanical  property  along  with  all  the  associated 
Input  data. 

(5)  Plot  all  the  implied  stress-strain  curves  for  each  mechanical 
property  approximation  along  with  the  associated  input  stress- 
strain  data. 

The  program  output  In  the  form  of  plots  is  inspected  to  see  if  the  model 
fits  for  the  specified  set  of  data  points  to  be  used  are  good  fits  of 
the  input  data.  If  not,  then  the  set  of  data  points  is  revised  based 
on  experience  and  reinput  to  another  run  of  the  program.  When  one  prop- 
erty is  satisfactorily  fit,  then  that  property  can  be  omitted  from  sub- 
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sequent  computer  runs.  That  is,  for  most  materials,  many  properties  are 
approximated  - some  better  and  sooner  than  others.  For  example,  a trans- 
versely Isotropic  material  with  different  moduli  in  tension  and  compres- 
sion under  axi symmetric  load  has  five  independent  mechanical  properties 

in  tension  (E  , E , v , v , and  E^  ) and  five  independent  prop- 
rt  zt  ™t  ztit  rzt  45 

erties  in  compression  (E  , Ez  , vrQ  , vzQ  , and  E^z  ).  Each  of  these 

rc  c °c  °c  c 

properties  is  plotted  as  in  Fig.  2-24  except  for  the  Poisson's  ratios 
which  are  plotted  as  \j  versus  strain  energy,  etransverse  versus  t'nomil, 

and  °normal  versus  "transverse  as  1n  Fig‘  2'26*  Plots  of  1ndiv1dual 
mechanical  properties  can  be  made  as  long  as  sufficient  data  are  avail- 
able, e.g.,  if  v . is  desired,  the  corresponding  direct  modulus  E^ 
rut  rt 

must  also  be  calculated  and  plotted.  Thus,  the  satisfactorily  approxi- 
mated mechanical  properties  can  be  omitted  from  the  subsequent  computer 
runs  so  less  time  and  money  is  expended.  The  computer  charges  are  very 
low  per  run  (typically  about  $10  on  the  SMU  CDC  CYBER  72),  but  the  elapsed 
time  to  generate  the  plots  is  significant.  At  SMU,  the  mechanical  prop- 


FIGURE  2-26  POISSON'S  RATIO  CURVES 


75 


erty  versus  strain  energy  curve  and  its  associated  stress-strain  curve 
are  plotted  in  five  minutes.  Thus,  a full  set  of  curves  for  a trans- 
versely isotropic  multimodulus  material  is  plotted  in  sixty  minutes  (5 

45  45 

minutes  each  for  E , E , E , E , E , and  E ' and  7 1/2  minutes 
rt  zt  rzt  rc  zc  rzc 

each  for  vrQ  , vz0  , vrQ  , and  vzQ  ).  Accordingly,  any  needed  changes 
t t c c 

are  made  more  quickly  (less  elapsed  calendar  time)  if  only  the  affected 
curves  are  replotted.  This  time  constraint  would  be  much  less  important 
for  an  oscilliscope-type  plotter  than  for  SMU's  pen  on  paper  (Zeta)  plot- 
ter. Moreover,  an  interactive  graphics  terminal  would  be  even  more  effec- 
tive because  changes  could  be  made,  their  effects  observed,  and  new 
changes  made  - all  in  rapid  sequence. 

In  all  calculations  from  input  data,  the  potential  data  pitfalls 
discussed  in  Section  2.2.2,  Implementation  of  the  Material  Model,  are 
taken  Into  account.  That  is,  the  program  has  "trapdoors"  or  "fixes"  as 
appropriate  for  each  of  the  pitfalls  mentioned.  In  fact,  several  of  the 
pitfalls  were  discovered  during  the  development  of  the  JNMDATA  program. 

One  such  pitfall,  a "low"  first  stress-strain  curve  data  point  is  seen 
in  Fig.  2-24  which  is  the  with  grain  tensile  stress-strain  curve  at  3000°F 
for  ATJ-S(WS)  Graphite  from  Fig.  5 of  the  SoRI  Graphite  Data  Book  [2-8]. 
Similar  low  first  points  are  seen  in  Figs.  3,  4,  6,  7,  8,  9,  and  10  of 
the  SoRI  Reference.  We  avoid  these  pitfalls  by  inputing  an  appropriate 
initial  slope  which  is  higher  than  the  slope  to  the  low  first  point  and 
in  fact  higher  than  the  slope  to  the  second  data  point. 

An  Interesting  result  from  the  JNMDATA  program  Is  that  we  now  realize 
that  the  actual  initial  slope  of  a stress-strain  curve  is  usually  higher 
than  the  slope  we  graphically  "pick  off"  the  plotted  stress-strain  curve. 
For  example,  the  slope  we  would  pick  off  the  or  versus  er  curve  at  70° 
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(Fig.  12  in  the  SoRI  Graphite  Data  Book  [2-8])  is 

SLOPE  = f = 1 . 752x1 06psi 

The  slope  from  the  properly  fit  mechanical  property  versus  energy  data  is 
3 . 6x1 O^psi ! Thus,  our  eyes  deceive  us  when  we  attempt  to  visualize  slope 
in  the  form  of  secant  modulus  to  a stress-strain  curve.  The  JNMDATA  pro- 
gram is  useful  in  remedying  this  situation  because  of  the  ability  to  plot 
the  mechanical  property  (slope)  directly. 

In  summary,  the  JNMDATA  computer  program  Is  an  exceptionally  useful 
adjunct  to  the  Jones-Nel son-Morgan  nonlinear  material  model.  The  tedious 
calculations  for  the  model  are  automated  and  the  associated  mechanical 
property  versus  strain  energy  curves  and  implied  stress-strain  curves 
are  plotted.  With  this  tool,  rapid  and  highly  visible  material  modeling 
can  be  performed. 


] 

J 

5 3.  MODELING  OF  GRAPHITIC  MATERIALS 

;i: 

The  use  of  the  Jones-Nelson  and  Jones-Nel son-Morgan  nonlinear  mate- 
rial models  for  analysis  of  the  behavior  of  graphitic  materials  is  divid- 
ed in  three  major  parts.  First,  the  Jones-Nelson  model  is  applied  to 
analysis  of  the  Southern  Research  Institute  thermal  stress  disk  test. 
There,  the  model  Is  used  to  predict  the  disk  deformations  which  are  then 
compared  with  the  measured  deformations.  Second,  the  Jones-Nel  son-Morgan 

. i 

■ model  is  applied  to  analysis  of  a nosetip  subjected  to  the  Air  Force 

; i 

Flight  Dynamics  Laboratory  50  MW  test  environment.  Finally,  the  appli- 
i.  cation  of  Jones-Nelson  and  Jones-Nel  son-Morgan  modeling  tc  graphitic 

<•  materials  is  summarized. 

- 3.1  THERMAL  STRESS  DISK  TEST  CORRELATION 

3.1.1  INTRODUCTION 

ill 

ATJ-S  graphite  is  a transversely  isotropic  granular  composite  mate- 
*:•  rial  which  is  made  in  cylindrical  billet  form  as  shown  in  Fig.  3-1  along 


with  the  coordinates  used  in  this  section.  The  flake-liKe  graphite  parti 
cies  are  aligned  in  planes  of  isotropy  during  billet  compaction  in  the  in- 
direction, The  resulting  material  stress-strain  behavior  is  highly  tem- 
perature-dependent as  displayed  for  the  r-0  plane  in  Fig.  3-2.  There, 
the  boxes  are  actual  experimental  data  reported  by  Starrett  and  Pears 
[3-1]  and  the  curves  are  the  Jones-Nelson  nonlinear  material  model  [3-2 
thru  3-6]  fits  to  the  data.  In  Fig.  3-2,  the  tension  behavior  becomes 
stlffer  as  the  temperature  approaches  2000°F  and  even  stiffer  at  3000°F. 
However,  the  stress-strain  curve  at  3000°F  is  slightly  lower  for  high 
strains  than  at  2000°F.  The  compression  behavior  for  the  r-e  plane  mono- 
tonically  becomes  more  flexible  as  the  temperature  Increases  from  3000°F 
to  5000°F.  Moreover,  at  every  temperature,  the  stress-strain  behavior  is 
different  under  tension  loading  than  under  compression  loading. 

The  objective  of  this  section  is  to  extend  the  qualification  of  the 
Jones-Nelson  nonlinear  material  model  [3-2  thru  3-6]  from  mechanical  load 
1 ng  problems  for  homogeneous  bodies  to  thermal  loading  problems  for  non- 
homogeneous  bodies.  The  previous  level  of  qualification  of  the  model  is 
to  uniaxial  off-axis  mechanical  loading  problems  [3-3  and  3-6]  and  to 
biaxial  mechanical  loading  problems  in  principal  material  directions  for 
ATJ-S  graphite  [3-2  and  3-6] . The  multimodulus  character  (different  modu 
li  or  stiffnesses  in  tension  than  in  compression)  is  treated  by  Jones  and 
Nelson  [3-5  and  3-6],  but  nonhomogeneous  bodies  due  to  temperature-depend 
ent  behavior  have  not  been  addressed  with  this  model.  Part  of  the  moti- 
vation for  this  extension  of  the  model  is  the  use  of  ATJ-S  graphite  In 
reentry  vehicle  nosetips  where  substantial  temperature  gradients  exist. 
These  temperature  gradients  lead  to  thermal  stresses  for  which  the  asso- 
ciated nonlinear  deformations  must  be  predicted  in  order  to  rationally 
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and  accurately  design  nosetlps  and  other  structural  elements. 

The  simplest  vehicle  to  accomplish  the  objective  of  qualifying  the 
model  for  thermal  loading  problems  is  the  Southern  Research  Institute 
(SoRI)  thermal  stress  disk  test  [.3-7].  In  this  test,  an  annular  disk  is 
rapidly  heated  at  its  outer  diameter  resulting  in  an  outside  to  inside 
diameter  temperature  gradient.  The  outside  portion  of  the  disk  tends  to 
expand  more  than  the  inside  so  thermal  stresses  are  generated.  In  the 
circumferential  direction,  the  stresses  are  compressive  near  the  outer 
diameter  and  tensile  near  the  inner  diameter.  The  level  of  both  the  com- 
pressive and  tensile  stresses  is  important  because  the  original  purpose 
of  the  SoRI  thermal  stress  disk  test  was  to  quickly  and  easily  fail  disks 
of  many  different  graphites  so  their  relative  merits  can  be  assessed. 

The  failures  almost  invariably  occur  in  the  region  of  tensile  stress. 
Accordingly,  the  compressive  stresses  must  be  high  enough  to  generate 
(self-equilibrating)  tensile  stresses  which  will  cause  failure.  The  ten- 
sile stresses  are  too  low  for  the  simple  flat  annular  disk  and  the  rim 
disk  in  Fig.  3-3.  That  is,  not  enough  disk  material  is  in  compression 
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to  force  high  enough  tensile  stresses  for  fracture  to  occur.  However, 
the  disks  with  a wedge-shaped  cross  section  in  Fig.  3-3  do  have  high 
enough  tensile  stresses  to  cause  fracture.  The  slanted  wedge  disk  has 
an  Inclined  inner  diameter  surface  to  provide  a well-defined  target  for 
the  laser  diameter  measuring  device  which  will  be  discussed  later.  The 
inner  diameter  of  the  wedge  disk  in  Fig.  3-3  is,  however,  sufficiently 
well-defined  to  provide  a good  target.  Thus,  the  simple  wedge  disk  is 
used  in  the  theoretical -experimental  correlation  studies. 

This  section  is  divided  in  two  major  subsections:  (1)  measurement 

of  temperatures  and  deformations  and  (2)  predicted  deformations,  stresses, 
and  strains.  In  the  measurements  section  (Section  3.1.2),  the  overall 
test  setup  is  described  along  with  the  character  of  the  individual  mea- 
surements which  are  made.  In  the  predictions  section  (Section  3.1.3), 
the  Jones-Nelson  nonlinear  material  model  and  its  extensions  for  the  pre- 
sent problem  are  described  along  with  the  modeling  of  ATJ-S  graphite. 
Moreover,  the  predicted  diameter  changes  are  compared  with  the  measured 
values.  Finally,  the  stress  state  In  the  disk  is  displayed  and  discussed. 
3.1.2  MEASUREMENTS  OF  TEMPERATURES  AND  DEFORMATIONS 
3. 1.2,1  Overall  Test  Setup 

The  SoRI  thermal  stress  disk  test  [3-7]  Is  basically  a tension  test 
in  which  the  loads  are  induced  solely  by  thermal  deformation.  Radial 
temperature  gradients  are  generated  in  the  disk  by  induction  heating  as 
shown  in  Fig.  3-4.  The  power  from  the  induction  heating  coils  is  con- 
centrated In  the  outer  periphery  of  the  disk  (about  85%  in  the  outer  .10 
inch)  while  the  inside  remains  relatively  cool.  Many  factors  enter  in 
the  proper  matching  of  the  power  supply  to  the  test  specimen,  but  the 
major  considerations  include  coil  design,  specimen  material,  specimen 
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FIGURE  3-4  SCHEMATIC  OF  SoRI  THERMAL  STRESS  DISK  TEST 
design,  and  speclmen-to-coll  proximity.  The  disk  Is  made  with  Its  dia- 
meter in  the  graphite  plane  of  Isotropy  shown  in  Fig.  3-1.  The  disk  Is 
not  fastened  to  the  test  fixture,  but  Is  positioned  horizontally  as  in 
Fig.  3-4  and  rests  on  a smooth  surface. 

Several  measurements  are  made  during  a single  test  run  for  an  Indi 
vidua!  specimen.  These  measurements  are  (1)  a temperature  measurement 


at  the  Inner  diameter,  (2)  a temperature  measurement  at  a point  near  the 
outer  diameter,  and  (3)  measurements  of  the  diametral  deformations  along 
two  orthogonal  axes  of  the  internal  hole.  Each  of  these  measurements  is 
discussed  in  the  following  subsections. 

3. 1.2. 2 Inner  Diameter  Temperature  Measurement 

The  temperature  at  the  disk  inner  diameter  is  measured  with  a thermo- 
couple  probe.  The  fine  thermocouple  wires  (10  x 10  Inches)  are  not 
welded  together  but  Instead  are  spring  pressed  on  the  disk  so  the  graphite 
makes  the  junction.  This  procedure  Is  used  to  avoid  a weld  bead  which 
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invariably  leads  to  distortion  of  the  parent  material  thermal  characteris- 
tics and  to  errors  of  several  hundred  degrees  at  the  heating  rates  used 
in  this  test. 

The  thermocouple  probe  is  calibrated  by  obtaining  a steady  state 
reading  between  1 200°F  and  2000°F  for  a graphite  specimen  and  comparing 
that  reading  with  the  reading  from  an  optical  pyrometer  which  is  corrected 
for  emittance.  The  error  due  to  the  graphite  junction  is  on  the  order  of 
70°F  at  1500°F.  Then,  typical  rapid  heating  runs  are  made  and  the  thermo- 
couple probe  reading  is  compared  to  a micro-optical  reading  which  Is  cal- 
ibrated to  a hand-held  optical  pyrometer.  Thus,  the  basic  calibration  is 
for  a range  from  the  upper  bound  at  the  steady  state  reference  to  the 
lower  bound  at  the  calibration  to  the  hand  optical  pyrometer  through  the 
micro-optical  recorder.  The  temperature  at  the  Inside  diameter  is  known 
to  within  5°F  on  any  run. 

3.1. 2. 3 Outer  Diameter  Temperature  Measurement 

The  temperature  at  the  disk  outer  diameter  is  measured  with  a micro- 
optical  recorder  as  a function  of  time.  Thus,  the  temperature  reading  is 
subject  to  the  errors  normally  associated  with  optically  measuring  the 
temperature  of  a gray  body  plus  the  errors  peculiar  to  the  particular 
specimen  and  to  this  system.  The  latter  errors  Include  those  peculiar  to 
transient  measurements  plus  those  resulting  from  the  existence  of  a 200°F 
temperature  gradient  over  the  field  of  view  of  the  recorder.  That  Is, 
the  recorder  is  focused  on  a spot  as  in  Fig.  3-4  on  the  top  surface  of 
the  disk  since  the  outer  diameter  surface  is  covered  by  the  Induction 
coils.  Since  the  spot  is  necessarily  of  finite  width  (about  45  x 10 
inches)  and  is  located  in  a region  of  high  radial  thermal  gradient  (about 
4°F  per  10-3  inch),  inaccuracies  in  temperature  measurement  Inevitably 
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occur. 


Various  calibration  efforts  for  this  measurement  are  discussed  by 
Pears  and  Starrett  [3-7].  They  conclude  that  the  micro-optical  recorder 
can  be  used  to  measure  temperatures  up  to  about  5500°F  on  the  disk  outer 
surface  and  could  be  extended  to  even  higher  temperatures.  The  inherent 
accuracy  of  the  readings  is  about  + 50°F  up  to  5000°F.  However,  the 
recorder  slewing  speed  is  approximately  the  same  as  the  temperature  rise 
in  many  tests  so  the  final  outer  diameter  temperature  measurements  are 
often  Ignored. 

3. 1.2. 4 Diametral  Deformation  Measurements 

The  changes  in  the  disk  inner  diameter  are  measured  along  x-  and 
y-axes  In  Fig.  3-4  by  use  of  a laser  tracking  device.  The  laser  beams 
are  focused  and  located  with  a complex  prism  system  described  by  Pears 
and  Starrett  [3-7].  Basically,  the  two  orthogonal  diameter  change  mea- 
surements are  different  because  the  electronic  circuitry  is  slightly  dif- 
ferent (each  corresponding  component  of  the  two  circuits  is  different 
within  the  manufacturers'  tolerances).  The  output  of  the  x-  and  y-mea- 
surements  is  recorded  on  the  same  recorder  with  offset  Initial  points. 
These  initial  points  are  located  by  use  of  a tick  mark  placed  on  the 
recorder  when  the  induction  coils  are  turned  on.  The  tick  mark  is  not, 
however,  recorded  precisely  when  the  coils  are  turned  on,  Pears  and 
Starrett  [3-7]  conclude  that  the  maximum  diameter  change  error  for  a sin- 
gle experiment  Is  72  x l(f®  inches  plus  the  timer  error  which  in  these 
correlations  is  eliminated  by  adjustment  of  the  zero  time  mark. 

3. 1.2. 5 Summary 

The  inner  diameter  temperature  measurements  are  used  to  determine 
the  approximate  temperature  distribution  In  the  disk.  That  Is,  the  the- 
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oretical  temperature  distributions  at  various  times  are  shifted  in  time 
until  the  measured  time  when  the  predicted  inner  diameter  temperature 
agrees  with  the  measured  inner  diameter  temperature.  The  outer  diameter 
temperature  is  recorded  at  such  a speed  that  the  results  are  not  meaning- 
ful. Otherwise,  they  could  be  used  to  further  adjust  the  temperature 
distributions  in  the  manner  discussed  for  the  inner  diameter  temperature 
measurements.  The  diametral  deformation  measurements  are  the  primary 
measurements  and  are  half  of  the  basis  for  comparison  of  measured  and 
predicted  deformations. 

3.1.3  PREDICTED  DEFORMATIONS,  STRESSES,  AND  STRAINS 

The  predicted  deformations,  stresses,  and  strains  of  the  wedge- 
shaped  annular  disk  are  obtained  with  the  Jones-Nelson  nonlinear  material 
model  which  Is  described  in  the  next  subsection.  Then,  the  Jones-Nelson 
model  of  ATJ-S  graphite  Is  displayed  and  discussed.  Next,  the  predicted 
inner  diameter  changes  are  compared  with  measured  values.  Finally,  the 
stresses  and  strains  in  the  disk  are  described. 

3. 1.3.1  Jones-Nelson  Nonlinear  Material  Model 

The  Jones-Nelson  nonlinear  material  model  Is  described  in  References 
3-2  thru  3-6.  Here,  we  merely  reiterate  the  details  essential  to  our  use 
of  the  model.  Basically,  the  various  secant  moduli  and  Poisson's  ratios, 
all  of  which  are  called  mechanical  properties,  in  the  orthotropic  stress- 
strain  relations  are  approximated  with 

Mechanical  Property^  » - B^(U/U0^)  ^ (3.1) 

where  the  A^  are  the  elastic  values  of  the  1^  mechanical  property,  the 
B.|  and  are  related  to  the  initial  curvature  and  initial  change  of  cur- 
vature, respectively,  of  the  stress-strain  curve  [3-3]  (slightly  different 
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interpretations  exist  for  Poisson's  ratios),  and  U is  the  strain  energy 
density  of  an  equivalent  elastic  system  at  each  stage  of  nonlinear  defor- 
mation: 


U - (°rEr  * Vz  + o0ce  * T„Yr2)/2 


(3.2) 


The  strain  energy  density  U in  Eq.  (3.1)  is  normalized  by  U so  that 

°1 

and  are  dimensionless. 

The  Jones-Nelson  model  is  actually  much  more  complicated  than  Eq. 
(3.1).  For  example,  In  states  of  mixed  tension  and  compression,  the 
strain  energy  could  be  the  total  strain  energy;  the  tension  strain  energy, 
Ut,  for  some  properties  and  the  compression  strain  energy,  U , for  other 
properties,  or  some  weighted  combination  of  the  tension  energy  and  the 
compression  energy.  The  strain  energy  density  used  in  this  section  is 


uw  " <ut  + Uc)/U 


(3.3) 


where  Uw  Is  the  weighted  strain  energy.  Moreover,  all  coefficients  in 
Eq.  (3.1)  have  different  values  in  tension  than  In  compression.  The  final 
combination  of  mechanical  properties  Is  determined  in  an  iteration  pro- 
cedure wherein  the  stresses  In  principal  stress  directions  are  a key  fac- 
tor [3.5]. 

3, 1.3.2  ATJ-S  Graphite  Mechanical  Properties 

The  nonlinear  stress-strain  behavior  in  the  plane  of  isotropy  is 
depicted  in  Fig.  3-2.  Since  the  disk  Is  stressed  primarily  in  the  r-0 
plane,  we  are  concerned  mainly  with  the  properties  In  that  plane.  How- 
ever, many  other  properties  are  necessary  to  perform  an  analysis.  The 
values  of  , and  are  given  as  a function  of  temperature  for  all 

pertinent  mechanical  properties  in  Table  3-1  (the  value  of  IL  Is  1 psl 

°1 

for  all  properties).  Also,  the  coefficients  of  thermal  expansion  are 
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TABLE  3-1 , continued 


TEMPERATURE 

MECHANICAL 

PROPERTY 

A 

B 

C 

S 

s 

2.1 0x1 06ps i 
1.45x106ps1 

.0510 

.0254 

.660 

.921 

vret 

.140 

0. 

1. 

v „ 

.116 

0. 

1. 

z0t 

E45 

rZt 

2.11x1 06psi 

.0380 

.790 

3000° F 

.198 

Er 

3.60x1 0°ps1 

.350 

C 

Ez 

C 

1 .80x1 06ps1 

.194 

.291 

H 

.090 

0. 

1. 

V /V 

.090 

0. 

1. 

z6c 

e"5 

3 . 27x1 06ps1 

.310 

.240 

S 

S 

2.40x1 0°psi 
1.40xl06psi 

.277 

.133 

~ .512 
.395 

H 

.145 

0. 

1. 

X 

.120 

0. 

1. 

E45 

rzt 

2.04xl06psi 

.158 

.388 

3500°F 

Er, 

3.60x10°psi 

.316 

.220 

c 

Ez, 

2.00x106psi 

.158 

.347 

c 

1. 

H 

.0950 

0. 

VZ9r 

.0950 

0. 

1. 

c 

-45 

trz1. 

3.27x106psi 

.270 

.290 

c 

h 
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given  in  Table  3-2.  Both  sets  of  data  are  obtained  from  Ref.  3-1.  The 


h 

ij 


direct  moduli  Er  and  Ez  are  shown  in  stress-strain  curve  form  in  Ref. 

3-1  along  with  tables  of  (constant)  Poisson's  ratios  and  coefficients  of 

45 

thermal  expansion.  The  values  for  E^z  are  obtained  by  observing  the  rela- 

45 

tlon  between  the  coefficients  in  Eq.  (3.1)  for  Erz»  Er,  and  Ez  measured 
for  ATJ-S  graphite  at  room  temperature  by  Jortner  [3-8]  and  applying  that 
relation  to  the  measured  data  for  Er  and  Ez  at  elevated  temperatures. 

The  boxes  In  Fig.  3-2  are  actual  data  from  Ref.  3-1,  and  the  curves 
are  Jones-Nelson  curves  through  the  data.  The  curves  are  obtained  by 
use  of  the  JNMDATA  computer  program  in  which  all  the  stress-strain  curve 
data  are  input  along  with  instructions  for  which  points  to  use  In  the 
calculation  of  A^ , B^,  and  C^.  Obviously,  very  good  fits  to  the  data  can 
be  generated  with  the  JNMDATA  computer  program.  Moreover,  all  stress- 


TABLE  3-2 

COEFFICIENTS  OF  THERMAL  EXPANSION  FOR  ATJ-S(WS)  GRAPHITE 
AS  A FUNCTION  OF  TEMPERATURE 


TEMPERATURE 

°F 

10"5  per  °F 

cx2 

10"5  per  °F 

70 

.1200 

.1720 

2000 

.1915 

.2471 

3000 

.2194 

.2727 

3500 

.2315 

.2864 

4000 

.2414 

.2997 

4500 

.2547 

.3174 
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strain  curves  and  mechanical  property  versus  strain  energy  curves  can  be 
plotted  easily. 

The  data  in  Tables  3-1  and  3-2  are  input  to  a version  of  the  SAAS 
III  finite  clement  stress  analysis  computer  program  [3-9]  which  has  been 
modified  to  use  the  Jones-Nelson  material  model.  In  the  use  of  the  new 
version  called  SAAS  HIM,  stress-strain  data  at  temperatures  between  the 
values  listed  in  Tables  3-1  and  3-2  must  be  generated.  If  the  Jones- 
Nelson  parameters  A^,  B^,  and  C^  for  T = 70°F  and  T = 2000°F  are  inter- 
polated for  a specific  temperature,  say  1403°F  which  is  the  inner  dia- 
meter temperature  at  t = 1.9  seconds,  then  the  stress-strain  curve  la- 
beled "PARAMETER  INTERPOLATION  FOR  1403°F"  in  Fig.  3-5  results.  However, 
this  stress-strain  curve  doesn't  lie  between  the  two  curves  from  which 
it  is  Interpolated.  On  the  other  hand,  if  the  nonlinear  mechanical  prop- 
erties calculated  from  Eq.  (3.1)  for  T = 70°F  and  T = 2000°F  are  inter- 
polated (i.e.,  E at  7Q°F  and  at  2000°F),  then  the  curve  labeled  " PROP- 
rt 

ERTY  INTERPOLATION"  in  Fig.  3-5  results.  Finally,  if  the  points  at  the 
same  strain  energy  on  the  two  stress-strain  curves  are  interpolated,  then 
the  curve  labeled  "STRESS-STRAIN  CURVE  INTERPOLATION"  in  Fig.  3-5  results. 
Obviously,  the  latter  two  approaches  coincide  for  all  practical  purposes. 
The  property  interpolation  approach  is  used  in  SAAS  HIM  because  it  is 
simpler  to  program  and  cheaper  to  execute  than  the  stress-strain  curve 
interpolation  approach. 

3, 1.3. 3 Inner  Diameter  Change  Predictions 

Times  late  in  the  test  run  are  selected  for  correlation  of  predicted 
and  measured  results.  Then,  the  disk  should  be  deforming  nonlinearly, 
i.e.,  the  stresses  should  be  Inelastic.  The  power  in  the  Induction  coils 
does  not  reach  the  nominal  levels  used  to  predict  the  temperature  distri- 
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FIGURE  3-5  TEMPERATURE  INTERPOLATION  OF  STRESS-STRAIN  BEHAVIOR 
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butions  at  various  times  with  the  MOATS  heat  transfer  computer  program 
[3-10]  (in  which  convection  is  not  considered).  Hence,  the  temperature 
gradients  are  lower  at  all  times  than  predicted.  We  account  for  this 
discrepancy  by  assuming  that  a predicted  temperature  distribution  at  one 
time  with  a specific  inner  diameter  temperature  is  the  actual  temperature 
distribution  at  another  time  when  the  predicted  inner  diameter  temperature 
is  the  same  as  the  measured  temperature.  The  validity  of  this  time-power 
shift  hypothesis  would  be  even  more  plausible  if  we  could  also  match  pre- 
dicted and  measured  outer  diameter  temperatures.  However,  the  outer  dia- 
meter temperature  measurement  with  e micro-optical  recorder  is  not  suffi- 
ciently accurate  to  be  relied  upon  because  the  temperature  rise  is  too 
fast  for  the  recording  device  speed  used  for  all  measurements.  Thus,  we 
depend  on  the  inner  diameter  temperature  correlation  to  determine  the 
times  at  which  the  predicted  temperature  distributions  are  applicable. 
These  times,  t ■ 1.2,  1.55,  and  1.9  seconds,  correspond  to  measured  inner 
diameter  temperatures  of  889°F,  1135°F,  and  1410°F  whereas  the  predicted 
temperatures  are  886°F,  1130°F,  and  1402°F,  respectively.  Revision  of 
these  times  to  force  the  predicted  temperatures  to  correspond  more  close- 
ly to  the  measured  temperatures  is  regarded  as  unnecessary.  The  actual 
calibrated  -inner  diameter  temperatures  are  plotted  as  a function  of  time 
in  Fig.  3-6.  There,  dots  are  used  to  denote  the  times  and  correspond- 
ing temperatures  at  which  the  deformations,  stresses,  and  strain  are  pre- 
dicted. The  predicted  deformations  should  be  slightly  too  low  because  the 
temperature  gradient  increases  with  time  and  the  predicted  temperatures 
are  lower  than  the  measured  temperatures. 

The  measured  inner  diameter  changes  are  shown  in  Fig.  3-7  for  two 
orthogonal  directions  as  a function  of  time.  The  change  of  diameter  in 
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FIGURE  3-6  INNER  DIAMETER  TEMPERATURE  VERSUS  TIME 
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the  two  directions  is  measured  with  two  electronically  equivalent  but 
not  identical  circuits.  The  measurement  in  the  y-direction  is  much  less 
noisy  than  in  the  x-direction.  However,  both  measurements  are  suffi- 
ciently accurate  for  the  present  correlation  effort  without  calibration. 
The  difference  in  measured  deformations  In  the  x-  and  y-directlons  can 
also  be  attributed  to  the  disk  hole  not  remaining  perfectly  circular.  Of 
course,  the  hole  should  remain  circular  since  the  material  Is  nominally 
Isotropic  in  the  plane  of  the  hole  and  the  temperature  distribution  Is 
axlsymmetrlc  about  the  z-axls  perpendicular  to  the  plane  of  the  disk. 
However,  the  material  does  not  have  perfect  transverse  isotropy  nor  Is 
the  temperature  distribution  perfectly  axlsymmetrlc. 

The  predicted  Inner  diameter  changes  are  plotted  In  Fig.  3-7  for 
t a 1.2,  1.55,  and  1.9  seconds.  These  predictions  are  obtained  with  the 
SAAS  IIIM  finite  element  computer  program  for  which  the  disk  finite  ele- 


ment mesh  Is  shown  in  Fig.  3-8.  The  SAAS  IIIM  program  Is  the  SAAS  III 


FIGURE  3-8  WEDGE-SHAPED  ANNULAR  DISK  AND  FINITE  ELEMENT  IDEALIZATION 
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program  [3-9]  with  the  Jones-Nelson  nonlinear  material  model  incorporated. 
As  is  obvious  from  Fig.  3-7,  the  predicted  deformations  are  quite  close 
to  the  measured  values.  The  actual  numerical  values  are  listed  in  Table 
3-3  wherein  we  observe  that  the  predicted  deformations  are  within  2.2  to 
3.3%  of  the  average  of  the  two  measured  diameter  changes. 

TABLE  3-3 

MEASURED  AND  PREDICTED  INNER  DIAMETER  CHANGES 


TIME 

teconds 

I.D.  TEMPERATURE 
op 

CHANGE  IN  INTERNAL  DIAMETER,  Inches 

MEASURED 

PREDICTED 

X ERROR 

MEASURES 

USED  IN 
CALCULATION 

X-CHANNEL 

Y-CHANNEL 

AVERAGE 

1.2 

869 

866 

.00230 

.00250 

.00240 

.00232 

3.3 

1.55 

1135 

1130 

.00300 

.00320 

.00310 

.00303 

2.2 

1.9 

1410 

1402 

.00365 

.00390 

.00378 

.00367 

2.9 

u 


3. 1.3. 4 Stress  and  Strain  Predictions 

The  predicted  stresses  or  and  oQ  are  shown  along  with  the  correspond- 
ing temperature  distribution  at  t = 1.9  seconds  in  Fig.  3-9.  Although  we 
expect  oQ  to  predominate,  substantial  values  of  ar  exist.  These  radial 
stresses  are  always  tensile  and  can  be  explained  with  the  aid  of  the  free 
body  diagrams  in  Fig.  3-10.  There,  the  view  Is  down  the  z-axis  of  the 
disk  perpendicular  to  the  plane  of  the  disk.  The  shear  stress  is  zero 
by  virtue  of  axial  symmetry  of  the  loading  and  geometry.  The  free  bodies 
extend  through  the  (variable)  thickness  of  the  disk  so  the  shearing  stress 
Trz  zer0,  In  ttie  body  next  to  the  outer  diameter,  the  circumferential 
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FIGURE  3-9  TEMPERATURE  AND  STRESS  DISTRIBUTIONS  AT  t - 1.9  second 
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FIGURE  3-10  FREE  BODY  DIAGRAMS  OF  INNER  DIAMETER 
AND  OUTER  DIAMETER  ELEMENTS 
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the  unloaded  Inner  and  outer  surfaces.  Despite  these  surfaces  of  zero 
radial  stress,  can  have  substantial  values  elsewhere  because  of  the 
small  disk  inner  diameter  and  the  high  circumferential  stresses  to  which 
ar  Is  Inversely  proportional  and  proportional,  respectively. 

The  stress-strain  behavior  becomes  more  and  more  nonlinear  when  the 


test  time  Increases  as  we  see  from  the  increasing  disparity  between  elas- 
tic and  nonlinear  stresses  in  Table  3-4.  The  disparity  Increases  to 
nearly  a factor  of  two  at  t - 1.9  seconds.  This  disparity  Is  depicted 


TABLE  3-4 

PREDICTED  CIRCUMFERENTIAL  STRESSES  AND  STRAINS 


1 

i 


I 


RV5 


graphically  for  circumferential  stresses  and  strains  in  Fig.  3-11. 

There,  the  elastic  stresses  are  too  high  by  nearly  a factor  of  two. 
However,  the  elastic  strains  are  not  as  different  from  the  nonlinear 
strains  as  the  corresponding  stresses.  Moreover,  the  elastic  strains 
are  not  always  higher  than  the  nonlinear  strains.  In  fact,  the  non- 
linear inner  diameter  strain  decreases  by  only  5%  and  the  nonlinear  outer 
diameter  strain  increases  by  about  5%  from  the  respective  elastic  strains! 
Thus,  a failure  criterion  based  on  elastic  strains  Is  much  more  accurate 
than  a failure  criterion  based  on  elastic  stresses.  However,  even  an 
elastic  strain  failure  criterion  is  inaccurate  because  of  the  5%  error 
In  predicting  elastic  strains  and  the  fact  that  the  error  Is  sometimes 
positive  and  sometimes  negative. 

The  stresses,  strains  and  displacements  converge  to  the  final  results 
in  5,  5,  and  6 iterations  at  t ■ 1.2,  1.55,  and  1.9  seconds,  respectively. 


FIGURE  3-11  DEGREE  OF  NONLINEAR  STRE5S-STRAIN  BEHAVIOR 
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Thus,  the  rapidity  of  convergence  depends  on  the  degree  of  nonlinearity 
which  increases  with  test  time.  The  Inner  diameter  stresses  at  t s 1.9 
seconds  in  Table  3-5  oscillate  with  decreasing  amplitude  about  the  final 
results.  However,  the  inner  diameter  strains  and  displacements  monoton- 
Ically  decrease  to  the  final  values. 

The  energies  in  the  elastic  state  are  well  above  the  energies  in  the 
actual  nonlinear  states.  Thus,  the  energies  In  the  iteration  procedure 
must  decrease  as  the  number  of  iterations  Increases.  Conceivably,  the 
elastic  energy  could  be  so  much  higher  than  the  actual  nonlinear  energy 
that  the  Jones-Nelson  model  would  not  converge  because  of  unfavorable 
numerical  characteristics  of  the  various  stress-strain  curves  (first  ris- 
ing, then  falling  as  for  3000°F  in  tension  In  Fig.  3-2).  However,  that 
unfortunate  numerical  behavior  Is  not  found  In  this  problem.  If  such  a 
difficulty  arises,  It  can  be  remedied  by  use  of  the  extended  stress-strain 

TABLE  3-5 

PREDICTED  CIRCUMFERENTIAL  STRESSES  AND  STRAINS  AT  I.D.  ELEMENT  1 
AND  RADIAL  DISPLACEMENTS  AT  I.D.  NODAL  POINT  2 AT  t » 1.9  seconds 


ITERATION 

dg.RSl 

e0 

ur, Inches 

1 

10,608 

.00527 

.001903 

2 

5,722 

.00503 

.001841 

3 

6,543 

.00500 

.001834 

4 

6,345 

.00500 

.001834 

5 

6,390 

.00500 

.001833 

6 

6,380 

.00500 

.001833 
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curve  approaches  due  to  Jones  and  Morgan  [3-11]. 
3.1.4  SUMMARY 


The  Oones-Nel son  nonlinear  material  model  is  extended  from  mechani- 
cal loading  problems  for  homogeneous  bodies  to  thermal  loading  problems 
for  nonhomogeneous  bodies.  The  nonhomogeneity  results  from  a temperature 
gradient  over  a body  with  temperature-dependent  mechanical  properties. 
Moreover,  the  model  is  shown  to  be  valid  for  materials  which  have  highly 
nonlinear  stress-strain  behavior  which  is  different  under  tension  load- 
ing than  under  compression  loading. 

The  vehicle  for  the  verification  of  the  model  extension  is  the  SoRI 
thermal  stress  disk  test.  The  Inner  diameter  changes  of  this  annular 
wedge-shaped  disk  made  of  ATJ-S  graphite  are  predicted  with  the  model  to 
within  about  3%.  The  mechanical  properties  of  ATJ-S  graphite  are  a 
strong  function  of  temperature  and  stress  level.  Thus,  the  present  re- 
sults are  a severe  test  of  the  model  and  are  an  Important  step  in  the 
qualification  of  the  model  for  general  use  in  nonlinear  material  defor- 
mation problems. 

3.2  50  MW  NOSETIP  STRESS  ANALYSIS 
3.2.1  INTRODUCTION 

The  stress  and  strain  field  is  analyzed  in  an  ATJ-S  graphite  nosetip 
subjected  to  the  Air  Force  Flight  Dynamics  Laboratory  50  MW  arc  jet  facil- 
ity environment.  The  nosetip  design  Is  supposedly  representative  of 
stress  states  that  exist  during  reentry.  Moreover,  the  nosetip  model 
is  supposed  to  fail  during  the  test.  The  specific  test  condition  analyzed 
is  50  MW  Run  R41-016  at  t = 1.60  seconds.  The  nosetip  is  of  the  shell 
type  (as  opposed  to  the  plug  type)  and  has  a .75"  nose  radius,  1"  over- 
hang, 10°  half  cone  angle,  and  ,15"  wall  thickness  as  shown  In  Fig.  3-12. 


The  temperature  distribution  at  t = 1.60  seconds  displayed  in  Fig.  3-lJ 
was  calculated  with  the  ASTHMA,  ARGEIBL,  and  BLIMP  computer  programs 
at  Aerotherm  Division  of  Acurex  Corporation  and  provided  to  Southern 
Methodist  University  by  AFML.  The  corresponding  surface  pressure  dis- 
tribution in  Fig.  3-14  was  also  provided  to  SMU. 

The  stresses  and  strains  predicted  with  the  Jones-Nel son-Morgan 
nonlinear  material  model  are  displayed  in  contour  plot  form  in  Section 

3.2.2  along  with  the  specific  ATJ-S  graphite  mechanical  properties  which 
are  used.  These  stresses  and  strains  are  compared  in  Section  3,2.3  with 
those  predicted  by  Baker,  Jackson,  Starrett,  and  Budde  [3-12]  who  used 
the  DOASIS  computer  program  [3-13],  The  basic  efforts  of  this  study  are 
summarized  In  Section  3.2.4. 

3.2.2  JONES-NELSON-MORGAN  NONLINEAR  MATERIAL  MODEL  PREDICTIONS 

The  Jones-Nel son-Morgan  model  predictions  are  based  on  the  mechani- 
cal properties  of  ATJ-S  graphite  given  In  Section  3. 2. 2.1.  Then,  In 
Section  3. 2. 2. 2,  the  predicted  elastic  stresses  and  strains  are  displayed. 
Finally,  the  predicted  nonlinear  stresses  and  strains  are  presented  In 
Section  3. 2. 2. 3. 

3. 2. 2.1  ATJ-S  Graphite  Mechanical  Properties 

The  basic  source  of  Information  for  ATJ-S  graphite  mechanical  prop- 
erties Is  the  SoRI  ATJ-S  data  book  [3-1],  However,  those  data  are  mod- 
eled with  the  Jones-Nel son-Morgan  nonlinear  material  model  Instead  of 
with  the  Jones-Nelson  nonlinear  material  model  as  In  the  Thermal  Stress 
Disk  Test  Correlation  In  Section  3.1.  The  reason  for  this  new  approach 
Is  that  the  stresses  and  consequently  the  energy  in  the  nosetlp  Is  too 
high  for  convergence  of  the  basic  Jones-Nelson  model.  That  is,  the  ener- 
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FIGURE  3-12  SHELL  NOSE! IP  GEOMETRY 


gles  are  high  enough  that  the  hump  of  the  stress-strain  curve  In  Fig. 

2-18  is  passed  so  the  model  is  unstable.  Accordingly,  we  must  use  a 
Jones-Nel son-Morgan  model  with  an  extended  stress-strain  curve  that 
doesn't  have  a hump  and  always  has  positive  slope  as  in  Fig.  2-22.  A 
typical  stress-strain  curve  and  mechanical  property  versus  energy  curve 
are  shown  along  with  experimental  data  from  the  SoRI  ATJ-S  data  book 
[3-1]  In  Fig.  3-15  and  3-16  for  behavior  in  the  radial  direction  (o 

rt 

versus  e and  E versus  U)  at  70°F  under  tension  loading.  The  Jones- 
rt  rt 

Nel son-Morgan  nonlinear  material  model  parameters  are  given  for  all 
available  temperatures  In  Table  3-6.  The  associated  coefficients  of 
thermal  expansion  are  given  In  Table  3-7. 

The  implied  slope  version  of  the  Jones-Nel son-Morgan  model  seems  to 
be  the  most  appropriate  of  those  available  although  the  Jones-Nelson- 
Crose  model  might  be  equally  applicable.  In  the  Implied  slope  model, 
the  slope  at  the  last  data  point  Is  arbitrarily  used  as  the  subsequent 
slope  of  the  stress-strain  curve.  In  contrast,  the  subsequent  slope  of 
the  Jones-Nel son-Crose  model  is  determined  by  fitting  a straight  line 
through  the  stress  and  strain  points  corresponding  to  5%  probability  of 
failure  and  50%  probability  of  failure.  Not  having  suitable  failure  data 
to  work  with,  we  chose  the  simpler  Jones-Nel son-Morgan  implied  slope  mod- 
el. However,  the  available  data  are  Insufficient  to  apply  the  implied 
slope  model  without  approximation.  Specifically,  the  data  are  available 

for  Er  and  Ez  in  tension  and  compression  for  all  temperatures  but  not 
45 

for  Erz!  Thus,  we  fit  the  implied  slope  model  to  Er  and  Ez  and  observe 

the  corresponding  values  of  A,  B,  C,  UQ,  U*,  E*,  and  aQ  for  each  case. 

45 

Then,  we  estimate  corresponding  values  for  Erz  based  on  the  relation  of 
^rz  to  ^r  anc*  Ez  at  r00m  temPerature  f°r  Jortner's  data  [3-8],  Simulta- 
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5 LINEAR  STRESS-STRAIN  CURVE  EXTRAPOLATION  WITH 
NONZERO  SLOPE  EQUAL  TO  SLOPE  AT  LAST  DATA  POINT 
FOR  RADIAL  DIRECTION  BEHAVIOR  AT  70°F  IN  TENSION 
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CURVE  EXTRAPOLATION  IN  FIGURE  3-15  FOR 
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TABLE  3-6 


JONES-NELSON-MORGAN  NONLINEAR  MATERIAL  MODEL  PARAMETERS 
FOR  ATJ-S(WS)  GRAPHITE  AS  A FUNCTION  OF  TEMPERATURE 


A 

B 

C 

U* 

psi 

E* 

1 06psi 

°o 

psi 

2.00x1 06psi 

.182 

.337 

11.8 

.708 

2050 

1 .50x1 06ps1 

.226 

.322 

13.1 

.354 

2230 

.no 

0. 

1. 

0. 

0. 

0. 

.090 

0. 

1. 

0. 

0. 

0. 

1.55xI06ps1 

.200 

.330 

12.1 

.480 

1950 

3.00x1 06ps1 

.477 

.125 

37.3 

.336 

4100 

2.10xl06psi 

.451 

.134 

30.3 

.304 

3010 

.065 

0. 

1. 

0. 

0. 

0. 

.065 

0. 

1. 

0. 

0. 

0. 

2.40x1 06p si 

.470 

.130 

29.7 

.310 

3230 

2.00xl06ps1 

.0651 

.583 

14.1 

.825 

2550 

1.60xl06psi 

.165 

.336 

12.3 

.645 

1710 

.130 

0. 

1. 

0. 

0. 

0. 

.107 

0. 

1. 

0. 

0. 

0. 

1.83xl06ps1 

.120 

.460 

11.6 

.660 

2210 

3.00x1 06ps1 

.333 

.203 

37.1 

.338 

5210 

1 . 90x1 06ps1 

.281 

.221 

36.0 

.337 

3830 

.082 

0. 

1. 

0. 

0. 

0. 

.082 

0. 

1. 

0. 

0. 

0. 

2.84x1 06ps1 

.350 

.210 

24.3 

.337 

4120 

TEMPERATURE 


MECHANICAL 

PROPERTY 


70°  F 


2000°F 


>0* 


Z9t 

-45 

Lrzt 


re. 


29  c 

-45 

rz. 


E., 


r9 . 


Z0t 

-45 

'rz* 


"re. 


z0c 

,45 

"rz. 


TABLE  3-6,  continued 


TEMPERATURE 


3000° F 


3500° F 


MECHANICAL 

PROPERTY 


A 

B 

c 

U* 

psi 

E* 

lO^psi 

ao 

psi 

2.10x106ps1 

.0510 

.660 

19.7 

.600 

3990 

1 .45xl0^psi 

.0254 

.921 

16.5 

.350 

3590 

.140 

0. 

1. 

0. 

0. 

0. 

.116 

0. 

1. 

0. 

0. 

0. 

2.11x1 06ps  1 

.0380 

.790 

19.5 

.400 

4830 

3.60x1 06psi 

.350 

.198 

45.8 

.243 

6730 

1.80x1 06ps1 

.194 

.291 

40.2 

.350 

4360 

.090 

0. 

1. 

0. 

0. 

0. 

.090 

0. 

1. 

0. 

0. 

0. 

3.27x1  06ps1 

.310 

.240 

29.9 

.275 

5500 

2.40x1 06psi 

.277 

.212 

70.0 

.293 

6430 

1 .40xl06psi 

.133 

.395 

36.1 

.226 

4360 

.145 

0. 

1. 

0. 

0. 

0. 

.120 

0. 

1. 

0. 

0. 

0. 

2.04xl06ps1 

.158 

.388 

30.0 

.232 

l 

5090 

3.60xl06psi 

.316 

.220 

41.4 

.293 

6560 

2.00xl06p5i 

.158 

.347 

42.2 

.300 

5440 

.0950 

0. 

1. 

0. 

i 

0. 

0. 

.0950 

0. 

1. 

0. 

0. 

0. 

3.27x1 0®psi 

.270 

.290 

24.5 

.240 

5480 

TABLE  3-6,  continued 


W 

l 

I 

fj 


9 


TEMPERATURE 

MECHANICAL 

PROPERTY 

A 

B 

C 

U* 

psi 

E* 

1 06p  si 

°o 

psi 

4000°F 

s 

s 

Vr0t 

H 

c45 

E«t 

S 

Ezc 

H 

X 

45 

r2c 

2.70x1 06psi 
1 .40xl0^psi 
.150 
.125 

1 . 97x1 06ps1 
2.80xl06psi 
1 .80xl06ps1 
.100 
.100 

3. 06xl06ps1 

.282 

.107 

0. 

0. 

.200 

.323 

.153 

0. 

0. 

.270 

.237 

.477 

1. 

1. 

.360 

.212 

.377 

1. 

1. 

.300 

69.5 
35.2 

0. 

0. 

30.0 

42.5 
31.7 

0. 

0. 

20.5 

.069 

.115 

0. 

0. 

.083 

.240 

.267 

0. 

0. 

.250 

8230 

5140 

0. 

0. 

5340 

5730 

4640 

0. 

0. 

4860 

4500° F 

% 

s 

H 

*z0t 

E45 

rzt 

E"c 

Ezc 

H 

X 

e45 

rzc 

2.50x1 06ps1 
1 .30x1 06ps1 
.170 
.150 

1 , 50x1 06ps1 
2.40x1 06ps1 
1 . 80x1 0®ps 1 
.150 
.150 

2,40x1 O^psi 

.300 

.100 

0. 

0. 

.200 

.275 

.244 

0. 

0. 

.300 

.250 

.500 

1. 

1 . 

.400 

.250 

.267 

1. 

1. 

.270 

46.1 

37.4 

0. 

0. 

21 .6 

38.0 
33.8 

0. 

0. 

20.1 

.030 

.060 

0. 

0. 

.035 

.227 

.260 

0. 

0. 

.220 

6700 

5420 

0. 

0. 

4190 

5330 

4160 

0. 

0. 

4030 
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TABLE  3-6,  concluded 


TEMPERATURE 


MECHANICAL 

PROPERTY 


2.09xl0°psi 


u I psi  10  psi  psi 


.250  46.8  .020  6050 


1 .20x1 0°p  s 1 .100  .500  39.2  .040  5400 


0.  0. 


5000° F 


0.  0. 


1 .40x1 0 psi  .200  .400  22.1  .025  4130 


2.00x10°ps1 


.237  31.8  .208  4280 


1 .80x1 0°psi  .408  .170  26.1  .214  3070 

.200  0.  1.  0.  0.  0. 

.200  0.  1.  0.  0.  0. 

1 . 74x1 06ps1  .410  .210  11.3  .210  2190 


TABLE  3-7 

COEFFICIENTS  OF  THERMAL  EXPANSION  FOR  ATJ-S(WS)  GRAPHITE 
AS  A FUNCTION  OF  TEMPERATURE 


TEMPERATURE 

°F 

u 

-5  r 

10  0 per  °F 

-s"z 

10  s per  °F 

70 

.1200 

.1720 

2000 

.1915 

.2471 

3000 

.2194 

.2727 

3500 

.2315 

.2864 

4000 

.2414 

.2997 

4500 

.2547 

.3174 

5000 

.2685 

.3384 

neously,  we  use  the  constant  values  for  the  Poisson's  ratios  given  In 
the  SoRI  ATJ-S  graphite  data  book  [3-1].  All  of  the  foregoing  modeling 
Is  performed  with  the  Invaluable  aid  of  the  JNMDATA  program  described  in 
Section  2.4. 


3. 2. 2. 2 Elastic  Stress  and  Strain  Predictions 


The  elastic  stresses  and  strains  are  predicted  with  the  SAAS  IIIM 
finite  element  stress  analysis  computer  program  which  is  a Jones-Nelson- 
Morgan  model  version  of  the  SAAS  III  program  [3-9],  Elastic  predictions 
are  obtained  In  the  first  Iteration  of  the  Jones-Nel son-Morgan  model. 
Specifically,  the  stresses  and  strains  are  predicted  with  the  elastic 
mechanical  properties,  i.e.,  the  values  of  A in  Eq.  (2.3)  and  Table  3-6. 
The  finite  element  mesh  Is  composed  of  398  elements  with  454  nodal  points 
as  shown  In  Fig.  3-17.  Some  element  numbers  of  special  interest  In  sub- 
sequent discussions  are  labeled. 


The  elastic  predictions  with  the  Jones-Nel son-Morgan  model  are 
different  from  elastic  predictions  with  any  other  commonly  used  material 
model.  This  situation  occurs  because  the  actual  high  Initial  slopes  are 
more  closely  approximated  with  the  Jones-Nel son-Morgan  model  than  by  "eye- 
ball" from  a stress-strain  curve.  Thus,  comparison  of  predicted  elastic 
stresses  and  strains  with  other  results  may  not  be  fruitful.  At  any  rate, 
the  elastic  strains  e^,  ea,  e_,  and  are  displayed  In  contour 
plot  form  in  Figs.  3-18  through  3-22. 


FIGURE  3-18  e - ELASTIC 
r 
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FIGURE  3-20  ez  - ELASTIC 


3. 2. 2.3  Nonlinear  Stress  and  Strain  Predict! on s 


The  nonlinear  stresses  and  strains  are  obtained  by  iteration  of  the 
SAAS  IIIM  program  until  the  energies  in  all  finite  elements  do  not  change 
more  than  a specified  small  percentage  (in  this  case^-%).  These  stresses 
and  strains  are  what  Is  predicted  to  actually  exist  in  the  nosetlp  under 
the  prescribed  environment.  Thus,  they  are  the  values  which  will  be  com- 
pared with  the  DOASIS  predictions  in  Section  3.2.3.  At  this  point,  we 
only  display  contour  plots  of  the  nonlinear  strains  er,  ce,  e2,  yrz,  and 
Gmax  ^9S‘  3“23  through  3-27. 

The  results  depicted  In  Figs.  3-23  through  3-27  are  obtained  after 
nine  Iterations  of  the  SAAS  IIIM  computer  program.  The  rate  of  conver- 
gence is  studied  for  two  specific  finite  elements:  element  232  along 
the  Inner  contour  of  the  nosetlp  in  Fig.  3-17  and  element  134  in  the 
hottest  region  on  the  outside  contour  of  the  nosetlp.  Values  of  repre- 


FIGURE  3-23  cr  - NONLINEAR 
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FIGURE  3-27  cmax  ~ NONLINEAR 
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sentative  stresses  and  strains  are  listed  in  Tables  3-8  and  3-9  for  the 
two  elements. 

The  stresses  and  strains  for  element  232  in  Table  3-8  appear  to  have 
converged  in  five  iterations  for  all  practical  purposes.  This  rapid  con- 
vergence occurs  because  the  stress-strain  behavior  at  a temperature  of 
1397°F  is  not  highly  nonlinear.  A reasonable  question  is  why  are  results 
obtained  for  nine  iterations  Instead  of  stopping  at  five  (or  six)  itera- 
tions? The  answer  is  related  to  the  fact  that  convergence  is  achieved 
in  the  SAAS  IIIM  program  only  when  every  single  element  has  a strain 
energy  change  of  less  than  1/2%  from  the  energy  in  the  previous  iteration. 
Thus,  we  would  expect  to  find  other  finite  elements  for  which  the  con- 
vergence is  slower  than  it  is  for  element  232.  Brute  force  of  tediously 
examining  the  convergence  in  all  elements  is  presently  the  only  way  to 
determine  the  elements  for  which  convergence  is  slowest.  However,  some 
good  estimates  can  be  made  as  to  which  elements  might  have  slower  con- 
vergence than  element  232.  For  example,  consider  element  134  on  the 
outer  nosetlp  contour  in  Fig.  3-17.  The  temperature  there  is  very  high 
(7768°F),  and  hence  the  stress-strain  behavior  Is  very  nonlinear.  In 
fact,  we  can  only  approximate  the  stress-strain  behavior  for  such  high 
temperatures.  The  convergence  of  the  stresses  and  strains  for  element 
134  in  Table  3-9  is  much  slower  than  for  element  232.  The  stresses  and 
strains  are  not  quite  different  enough  in  the  last  few  iterations  to  be 
the  governing  factors  in  the  overall  nosetlp  iteration  procedure.  How- 
ever, we  need  only  be  assured  that  elements  exist  for  which  convergence 
is  slower  than  for  element  232  to  justify  having  more  than  five  or  six 
iterations  in  a nosetlp  stress  analysis. 
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TABLE  3-8 


PREDICTED  STRESSES  AND  STRAINS  IN  ELEMENT  232 


or,psi 

a^.psi 

4950 

5879 

2359 

2775 

2798 

3300 

2674 

3149 

2705 

3188 

2697 

3177 

2699 

3180 

2699 

3179 

2699 

3180 

.00166 

.00166 

.00166 

.00166 

.00166 

.00166 


.00600 

.00600 

.00600 

.00600 


TABLE  3-9 

PREDICTED  STRESSES  AND  STRAINS  IN  ELEMENT  134 


am1n*ps1 


-19857 

-3900 

-5698 

-4966 

-5199 

-5116 

-5145 

-5134 

-5138 


VPSi 

c^.psl 

Og.PSl 

■11363 

-9235 

-25737 

-1780 

-1549 

-4305 

-2842 

-2728 

-7225 

-2460 

-2251 

-5916 

-2580 

-2402 

-6321 

-2538 

-2348 

-6178 

-2553 

-2367 

-6227 

-2547 

-2360 

-6210 

-2549 

-2362 

-6216 
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3.2,3  COMPARISON  OF  JONES-NELSON-MORSAN  AND 


DOASIS  STRESS  AND  STRAIN  PREDICTIONS 
The  elastic  and  nonlinear  strains  plotted  In  contour  form  in  Figs. 
3-18  through  3-27  are  quite  similar  in  appearance  to  corresponding 
DOASIS  results  obtained  by  Baker,  Starrett,  and  Budde  [3-12],  The  only 
apparent  differences  are  somewhat  higher  strains  for  SAAS  HIM  than  for 
DOASIS  and  a difference  in  sign  for  yrz- 

The  SAAS  IIIM  sign  convention  for  shear  stress  and  shear  strain  is 
opposite  to  that  In  DOASIS,  This  difference  was  noted  after  the  rather 
expensive  computer  runs  were  made  so  the  wrong  contours  were  requested 
for  SAAS  IIIM.  That  Is,  a request  for  the  same  contours  as  plotted  from 
DOASIS  resulted  in  only  those  contours  which  were  common  to  both  sign 
conventions.  For  example,  yT1  contours  .005,  0.,  -.U05,  -.010,  and  -.015 
were  requested,  but  in  the  SAAS  IIIM  sign  convention  only  .015,  .010, 

.005,  0.,  and  -.005  exist  in  the  results.  Thus,  only  results  fur  .005, 
0.,  and  -.005  are  plotted.  Contours  for  .010  and  .015  do  exist  as  is 
easily  verified  by  Inspection  of  the  computer  output  listings,  but  simply 
are  not  plotted.  These  comments  apply  to  both  the  elastic  vrz  In  Fig. 
3-21  and  to  the  nonlinear  y , in  Fig.  3-26. 

The  higher  SAAS  IIIM  strains  In  Figs  C-18  through  3-27  than  DOASIS 
strains  [3-12]  are  more  easily  quantified  if  we  examine  the  strains  for 
a specific  element.  We  choose  element  number  232  on  the  inside  surface 
in  Fig.  3-17.  This  element  has  the  highest  predicted  maximum  strain  in 
the  cool  region  of  the  nosetip  and  is  the  element  about  which  a clip-on 
extensomett.'  Is  centered  in  the  50  MW  tests.  The  SAAS  IIIM  elastic 
stresses  for  element  232  are  seen  in  Table  3-10  to  be  nearly  50%  higher 
than  the  DOASIS  stresses.  Moreover,  the  SAAS  IIIM  strains  in  Table  3-11 
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TABLE  3-10 


are  15  to  30%  higher  than  the  DOASIS  strains.  However,  comparison  of 
these  two  elastic  predictions  is  not  entirely  valid  because  the  SAAS  HIM 
elastic  properties  are  much  higher  than  the  DOASIS  elastic  properties 
(higher  elastic  properties  than  in  conventional  analyses  is  characteristic 
of  the  Jones-Nelson  and  Jories-Nelson-Morgan  nonlinear  material  models 
as  noted  in  Section  2.2.2).  Thus,  we  do  not  make  an  issue  of  the  lack 
of  comparison  of  SAAS  IIIM  and  DOASIS  elastic  results  because  they  should 
not  be  the  same. 

However,  the  SAAS  IIIM  and  DOASIS  nonlinear  stresses  and  strains 
are  comparable.  In  Table  3-12,  we  see  that  the  SAAS  IIIM  nonlinear 
stresses  are  about  13%  lower  than  the  DOASIS  nonlinear  stresses.  At 
the  same  time,  the  SAAS  IIIM  nonlinear  strains  in  Table  3-13  are  nearly 
20%  higher  than  the  DOASIS  nonlinear  strains  in  element  232.  The  signif- 
icance of  this  latter  conclusion  will  be  revealed  In  Section  3.2.4. 

The  linear  and  nonlinear  a -e„  and  a -e,  results  are  shown  for  ele- 

r r z z 

ment  232  from  both  the  DOASIS  program  and  the  SAAS  IIIM  program  with  the 
Jones-Nelson-Morgan  (ONM)  material  model  in  Fig.  3-28.  There,  the  DOASIS 
uniaxial  bilinear  stress-strain  curves  are  plotted  with  a medium  solid 
line  (and  the  extension  of  the  elastic  slope  is  shown  with  a long  dash  - 
short  dash  line).  Also,  the  actual  nonlinear  stress-strain  curves  are 
shown  with  a heavy  line  along  with  the  elastic  slope  (the  tangent  to 
the  stress-strain  curve  at  o*0)  as  a dashed  line.  The  DOASIS  elastic 
results  are  plotted  as  a circle  which  lies  to  the  left  of  the  elastic 
uniaxial  stress-strain  curve  because  of  Poisson  effects  in  the  multlaxlal 
stress  state.  The  DOASIS  nonlinear  results  are  plotted  as  a dot  which 
Is  similarly  to  the  left  of  the  uniaxial  stress-strain  curve.  The  JNM 
elastic  multiaxial  results  are  shown  as  open  squares  and  are  to  the  left 
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TABLE  3-12 

NONLINEAR  STRESSES  IN  ELEMENT  232 
CALCULATED  WITH  DOASIS  AND  SAAS  IIIM 


m 

SAAS  IIIM 
psi 

A 

3,133 

2,699 

m 

1,647 

1,327 

BS 

HI 

3,666 

3,180 

BE9 

3,257 

-2,821 

-13% 

6,667 

5,771 

-13% 

TABLE  3-13 

NONLINEAR  STRAINS  IN  ELEMENT  232 
CALCULATED  WITH  DOASIS  AND  SAAS  IIIM 


STRAIN 

DOASIS 

SAAS  IIIM 

m 

Er 

.00140 

.00166 

+ 19% 

e0 

.00026 

.00042 

+62% 

ez 

.00284 

.00330 

+16% 

Yrz 

.00576 

-.00685 

+ 19% 

.00509 

. 00600 

+18% 

max 
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.001  .002  .003  .004 


FIGURE  3-28  LINEAR  VS.  NONLINEAR  BEHAVIOR  AT  NOSETIP  EL.  232 


of  the  elastic  uniaxial  stress-strain  curve.  Finally,  the  JNM  nonlinear 
multlaxial  results  are  shown  as  solid  squares  and  are  to  the  right  of  the 
actual  nonlinear  stress-strain  curve  as  Is  typical  of  graphitic  materials 
with  biaxial  softening.  Obviously,  the  JNM  linear  and  nonlinear  results 
are  drastically  different  from  each  other  even  for  the  relatively  modest 
nonlinear  behavior  at  this  temperature.  More  drastic  differences  between 
linear  and  nonlinear  results  exist  for  elements  In  hotter  and  hence  more 
nonlinear  regions  of  the  nosetip.  Note  that  the  elastic  strains  are  larg' 
er  than  the  nonlinear  strains  for  this  element.  In  contrast,  the  elastic 
strains  are  less  than  the  nonlinear  strains  for  element  134  In  a hotter 
region  of  the  nosetip.  These  two  results  are  similar  to  the  results  for 
the  Inside  and  outside  elements  In  the  SoRI  thermal  stress  disk  test  anal 
ysls.  That  is,  sometimes  the  elastic  strains  are  higher  than  the  non- 
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linear  strains  and  sometimes  they  are  lower. 

The  DOASIS  results  were  obtained  in  ten  Iterations  in  contrast  to 
the  nine  iterations  required  for  the  same  problem  in  SAAS  HIM.  No 
claim  is  made  that  SAAS  HIM  is  "faster"  than  DOASIS  on  the  basis  of 
the  observed  one  less  iteration.  As  a matter  of  fact,  DOASIS  is  prob- 
ably faster  for  this  nosetip  problem  because  DOASIS  is  faster  per  itera- 
tion than  SAAS  HIM.  The  meaningful  comparative  measure  is  the  accuracy 
of  the  results,  and  this  issue  will  be  addressed  in  the  next  section. 

3.2.4  SUMMARY  OF  50  MW  NOSETIP  STRESS  ANALYSIS 

The  objective  of  this  50  MW  nosetip  stress  analysis  Is  to  compare 
the  JNM  model  results  with  those  obtained  with  the  DOASIS  computer  pro- 
gram. That  comparison  Is  made  in  Section  3.2.3  where  the  JNM  model 
strains  are  seen  to  be  about  20%  greater  than  the  DOASIS  nonlinear 
strains  in  element  232  about  which  the  extensometer  in  the  50  MW  tests 
is  centered.  The  scope  of  this  contract  does  not  include  a correlation 
of  JNM  model  predictions  with  the  actual  50  MW  nosetip  strain  measure- 
ments. However,  a cursory  correlation  will  be  performed  to  put  the  JNM 
model  predictions  Into  some  perspective. 

The  elongation  between  two  points  on  the  nosetip  inside  contour 
was  measured  in  the  50  MW  tests  in  order  to  estimate  the  strain  tangent 
to  the  inside  contour.  This  elongation  between  nodal  points  199  and  304 
was  measured  with  an  extensometer  centered  about  element  232  In  Fig.  3-17. 
For  run  R41-016,  the  elongation  was  .0028  inches  [3-12],  The  mechanical 
strain  is  obtained  by  dividing  the  elongation  by  the  gage  length  of  .36 
Inches  and  subtracting  aAT  = ,0031.  The  resulting  mechanical  strain  is 
.0047.  The  maximum  mechanical  strain  in  element  232  is  predicted  to  be 
.0051  (Baker,  Starrett  and  Budde  report  .0038,  but  .0051  appears  in  the 
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DOASIS  output  provided  to  SMU  by  Captain  Budde)  with  DOASIS  and  .0060 
with  SAAS  HIM.  However,  the  maximum  mechanical  strain  in  element  232 
Is  not  tangent  to  the  inner  contour  of  the  nosetip.  This  fact  Is  readily 
observed  when  we  realize  that,  although  principal  stress  direction  Is 
tangent  to  the  Inner  contour  because  the  shear  stress  Is  zero  at  the 
Inner  contour,  the  principal  strain  direction  Umax)  Is  not  necessarily 
In  the  same  direction.  That  Is,  for  orthotropic  materials  In  nonprincipal 
material  directions  (i.e.,  not  In  the  present  r-6-z  principal  material 
coordinates),  principal  stress  directions  do  not  coincide  with  principal 
strain  directions!  The  tangent  to  the  Inner  contour  Is  at  47.7°  to  the 
r-dlrectlon.  Also,  the  angle  to  the  maximum  principal  stress  Is  47.4° 
(47.3°  In  DOASIS),  but  the  angle  to  the  maximum  principal  strain  Is 
51.7°  (52.0°  In  DOASIS).  A strain  transformation  relation  for  the  correct 
angle  could  be  used,  but  a more  direct  comparison  Is  available  - compare 
the  predicted  elongations  themselves, 

The  elongation  between  elements  199  and  304  Is  predicted  with  the 
SAAS  IIIM  program  to  be  .0033  Inches  and  with  DOASIS  to  be  .0030  Inches. 
Both  values  are  obtained  from  the  equation 


AL  = Lf Inal  " initial 


where 


Lf1nal  ” {pR  + ur^304  ' + ur^l99-*2  + + uz^304  “ + uz^l99-*2) 


L 


Initial 


C^304  " ^199^  + ^"»na  “ ^iqq)  3 


1/2 


304  ‘l 99' 


Both  the  DOASIS  and  JNM  model  elongation  predictions  are  higher  than  the 


measured  value.  On  the  other  hand,  the  DOASIS  prediction  reported 

ma  X 
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by  Baker,  Starrett,  and  Budde  [3-12]  Is  less  than  the  estimated  measured 

em  whereas  the  JNM  model  e„„u  prediction  is  greater  than  the  estimated 

measured  value.  One  reason  for  the  strain  differences  might  be  the  fact 

that  the  temperature  of  1600°F  used  by  Baker,  Starrett,  and  Budde  is 

higher  than  the  element  temperature  of  1 397° F used  in  the  analysis.  The 

extensometer  is  exposed  only  to  the  temperature  at  the  inner  contour  of 

the  nosetip.  Thus,  the  extensometer  temperature  should  be  less  than  the 

adjacent  element  temperature.  A lower  temperature  would  lead  to  a lower 

aAT  in  the  estimated  measured  eMau  calculation  and  hence  to  a higher 

max  ’ max 

These  kinds  of  speculations  and  questions  can  be  cleared  up  only  with  a 
detailed  correlation  study  (which,  as  mentioned  before,  is  not  the  objec- 
tive of  the  current  contract). 

Basically,  the  Jones-Nel son-Morgan  model  nosetip  predictions  are 
higher  in  strains  and  lower  In  stresses  than  the  DOASIS  predictions. 

The  key  elements  of  difference  are  the  manner  of  representation  of  uni- 
axial stress-strain  curve  behavior  and  the  manner  in  which  uniaxial 
stress-strain  curves  are  used  to  obtain  stresses  and  strains  under  multi- 
axial  loading.  The  correctness  of  approximation' for  both  models  can 
be  measured  only  by  meaningful  comparison  with  experimental  results. 

How  the  Jones-Nel son-Morgan  material  model  has  been  applied  in  such  a 
comparison  with  experiment  is  described  in  the  next  section. 


3.3  SUMMARY  OF  GRAPHITIC  MATERIAL  MODELING 

The  Jones-Nel son-Morgan  nonlinear  material  model  has  been  validated 
in  a meaningful  hierarchy  of  comparisons  with  experimental  results.  First, 
the  basic  mechanical  property  versus  energy  equations  are  curve-fit  to 
data  from  uniaxial  tests  In  principal  material  directions.  In  this  fash- 
ion, the  basic  model  is  defined.  The  model  is  then  validated  by  compari- 
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son  of  predicted  and  measured  strains  under  the  following  loading  condi- 
tions: 

(1)  uniaxial  mechanical  loading  in  other  than  principal  material 
directions  [3-3,  3-6,  3-14]. 

(2)  biaxial  mechanical  loading  in  principal  material  directions 
[3-2,  3-6,  3-14]. 

(3)  biaxial  thermal  loading  in  principal  material  directions  (Sec- 
tion 3.1 ). 

(4)  multlaxlal  thermal  and  mechanical  loading  In  other  than  principal 
material  directions  (Section  3.2). 

The  predicted  and  measured  strains  or  deformations  were  In  very  good  to 
excellent  agreement  In  the  first  three  loading  conditions.  The  agree- 
ment was  not  particularly  good  for  the  nosetlp  problem  In  condition  four, 
but  a meaningful  correlation  of  the  data  was  out  of  scope  of  the  contract 
so  only  cursory  results  are  available.  A more  complete  effort  to  corre- 
late the  nosetlp  results  should  be  made.  Then,  the  graphite  thermal  and 
mechanical  stress  analysis  problem  for  nosetlps  would  be  well  under  con- 
trol because  of  this  modeling  work  and  the  failure  stress  and  failure 
strain  work  of  Crose  [3-15]  and  Batdorf  [3-16].  The  latter  efforts 
depend  on  the  present  modeling  efforts  because  of  the  need  for  accurate 
nonlinear  stresses  and/or  strains. 
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4.  MODELING  OF  CARBON-CARBON  MATERIALS 


4.1  INTRODUCTION 

The  nonlinear  modeling  work  for  carbon-carbon  materials  was  initi- 
ated with  this  contract.  However,  only  a few  solid  accomplishments  were 
achieved  within  ne  time  and  funding  constraints.  These  results  are 
reported  in  this  section  in  the  following  order.  First,  the  basic  char- 
acteristics of  carbon-carbon  are  reviewed  and  examined  relative  to  their 
Importance  for  nonlinear  material  modeling  in  Section  4.2.  Then,  the 
flexural  modulus  and  strength  studies  for  multimodulus  materials  are  dis- 
cussed in  Section  4.3.  Finally,  the  essential  future  modeling  work  is 
described  In  Section  4.4. 

4.2  CHARACTERISTICS  OF  CARBON-CARBON 

Carbon-carbon  Is  not  a single  material,  but  is  Instead  a broad  class 
of  materials.  Many  types  of  carbon-carbon  are  manufactured  including 
various  kinds  of  woven  fibers  or  felt  materials  both  of  which  are  coated 
or  Impregnated  with  a matrix  material  In  a vapor  deposition  process  or  a 
pressure  Impregnation  process. 

The  steps  In  a typical  pressure  Impregnation  process  are  Illustrated 
In  Fig.  4-1.  In  the  first  step,  woven  graphite  fabric  In  the  horizontal 
plane  is  pierced  with  metal  rods  In  the  vertical  direction  while  being 
built  up  layer  by  layer  In  the  vertical  direction.  Next,  the  metal  rods 
are  replaced  with  graphite  yarn  or  graphite  fibers.  This  somewhat  loose 
assembly  of  layers  Is  Impregnated  with  phenolic  resin  In  step  three.  The 
impregnation  process  takes  place  first  by  evacuation  of  the  chamber  in 
which  the  material  Is  located  and  then  by  pressurization  to  force  the 
phenolic  resin  into  the  voids  of  the  woven  and  pierced  fabric.  Then,  In 
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FIGURE  4-1  MOD-3  FABRICATION  PROCESS 

step  four,  the  material  Is  cured  in  an  oven  at  3P0°F.  Next,  the  material 
Is  carbonized  at  2000°F  In  a kiln  In  step  five.  Finally,  the  material  Is 
graphltlzed  to  4900°F  In  an  Induction  furnace  in  step  six.  At  this  stage, 
the  billet  of  carbon-carbon  does  not  have  the  degree  of  density  possible, 

1 . e . , not  all  the  voids  are  filled  with  a graphi tl zed  form  of  the  phenolic 
resin.  Accordingly,  steps  three  through  six  are  repeated  until  the  de- 
sired density  is  attained. 

Carbon-carbon  Is  made  with  several  different  geometries.  Each  of 
the  layers  In  Fig.  4-1  can  be  either  woven  in  the  horizontal  (x-y)  plane 
as  In  Fig.  4-2a  or  made  of  parallel  fibers  as  In  Fig.  4-2b.  Avco  3D  Is 
of  parallel  fiber  construction  and  Is  shown  schematically  In  Fig.  4-3. 


./f  A*  -.v* 


a,  ORTHOGONALLY  REINFORCED 
PARALLEL  YARNS 


b.  ORTHOGONALLY  REINFORCED 
WOVEN  FABRIC 


FIGURE  4-2  GEOMETRY  OF  LAYERS  IN  x-y  PLANE 


FIGURE  4-3  AVCO  3D  CONSTRUCTION 
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One  variation  on  a three-dimensional  weave  or  three-dimensional  construc- 
tion of  orthogonal  fibers  is  to  add  fibers  at  45°  angles  to  the  x-,  y-, 
and  z-directions.  The  resulting  "7D"  construction  is  shown  schematically 
in  Fig.  4-4  for  fibers  with  prismatic  cross  sections.  A much  simpler 


BIAt 


FIGURE  4-4  PACKING  MODEL  OF  PRISMS  HAVING  EQUAL  CROSS  SECTIONAL  AREA 
IN  7-D  CUBIC  GEOMETRY 

carbon-carbon  structure  Is  obtained  after  carbon  vapor  deposition  on  car- 
bonized felt.  The  fibers  of  uncarbonized  felt  are  shown  in  the  scanning 
electron  photomicrograph  in  Fig.  4-5  where  obviously  the  fibers  have  no 
preferred  orientation.  The  scale  of  all  these  possible  carbon-carbon 
constructions  is  revealed  in  Fig.  4-6  where  the  approximate  fiber  spac- 
Ings  in  7-D  carbon-carbon  are  shown  In  a nosetip.  Obviously,  the  micro- 
scale of  carbon-carbon  materials  is  not  negligible  In  comparison  to  the 
nosetip  dimensions.  Generally,  the  microscale  dimensions  of  graphite, 
namely  the  particle  size,  etc.,  are  negligible  for  nosetip  stress  analysis 
problems.  However,  the  characteristic  material  dimensions  of  carbon- 
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FIGURE  4-5  UNCARBONIZED  VISCOSE-RAYON  FELT 


FIGURE  4-6  7-U  CARBON-CARBON  NOSETIP  WITH 
REPRESENTATIVE  FIBER  SPACINGS 
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carbon  are  not  clearly  negligible  for  nosetip  problems  although  perhaps 
the  answer  to  this  question  is  dependent  on  the  specific  carbon-carbon 
material  to  be  considered. 

The  anisotropy  of  the  various  carbon-carbon  materials  depends  pri- 
marily on  the  geometry  of  the  fiber  construction  and  secondarily  on  the 
manner  of  incorporating  a carbon  matrix.  For  example,  the  orthogonally 
reinforced  carbon-carbons  with  parallel  fibers  such  as  AVCO  3D  (Figs. 

4-2a  and  4-3)  and  with  woven  fabric  such  as  AVCO  MOD-3  (Fig.  4-2b)  are 
orthotropic  with  principal  material  directions  in  the  three  fiber  direc- 
tions. In  contrast,  GE  7-D  Is  not  only  orthogonally  reinforced,  but  also 
has  fibers  at  45°  to  the  x-,  y-,  and  z-directions  in  which  the  orthogonal 
fibers  are  placed.  Accordingly,  GE  7-D  is  anisotropic  (although  with  so 
many  fiber  directions,  the  degree  of  anisotropy  for  GE  7-D  is  lower  than 
the  degree  of  orthotropy  for  AVCO  3D,  i.e.,  GE  7-D  is  more  like  an  iso- 
tropic material  than  is  AVCO  3D).  The  manner  of  incorporating  a carbon 
matrix  in  carbon-carbon  can  be  an  influencing  factor  for  anisotropy  only 
when  the  process  of  constructing  or  depositing  the  matrix  has  some  direc- 
tional dependence.  For  example,  a matrix  formed  by  chemical  vapor  deposi- 
tion in  a shell  structural  element  can  have  a directional  dependence  of 
mechanical  properties  because  of  vapor  penetration  perpendicular  to  the 
shell  surface  and  no  penetration  parallel  to  the  surface.  This  effect 
is  similar  to  effects  obtained  with  pyrolytic  graphite  because  of  align- 
ment of  particles  due  to  direction  preferential  deposition  of  material. 

In  summary,  the  most  obvious  clue  to  the  degree  of  anisotropy  remains  the 
fiber  geometry,  but  anisotropy  is  also  influenced  by  the  matrix  construc- 
tion. 
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At  least  two  difficulties  arise  in  analysis  of  carbon-carbon  mate- 
rials, The  first  obvious  difficulty  is  the  analytical  complication  due 
to  the  many  more  mechanical  properties  necessary  to  characterize  carbon- 
carbons  as  compared  to  simpler  materials.  This  analytical  complication 
is  handled  with  relative  ease  in  comparison  to  the  more  troublesome  prob- 
lem of  measuring  the  additional  mechanical  properties.  Basically,  the 
problem  is  that  simple  strain  states  do  not  result  from  simple  Imposed 
stress  states.  For  example,  a uniaxial  stress  at  45°  to  the  fiber  direc- 
tions in  Fig.  4-7  leads  to  both  extension  in  the  direction  of  the  load 
(and  the  usual  Poisson  contraction  perpendicular  to  the  load)  and  In- 
plane shearing  if  the  ends  of  the  specimen  are  not  restrained  in  any  way. 
On  the  other  hand,  if  the  specimen  ends  are  restrained  to  remain  perpen- 
dicular to  the  load,  then  a complicated  nonuniform  shearing  and  exten- 
slonal  response  results.  The  complicated  response  occurs  if  the  test 
specimen  is  short  and  wide  whereas  the  simpler  response  occurs  if  the 
test  specimen  is  long  and  narrow.  In  fact,  If  the  specimen  Is  relatively 
short  and  wide,  the  stiffness  being  measured  is  not  the  Young's  modulus 


NO  END  EFFECT  RESTRAINED  ENDS 

FIGURE  4-7  DEFORMATION  OF  A UNIDIRECTIONALLY  REINFORCED  LAMINA 

LOADED  AT  45°  TO  THE  FIBER  DIRECTION 
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in  the  x-direction  of  the  sketch  In  Fig.  4-8,  but  the  two-dimensional 
(transformed  reduced)  stiffness  . The  reason  for  this  discrepancy  Is 
that  the  geometrically  admissible  state  of  strain  in  the  specimen  depends 
strongly  on  the  geometry.  If  the  specimen  is  long  and  slender,  then  the 
boundary  conditions  at  the  specimen  end  grips  are  of  no  consequence  5 la 
Saint  Venant.  Accordingly,  a pure  uniaxial  strain  is  obtained  and 

°x  ■ E* «,  o-n 

However,  for  a short,  wide  specimen,  the  end  restraint  of  ax  t 0, 

a Yvu  3 0 leads  to  the  stress-strain  relation 
y *y 

°x  ■ *n  EX 

Similar  results  are  shown  in  Fig.  4-8  for  an  applied  shear  stress  result- 
ing In  either  Gxy  as  desired  or  *66  If  the  test  specimen  geometry  is  not 


FIGURE  4-8  STIFFNESSES  Qn  AND  ^ VERSUS  MODULI  Ex  AND  G/y 
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The  foregoing  situation  Is  but  one  example  of  many  complicated 
mechanics  problems  which  arise  In  the  design  of  "proper"  test  specimens. 

By  "proper"  test  specimens  is  meant  specimens  for  which  the  desired  re- 
sponse can  be  measured  without  undue  influence  of  some  unwanted  response. 
For  example,  even  with  end  restraint  In  Fig.  4-7,  a specimen  can  be  made 
long  enough  that  the  resulting  deformation  looks  like  the  unrestrained 
end  response.  From  a practical  standpoint,  If  the  response  can  be  made 
predominantly  simple  by  suitable  choice  of  specimen  geometry,  then  the 
goal  of  a proper  test  specimen  Is  attained.  Analysis  of  specimen  re- 
sponse for  various  geometries  to  determine  the  geometry  for  which  simple 
response  occurs  Is  the  objective  In  many  mechanics  efforts.  However,  to 
date,  little  has  been  done  to  apply  the  principles  of  mechanics  in  ra- 
tional treatment  of  carbon-carbon  materials. 

Some  of  the  difficulties  In  rationally  analyzing  the  mechanical  be- 
havior of  carbon-carbon  materials  are  related  to  even  more  complicated 
response  characteristics  than  just  discussed  for  a unidirectlonally  rein- 
forced lamina.  For  example,  anisotropic  carbon-carbon  under  shear  loading 
In  one  principal  material  plane  will  also  have  shear  deformation  In  the 
two  principal  material  planes  which  are  perpendicular  to  the  stressed 
principal  material  plane.  This  shear  coupling  obviously  results  In  very 
complicated  strain  response.  Moreover,  the  multimodulus  characteristic 
of  carbon-carbon  materials  Is  yet  another  complicating  factor  In  addition 
to  the  nonlinear  stress-strain  behavior  in  every  principal  material  direc- 
tion. 

If  we  Ignore  for  a moment  both  the  multimodulus  character  and  the 
nonlinear  character  of  carbon-carbon  materials,  the  strain-stress  equa- 
tions for  linear  elastic  behavior  of  orthotropic  materials  are  (In  ortho- 
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gonal  x,  y,  z coordinates  aligned  with  the  principal  material  directions): 


Obviously,  nine  independent  constants 

E , E , E - Young's  moduli  in  x-,  y-,  and  z-directlons,  respec- 
x y z 

tively 


■ Poisson's  ratio  for  transverse  strain  in  the  j- 
dlrectlon  when  stress  exists  in  the  1-dlrectlons, 

1 .e., 

vij  * ■ 4 

for  cr,|  » o and  all  other  stresses  are  zero. 

G , G , G » shear  moduli  In  the  y-z,  z-x,  and  x-y  planes,  respec- 

y z z a 

tively. 

must  be  measured. 

We  now  recognize  that  the  strain-stress  relations  in  Eq.  (4.3)  must 


45  45  45 
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Jones  [4-1],  the  shear  moduli  are  replaced  with  Young's  moduli  at  45° 
to  principal  material  directions.  We  then  have  the  following  nine  inde- 
pendent mechanical  properties: 

/ic  a c A e 

(4.4) 

when  the  carbon-carbon  composite  has  three  orthogonal  principal  material 
directions  with  unequal  fiber  volumes  or  unequal  fiber  sizes  in  the  three 
directions.  If  the  fiber  volumes  or  fiber  sizes  are  equal  in  two  direc- 
tions (as  In  AVCO  3D  shown  In  Fig.  4-3),  the  six  independent  mechanical 
properties  are: 

Ex*  Ey  = Ex’  Ez’  vyz*  vzx  * vzy*  vxy*  Eyz’  Ezx  " Ezy’  Exy 
Finally,  If  equal  fiber  volumes  or  fiber  sizes  exist  in  all  three  direc- 
tions, the  three  independent  mechanical  properties  are: 

A C 

(4.6) 

(Note,  however,  that  because  of  the  three  fiber  directions  the  modulus 
at  45°  Is  Independent  of  the  direct  moduli  and  Poisson's  ratios  because 
the  shear  behavior  is  Independent.)  All  of  the  properties  In  Eqs.  (4.4), 
(4.5),  and  (4.6)  must  be  measured  both  In  tension  and  In  compression. 
Moreover,  these  properties  must  be  measured  at  many  stress  levels  to 
account  for  the  nonlinear  stress-strain  behavior  In  the  manner  of  the 
Jones-Nelson-Morgan  material  models  discussed  In  Section  2. 

At  this  point,  we  have  not  yet  examined  the  equations  for  off-axis 
behavior  of  such  orthotropic  materials  nor  have  we  considered  any  behavior 
of  anisotropic  carbon-carbon  materials.  We  defer  treatment  of  these  com- 
plicated questions  to  subsequent  reports  on  carbon-carbon  material  behav- 
ior. 
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Instead,  we  examine  orthotropic  and  anisotropic  carbon-carbon  me- 
chanical behavior  which  Is  representative  of  that  found  in  References 
4-2  through  4-11.  Those  references  are  the  primary  body  of  information 
on  carbon-carbon  mechanical  properties  although  other  Information  exists 
in  the  Journal  of  Composite  Materials  and  in  other  classified  and  unclas- 
sified reports.  Typical  stress-strain  curves  are  shown  in  Figs.  4-9  and 
4-10  for  the  materials  examined  in  the  CCAP  program  (Carbon-Carbon  Assess- 
ment Program).  All  of  the  stress-strain  curves  displayed  have  some  degree 
of  nonlinearity  ranging  from  mild  to  strong.  Notice  how  the  shapes  of  the 
actual  stress-strain  curves  are  similar  to  the  basic  shape  of  the  Jones- 
Nel son-Morgan  nonlinear  material  model  stress-strain  curve  In  Fig.  2-21  a 
on  p.  67.  The  Jones-Nel son-Morgan  model  was  developed  to  treat  carbon- 
carbon  behavior  in  addition  to  ATJ-S  graphite  as  well  as  boron/epoxy, 
graphite/ epoxy,  and  boron/aluminum.  However,  the  model  has  not  been 
applied  to  carbon-carbon. 

The  multimodulus  characteristic  of  carbon-carbon  is  a manufacturer- 
dependent  phenomenon.  Early  materials  had  drastic  diffc, dices  in  tension 
and  compression  moduli.  For  example,  the  MDAC  low  modulus  block  course 
weave  material  [4-2]  has  a ratio  of  tension  modulus  to  compression  modulus 
In  the  z-directlon  ranging  from  4 to  5!  More  recent  materials  have  closer 
tension  and  compression  moduli  - for  example,  AVCO  MOD  3 [4-5]  has  a ten- 
sion modulus  greater  than  the  compression  modulus  by  about  60 % In  both 
the  x-dlrectlon  and  the  z-directlon. 

Another  behavioral  characteristic  worthy  of  note  is  the  classical 
viscoplastic  behavior  of  carbon-carbon  at  high  temperatures.  For  example, 
AVCO  MOD  3a  exhibits  this  behavior  at  5000°F  in  Fig.  4-11.  The  stress 
rate  is  increased  by  factors  of  ten  and  the  stress-strain  curve  rises  in 
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FIGURE  4-9  PROBABLE  VALUE  TENSION  STRESS-STRAIN  CURVES  FOR 
CCAP  MATERIALS  AT  70°F  IN  THE  Z-DIRECTION  [4-6] 
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FIGURE  4-10  PROBABLE  VALUE  COMPRESSION  STRESS-STRAIN  CURVES  FOR 
CCAP  MATERIALS  AT  5000°F  IN  THE  Z-DIRECTION  [4-6] 
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FIGURE  4-11  COMPRESSION  STRESS-STRAIN  CURVES  FOR  AVCO  MOD  3a 
AT  5000° F IN  THE  Z-DIRECTION 
FOR  VARIOUS  STRESS  RATES  [4-5] 
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the  usual  viscoplastic  manner. 

Carbon-carbon  has  a characteristic  of  noncyl indrical  orthotropy  which 
is  quite  different  from  the  transverse  isotropy  of  ATJ-S  graphite.  The 
practical  significance  of  this  difference  is  that  in  a reentry  vehicle 
nosetip  application,  ATJ-S  graphite  Is  axisymmetric  whereas  carbon-carbon 
is  not.  This  Important  distinction  Is  significant  from  the  point  of  view 
of  the  nosetip  stress  analyst,  but  will  be  Ignored  In  this  report. 

In  summary,  carbon-carbon  materials  have  mechanical  behavioral  char- 
acteristics which  are  a substantial  escalation  in  degree  of  complication 
over  the  characteristics  of  graphite  materials.  The  carbon-carbon  analy- 
sis problem  is  very  difficult  and  complex.  Hence,  the  progress  will  be 
slow  because  of  the  many  mechanics  problems  that  must  be  addressed. 

4.3  APPARENT  FLEXURAL  MODULUS  AND  FLEXURAL  STRENGTH  OF 

MULTIMODULUS  MATERIALS 
4.3.1  INTRODUCTION 

The  ASTM  flexure  test  [4-12]  Is  commonly  used  to  measure  the  flexural 
modulus  and  flexural  strength  of  materials  other  than  the  plastics  for 
which  It  was  originally  developed.  The  officially  designated  equations 
for  reporting  the  test  results  are  applicable  only  to  materials  that  are 
linear  elastic  to  failure.  However,  the  test  is  often  used  for  carbon- 
carbon  and  other  composite  materials  which  have  different  (and  sometimes 
nonlinear)  stress-strain  curves  under  tension  loading  than  under  compres- 
sion loading. 

In  the  ASTM  flexure  test  for  plastic  materials  [4-12],  a rectangular 
cross  section  beam  Is  subjected  to  either  3-polnt  or  4-polnt  transverse 
loading  as  shown  In  Fig.  4-12.  The  slope  of  the  measured  load-deflection 
curve  is  related  to  the  apparent  flexural  modulus  which  is  defined 
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FIGURE  4-12  ASTM  FLEXURE  TEST  LOADING  SETUP 

with  the  usual  moment-curvature  relation: 

M » Ef  I Vq  (4.7) 

If  the  tension  modulus  of  the  material  Is  the  same  as  the  compression 
modulus,  then  the  flexural  modulus  measurement  Is  just  a (redundant)  way 
of  measuring  Young's  modulus.  The  flexural  failure  stress  Is  related  to 
the  failure  load  through  the  simple  equations  [4-12]: 

, 3-point  loading 

2bd‘ 

Or  B < (4.8) 

max  pL 

-=7  , 4-point  loading 

L bd2 

If  the  tension  strength  of  the  material  Is  the  same  as  the  compression 
strength,  then  the  flexural  strength  measurement  is  ju^.  a (redundant) 
way  of  measuring  the  strength.  However,  most  materials  have  different 
strengths  In  tension  than  In  compression  even  when  the  tension  and  com- 
pression moduli  are  the  same.  Accordingly,  the  apparent  flexural  strength 
is  the  lesser  of  the  tension  and  compression  strengths.  If  the  moduli 
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are  different  in  tension  than  in  compression,  then  the  flexural  strength 
calculated  with  the  ASTM  equations,  Eq.  (4.8),  is  neither  the  compres- 
sion strength  nor  the  tension  strength.  Thus,  the  Interpretation  of  the 
apparent  flexural  strength  is  difficult  for  many  materials.  The  key 
limitation  of  the  ASTM  flexure  test  is  its  applicability  only  to  linear 
elastic  materials  with  the  same  modulus  under  tension  loading  as  under 
compression  loading. 

The  ASTM  flexure  test  is,  however,  commonly  used  for  materials  such 
as  fiber-reinforced  composite  materials  which  have  a different,  modulus 
Et  under  tension  loading  than  the  modulus  Ec  under  compression  loading. 
These  multimodulus  materials  typically  have  the  stress-strain  behavior 
depicted  In  Fig.  4-13.  There,  the  actual  nonlinear  stress-strain  behav- 
ior is  approximated  with  a bilinear  stress-strain  curve  for  which  Et  is 
different  from  Ec<  For  this  behavior,  the  axial  strain  in  a beam  under 
bending  varies  linearly  through  the  depth  of  the  beam  as  in  Fig.  4-14. 


WITH  DIFFERENT  MODULI  IN  TENSION  AND  COMPRESSION 
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LINEAR  STRAIN 


BILINEAR  STRESS 


FIGURE  4-14  STRESS  AND  STRAIN  VARIATION  FOR  A BEAM  SUBJECTED  TO  MOMENT 


Moreover,  the  stresses  vary  In  a bilinear  fashion  through  the  depth  of 
the  beam.  The  significant  factor  is  that  the  neutral  axis,  i.e.,  the 
point  of  zero  stress  and  zero  strain,  is  not  at  the  mid-depth  of  the  beam. 
Thus,  because  of  the  mixed  state  of  tension  and  compression,  the  flexural 
modulus  for  multimodulus  materials  is  neither  nor  E_,  but  must  be 
interpreted  differently.  Also,  the  flexural  strength  Is  neither  the  ten- 
sion strength  or  the  compression  strength. 

The  objective  in  this  section  Is  to  assess  the  meaning  of  the  ASTM 
flexure  test  when  applied  to  multimodulus  fiber-reinforced  composite 
materials.  The  effect  of  different  moduli  in  tension  and  compression 
on  the  apparent  flexural  modulus  is  predicted.  The  predicted  flexural 
moduli  are  then  compared  with  measured  values  for  two  carbon-carbon  com- 
posite materials,  Sandia  CVD  carbon  felt  and  Avco  3D.  Sandla  CVD  carbon 
felt  [4-13]  Is  a chemical  vapor  deposited  (CVD)  carbon  matrix  in  carbon 
felt  and  Is  an  orthotropic  material  because  of  the  directional  nature  of 
vapor  deposition.  Avco  3D  [4-14]  has  a phenolic  matrix  with  high  modulus 
graphite  (Thornel  40  and/or  Thornel  50)  fibers  as  woven  reinforcement  in 
one  plane  and  quartz  filaments  as  reinforcement  in  an  orthogonal  plane. 
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Thus,  Avco  3D  is  an  orthotropic  material.  The  effect  of  different  moduli 
as  well  as  different  strengths  in  tension  and  compression  on  the  apparent 
flexural  strength  is  also  predicted,  and  the  results  are  compared  with 
measured  values  for  Sandia  CVD  carbon  felt  and  Avco  3D. 

The  emphasis  : this  section  is  on  proper  analysis  of  the  apparent 

flexural  modulus  and  strength  obtained  for  mul tiinodulus  materials  with 
the  ASTM  flexure  test.  An  important  purpose  in  this  section  is  to  identi- 
fy deficiencies  in  the  current  ASTM  flexure  test,  to  stimulate  work  to 
resolve  these  deficiencies,  and  to  account  for  all  pertinent  behavioral 
aspects  of  multimodulus  composite  materials.  The  results  of  this  analysis 
will  enable  a materials  scientist  to  determine  whether  experimental  data 
are  consistent  In  the  same  manner  that  the  usual  reciprocal  relation 
v12/,E1  = v21,/E2  usec*  t0  c^ec^  experimental  data  for  linear  elastic 
orthotropic  materials.  The  usefulness  of  this  effort  is  also  related  to 
getting  what  we  can  out  of  present  multimodulus  material  data  and  to 
generating  appropriate  data  for  multimodulus  materials  In  the  future. 

4.3ft  APPARENT  FLEXURAL  MODULUS 

Consider  a beam  of  a mul timodul us  material  subjected  to  pure  moment 
as  in  Fig,  4-14.  The  axial  strain  is  a linear  function  of  y and  depends 

II  II 

on  the  curvature  of  the  neutral  surface,  vQ  (v  evaluated  at  y » 0): 


e = vQ  y 


(4.9) 


The  uniaxial  stress-strain  relation  is 


a = ' 


e,  a > 0 
Ec  e,  a < 0 


(4.10) 
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We  first  locate  the  neutral  axis  since  it  is  not,  in  general,  at  the 
beam  mid-depth.  The  beam  is  subjected  only  to  moment  so  the  resultant 


axial  force  must  be  zero: 


■/. 


(4.11) 


Upon  substitution  of  the  stress-strain  relation,  Eq.  (4.10),  and  the 
strain-displacement  relation,  Eq.  (4.9),  the  force  on  a rectangular  cross 
section  beam  is 


rn  „ r "dn  „ 

J Et  vQ  y b dy  - J Ec  vQ  y b dy  = 0 


(4.12) 


where  d is  the  distance  from  the  bottom  of  the  beam  to  the  neutral  axis 
n 

as  defined  in  Fig.  4-14.  After  integration. 


[Etdn  - Ec'dV_  * 


(4.13) 


e’«  d„  = E 1/2  (d-d) 
t n c n 


(4.14) 


whereupon  the  neutral  axis. is  located  at 


dn  = 7 1/2  1/2  d 

Lt  Lc 


(4.15) 


Next,  the  bending  moment  about  the  neutral  axis  is 


n a n C n H p 

Et  v0  y b dy  + J Ec  vQ  y b dy  (4.16) 

A 
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which  can  be  Integrated  to  obtain 


II 

„ vo  b 


[Vi  + Ec'd-dn>33 


M « \^d1  + E_(d-dJJ  I (4.17) 

Now  substitute  the  location  of  the  neutral  axis,  Eq.  (4.15),  to  get 


M 


" 3 

v bdJ 
o 


J hvl  \\t(  hvl  ' 


(4.18) 


which  can  be  simplified  to 


11  3 


EtEc 


(4.19) 


Compare  this  moment  expression  with  Eq.  (4.7)  and  realize  that  I = bdJ/12 
for  a rectangular  cross  section  to  obtain 


4E*E 
t c 


(4.20) 


This  apparent  flexural  modulus  can  be  normalized  with  respect  to  the 
compression  modulus  in  the  form 


it > 

Ec  ' 1 + Z(Et./Ec>1^  + Et/Ec 


(4.21) 


and  plotted  along  with  the  normalized  average  modulus 


(Et  + EJ/2 


E, 


= (Et/Ec  + D/2 


(4.22) 
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and  the  normalized  tension  and  compression  moduli  versus  the  multimodulus 
ratio  Et/Ec  in  Fig.  4-15,  There,  we  see  that 

(1)  is  always  less  than  Eavg  and 

(2)  E^  is  always  closer  to  the  smaller  of  and  Ec> 

These  observations  can  be  used  to  test  the  quality  of  experimental 
flexural  modulus  data.  That  is,  if  the  experimental  data  do  not  exhibit 
those  characteristics,  then  we  would  have  to  be  suspicious  of  their  valid- 
ity. This  test  Is  analogous  to  applying  the  reciprocal  relation  v12/E1  = 
v21^E2  t0  test  exPerimental  modulus  and  Poisson's  ratios  data  for  ortho- 
tropic materials. 

Experimental  results  for  materials  with  substantially  different  mod- 
uli In  tension  and  compression  are  plotted  In  Fig.  4-16  which  is  an  expan- 
sion of  Fig.  4-15  near  the  region  where  Et/Ec  ■ 1.  For  Sandla  CVD  carbon 
felt  [4-13],  an  orthotropic  material,  the  tension  modulus  is  as  much  as 
42%  higher  than  the  compression  modulus  depending  on  the  vapor  deposition 
pressure  in  Table  4-1.  For  Avco  30  [4-14],  the  tension  moduli  are  up  to 
25%  lower  than  the  compression  moduli  as  shown  in  Table  4-2.  Although 
Avco  3D  is  an  orthotropic  material,  the  orthotropy  Is  excited  only  in  the 

form  of  G since  the  material  In  the  beam  is  oriented  such  that  E * Eu 
xy  x y 

(1.e.»  the  x-  and  y-coordinates  of  the  beam  coincide  with  the  principal 
material  directions  in  the  plane  of  the  graphite  fibers).  The  flexural 
moduli  are  calculated  from  Eq . (4.20)  with  the  known  tension  and  compres- 
sion moduli  and  are  also  shown  in  Tables  4-1  and  4-2.  There,  the  calcu- 
lated moduli  are  5%  and  7%  too  high  and  3%  too  low  for  Sandla  CVD  carbon 
felt  and  2%  and  5%  too  high  for  Avco  3D.  This  agreement  between  theory 
and  experiment  Is  rather  good  considering  the  inherent  accuracy  of  the 
experiments  and  the  bilinear  stress-strain  curve  approximation.  Thus, 
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FIGURE  4-16  NORMALIZED  FLEXURAL  MODULUS 

- EXPERIMENTAL  AND  THEORETICAL  RESULTS 
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DEPOSITION 
PRESSURE 
mm  Hg 


TABLE  4-1 

STIFFNESSES  OF  SANU1A  CVD  CARBON  FELT 


Ec 

lOSpsi 

(GN/m2 


Ef 

MEASURED 

106ps1 

(GN/m2) 

CALCULATED 

106psl 

(GN/m2) 

ERROR 

M— 

HuuHHfi&fli 

+5.0% 

B&fiHHSSfH 

+6.7% 

4.0 

(27.6) 

BOScSI 

-2.5% 

TABLE  4-2 

STIFFNESSES  OF  AVCO  3D  THORNEL/PHENOLIC 


THORNEl. 

VOLUME 

FRACTION 


MEASURED 


CALCULATED 

lO^psI 

(GN/m2) 


a .i|r( . j u <*» 


we  can  claim  the  ability  to  calculate  the  flexural  modulus  for  these 
materials  without  having  to  resort  to  the  expense  and  time  of  an  ASTM 
flexure  test  (not  to  mention  the  high  cost  of  the  material  for  the  test 
specimens! ) . 

4.3.3  APPARENT  FLEXURAL  STRENGTH 

Irrespective  of  whether  3-point  or  4-point  loading  is  used,  the  ASTM 
maximum  stress  versus  moment  relation  is 


ASTM, 


max 


Me 

I 


Md/2 

bdVl2 


6M 

17 


(4.23) 


In  which  the  neutral  axis  is  presumed  to  be  at  the  mid-depth  of  the  beam. 
If  the  maximum  moment  versus  load  relation  for  each  type  of  loading  is 
considered,  then  the  ASTM  maximum  stress  versus  load  relations  are 


P . L 

? 


ASTM, 


= < 


max 


P . L 
3 


PL 


, 3-point  loading 


4-point  loading 


(4.24) 


For  multimodulus  materials,  we  will  derive  the  correct  maximum  stress 
versus  moment  relation  in  which  the  difference  between  the  neutral  axis 
and  the  beam  mid-depth  Is  accounted  for. 

The  maximum  stresses  in  the  beam  cross  section  are 


a*  = E.  e*  a E*  v.  y+  * E.  v d 
t t ...  t o t ton 


max  max 


max 


ac  = Ec  ec  " Ec  vo  yc  = Ec  vo  (d_dn) 
max  max  max 


(4.25) 
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Irrespective  of  whether  Et  Ec  or  Et  Ec . We  solve  for  vQ  from  Eq. 
(4.19)  and  substitute  for  dn  from  Eq.  (4.15)  to  get 


. c 3M 

°t  ' Et 
max  c bd 


¥c 


Et1/? 

TrrrrjTn 

,Et  + Ec 


3M. 

'max  v bd" 


°c_.  = Ec  73  i 


(4.26) 


EtEc 


. ■> 
which  can  be  simplified  to 

r, 

(Et/Ec)1/2  + 1 
>1/2 


Lt 


1/2 


°t  = < 

max 


a “ < 

max 


1 + (E^)1 


? S 6M 
2 } — j- 

b<r 


2 

bd^ 


(4.27) 


where  the  terms  in  braces  are  the  corrections  due  to  noncoincidence  of 
the  neutral  axis  and  the  beam  mid-depth.  Obviously,  when  Et  a Ec,  the 
stresses  In  Eq.  (4.27)  are  equal  to  each  other  and  to  the  ASTM  maximum 
stress  In  Eq.  (4.23). 

The  actual  maximum  tensile  and  compressive  stresses  for  multimodulus 
materials  In  the  ASTM  flexure  test  can  be  summarized  In  the  form 


o\  a C* 
tmax  1 maxASTM 


o„  B C a 
cmax  c ™XASTM 


(4.28) 


where  the  correction  factors  Ct  and  Cc  are  multipliers  on  the  calculated 
ASTM  stresses  used  to  obtain  the  actual  maximum  stresses.  These  correc- 
tion factors  are  plotted  in  Fig.  4-17  as  a function  of  the  multimodulus 
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FIGURE  4-17  ACTUAL  MAXIMUM  TENSILE  AND  COMPRESSIVE  STRESSES 


VERSUS  MULTIMODULUS  RATIO 
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15%  lower  than  the  ASTM  value. 


Note  that  simultaneous  failure  In  tension  and  compression  is  not  im- 
plied in  these  calculations.  That  is,  neither  ot  nor  a is  neces- 

cmax  max 

sarily  the  respective  strength  St  or  S . Thus,  if  we  are  given  the  load 
or  moment  at  which  a beam  falls,  we  calculate  the  maximum  tensile  and 
compressive  stresses  from  Eq.  (4.27)  and  compare  them  with  the  respective 
strengths.  However,  the  analysis  is  restricted  to  materials  which  can 
be  approximated  with  the  bilinear  stress-strain  curve  In  Fig.  4-13  and 
which  are  linear  elastic  to  failure.  This  latter  restriction  is  probably 
the  toughest  to  satisfy  for  most  materials. 

Four  subcases  of  behavior  occur  based  on  whether  Et  is  larger  or 
smaller  than  Ec  and  whether  St  is  larger  or  smaller  than  $c: 

O)  Et>Ec,  St>Sc 

(2)  E,  > Ec,  St<Sc 

(3)  Et<Ec,  St>Sc 

t4>  Et<Ec-  St<Sc 

In  both  cases  (1)  and  (2),  <jt  Is  underpredicted  and  a is  ovarpre- 

nnax  cmax 

dieted  with  the  ASTM  equations,  Eq.  (4.8),  as  seen  in  Fig.  4-18.  On  the 

other  hand,  in  cases  (3)  and  (4),  ot  is  overpredicted  and  a is 

tmax  Sax 

underpredicted.  In  case  (2),  the  beam  fails  in  tension  because  ot  > 

Sax 

o and  St  < S . Thus,  the  apparent  flexural  strength  is 
max 


S = — 


(4.29) 


which  is  always  less  than  the  tension  strength  because  Ct  > 1 . In  case 

(3),  the  beam  fails  in  compression  because  o„  > u.  and  S < S. . 

max 


Thus,  the  apparent  flexural  strength  is 


■ ■<**-. A-  j Hi  *.  m»n t # *»■  -■  . »*  .... 


(4.30) 


which  is  always  less  than  the  compression  strength  because  Cc  > 1 . No 
definite  failure  mode  exists  for  the  inequalities  in  cases  (1)  and  (4); 
instead,  the  apparent  flexural  strength  is 

S,  S 

S,  = lowest  of  and  (4.31) 

T 4 4 

which  is  always  between  the  tension  and  compression  strengths.  Obviously, 
the  apparent  flexural  strength  never  coincides  with  either  St  or  Sc  for 
multimodulus  materials.  In  fact,  is  sometimes  lower  than  St  [case  (2)], 
sometimes  lower  than  Sc  [case  (3)],  and  sometimes  between  St  and  Sc  [cases 
(1  )■  and  (4)]. 

Sandla  CVD  Carbon  Felt  has  Et  > Ec  and  Sc  > St  as  seen  In  Table  4-3 
so  Is  a case  (2)  material  which  always  fails  In  tension.  The  apparent 
flexural  strengths  are  calculated  from  Eq.  (4.29)  and  listed  In  Table 
4-3.  Note  that  the  calculated  flexural  strengths  are  below  the  values 
of  St  (as  we  calculated  them  to  be),  but  the  measured  flexural  strengths 
are  above  the  Sfc  values. 

Avco  3D  has  E_  > E*  and  S_<  S,.  as  seen  In  Table  4-4  so  Is  a case 
(3)  material  which  always  falls  in  compression.  The  apparent  flexural 
strengths  are  calculated  from  Eq.  (4.30)  and  listed  In  Table  4-4.  How- 
ever, the  calculated  values  are  lower  than  Sc,  but  the  measured  values 
of  5^  are  even  higher  than  the  St  values  (which  In  turn  are  higher  than 
the  Sc  values). 

The  measured  behavior  Is  significantly  different  from  the  predicted 
behavior  for  both  example  materials.  Thus,  we  must  question  either  the 
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TABLE  4-3 

STRENGTHS  OF  SANDIA  CVD  CARBON  FELT 


DEPOSITION 
PRESSURE 
mm  Hg 

St 

103ps1 

(MN/m2) 

Sc 

103ps1 

(MN/m2) 

Et 

*7 

sf 

MEASURED 

103ps1 

(MN/m2) 

CALCULATED 

103ps1 

(MN/m2) 

35 

4.78 

(33.0) 

23.5 

062) 

1.10 

6.18 

(42.6) 

4.67 

(32.2) 

100 

6.55 

(45.2) 

35.5 

(245) 

1.32 

9.07 

(62.5) 

6.10 

(42.0) 

760 

6.26 

(43.2) 

51 .8 
(357) 

1.42 

9.26 

(63.9) 

5.89 

(40.6) 

TABLE  4-4 

STRENGTHS  OF  AVCO  3D  THORNEL/PHENOLIC 


THORNEL 

VOLUME 

FRACTION 

St 

103ps1 

(MN/m2) 

Sc 

103ps1 

(MN/m2) 

Et 

Sf 

MEASURED 

103ps1 

(MN/m2) 

CALCULATED 

103ps1 

(MN/m2) 

50 

19.3 

(133) 

13.3 

(91.7) 

.855 

22.1 

(152) 

12.8 

(88.3) 

66 

23.6 

(163) 

17.6 

(121) 

.750 

30.9 

(213) 

16.3 

(112) 
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theory,  its  applicability  to  the  data,  or  the  data  themselves.  Presuming 
the  data  are  correct,  we  examine  the  theory  to  see  if  any  deficiencies 
exist.  The  most  obvious  deficiency  of  this  theory  as  applied  to  ortho- 
tropic materials  is  the  neglect  of  shear  behavior.  However,  Sandia 
CVD  carbon  felt  has  only  mild  orthotropy  and  has  the  same  kind  of  dis- 
agreement between  theory  and  experiment  as  Avco  3D  so  a shear-related 
deficiency  is  not  likely. 

Perhaps  another  possible  reason  for  the  disagreement  is  that  the 
stress-strain  behavior  of  both  example  materials  is  not  linear  elastic 
to  failure.  Thus,  since  the  theory  is  not  applicable  to  nonlinear  stress- 
strain  behavior,  it  is  not  applicable  to  the  strength  analysis  of  these 
materials.  However,  the  ASTM  prescribed  maximum  stress  calculations  are 
not  applicable  to  these  kinds  of  materials  either.  Thus,  the  meaning  of 
the  reported  flexural  strength  data  is  not  clear.  What  Is  clear  is  that 
the  ASTM  flexure  test  and  calculations,  as  presently  prescribed,  are  not 
applicable  to  multimodulus  materials  or  to  materials  with  nonlinear  stress- 
strain  behavior.  However,  these  expensive  tests  continue  to  be  run  and 
reported  without  sufficient  information  for  proper  interpretation  and 
utilization. 

4.3.4  SUMMARY 

Flexural  modulus  and  flexural  strength  data  for  multimodulus  mate- 
rials cannot  be  used  to  extract  what  happens  Individually  in  tension  and 
in  compression.  Instead,  flexure  behavior  is  simply  some  kind  of  averag- 
ing of  those  two  behaviors  with  shear  effects  thrown  in  for  confusion. 
Neglecting  shear  is  ignoring  an  Important  behavioral  aspect  of  most  com- 
posite materials.  However,  the  Inclusion  of  shear  effects  in  the  flexure 
test  analysis  may  not  be  warranted  because  flexure  data  are  not  basic 
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data  (i.e.,  not  independent  of  other  parameters)  and  therefore  need  not 
be  refined.  Accordingly,  the  ASTM  flexure  test  is  very  misleading  when 
applied  to  materials  with  distinctly  different  stiffnesses  and  strengths 
in  tension  than  in  compression.  Moreover,  the  typically  nonlinear  stress- 
strain  behavior  of  fiber-reinforced  composite  materials,  especially 
carbon-carbon,  is  another  invalidating  factor  for  the  ASTM  flexure  test. 
Designers  who  seem  to  want  these  flexure  data  are  actually  avoiding  or 
masking  the  real  problems  of  nonlinear  stress-strain  behavior  as  well 
as  sometimes  low  and  sometimes  high  shear  modulus  and  shear  strength  of 
many  composite  materials.  (In  fact,  the  span-to-depth  ratio  for  the 
flexure  test  specimen  is  probably  so  different  from  the  actual  designer's 
application  that  shear  plays  entirely  different  roles  in  the  specimen 
and  the  design  application.)  The  proper  course  of  action  is  to  develop 
and  use  more  rational  analyses  of  nonlinear  behavior  and  shear  effects 
coupled  with  better  shear  modulus  and  shear  strength  data  that  are  re- 
ported only  in  conjunction  with  the  pertinent  stress-strain  curves  to 
fail ure. 

4.4  FUTURE  MODELING  WORK 

The  Jones-Nel son-Morgan  nonlinear  material  model  must  be  fit  to 
representative  carbon-carbon  materials.  Subsequently,  the  model  must 
be  exercised  in  various  theoretical -experimental  correlation  efforts  to 
validate  the  model. 

The  first  step  for  each  material  is  to  fit  the  model  to  the  mechani- 
cal behavior  in  the  principal  material  directions  and  at  45°  to  the  prin- 
cipal material  directions  for  an  orthotropic  material  (more  complicated 
mechanical  behavior  must  be  treated  for  anisotropic  materials).  The  only 
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practical  way  to  perform  these  material  model  fits  is  to  use  the  JNMDATA 
computer  program  (discussed  in  Section  2.4),  to  ease  the  burden  of  plot- 
ting the  actual  and  approximate  stress-strain  and  mechanical  property 
versus  energy  curves.  However,  before  the  JNMDATA  program  can  be  uti- 
lized, it  must  be  extended  to  orthotropic  materials  and  possibly  to 
anisotropic  materials.  These  extensions  are  relatively  straightforward. 
With  the  extensions,  the  JNMDATA  program  will  be  an  essential  element 
in  the  procedure  for  modeling  carbon-carbon  properties. 

The  next  logical  step  in  applying  the  material  model  is  a series 
of  validations  of  the  model  by  comparison  of  predicted  and  measured 
strains  for  a representative  set  of  loading  conditions  ranging  from 
simple  to  complex  and  including  both  mechanical  and  thermal  loads.  The 

simplest  such  situation  Is  uniaxial  mechanical  loading  at  some  angle  to 

45 

the  principal  material  directions  (other  than  the  45°  angle  at  which  E 
Is  measured.  The  next  simplest  loading  condition  is  a biaxial  stress 
state  such  as  occurs  In  the  MDAC  tubular  test  specimen  for  graphite  [4-15]. 

Various  carbon-carbon  materials  have  been  examined  in  the  Southern 
Research  Institute  thermal  stress  disk  test  which  is  described  for  graphi- 
tic materials  in  Section  3.2.  However,  the  response  of  a carbon-carbon 

disk  is  not  axlsymmetric  because  of  the  orthotropic  character  of  carbon- 

/ 

carbon,  i . e. , the  disk  does  not  have  cylindrical  orthotropy  much  less 
the  in-plane  Isotropy  of  a graphite  disk.  Accordingly,  the  analysis  of 
a carbon-carbon  disk  is  a substantial  escalation  in  difficulty  over  that 
of  a graphite  disk.  For  example,  the  Jones-Nelson-Morgan  nonlinear  mate- 
rial model  could  be  incorporated  in  the  ASAAS  finite  element  program  [4-16] 
with  modification  for  noncyl indrlcal  orthotropy  or  some  three-dimensional 
finite  element  computer  program.  Obviously,  such  analyses  would  Involve 
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considerable  escalation  in  cost  for  computer  time  over  the  cost  for 
graphitic  material  analyses. 

Jortner  obtained  preliminary  measurements  of  strains  in  carbon-carbon 
disks  under  a biaxial  stress  field  due  to  high  speed  rotation  [4-17]. 
Experimental  results  such  as  these  could  be  used  in  yet  another  correla- 
tion of  theoretical  and  experimental  results. 

The  next  logical  step  in  a hierarchy  of  correlations  between  theory 
and  experiment  is  to  consider  the  thermal  and  mechanical  loading  for  an 
actual  reentry  vehicle  nosetlp.  The  level  of  analysis  needed  Is  identical 
to  that  for  the  thermal  stress  disk  test  because  of  the  noncyl indr leal 
orthotropy  of  carbon-carbon.  Accordingly,  the  nosetlp  analyses  would  be 
quite  expensive  relative  to  those  for  graphitic  materials. 

Obviously,  the  foregoing  series  of  correlations  between  theory  and 
experiment  Is  much  more  complicated  than  the  series  just  completed  for 
ATJ-S  graphite  (completed  with  the  exception  of  the  nosetlp  strain  corre- 
lation), Moreover,  because  of  the  relatively  high  expense  and  variability 
of  the  mechanical  properties  of  carbon-carbon,  the  basic  mechanical  proper- 
ties will  not  be  available  with  the  statistical  confidence  of  graphitic 
materials.  Thus,  there  is  no  room  in  the  continuing  development  of  carbon- 
carbon  for  time,  scarce  material,  and  money  for  tests  such  as  the  flexural 
modulus  and  strength  test  which  are  not  justified  or  useful.  We  must 
strive  to  obtain  the  most  and  highest  quality  Information  from  each  and 
every  measurement. 
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5.  CONCLUDING  REMARKS 


The  major  graphitic  material  analysis  problem  areas  of  biaxial 
softening  and  multimodulus  nonlinear  behavior  have  been  successfully 
treated  with  the  various  Jones-Nel  son-Morgan  nonlinear  material  models. 
The  models  have  been  validated  by  comparison  of  predicted  response  with 
measured  response  for  a series  of  meaningful  physical  problems.  These 
problems  include  uniaxial  off-axis  and  biaxial  mechanical  loading  as  well 
as  biaxial  thermal  loading  of  an  annular  disk  and  multiaxial  thermal  and 
mechanical  loading  of  a reentry  vehicle  nosetlp.  The  capability  to 
accurately  predict  stresses  and  strains  in  ATJ-S  graphite  under  complex 
mechanical  and  thermal  loading  has  been  achieved. 

On  the  other  hand,  the  application  of  the  Jones-Nel  son-Morgan  non- 
linear material  model  to  carbon-carbon  materials  is  just  beginning.  The 
model  is  apparently  well-suited  for  treatment  of  the  stress-strain  re- 
sponse typical  of  carbon-carbon  materials.  However,  this  supposition 
must  be  verified  in  a meaningful  series  of  correlations  between  Jones- 
Nelson-Morgan  material  model  predictions  and  carefully  measured  response 
for  various  mechanical  and  thermal  loading  problems. 
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APPENDIX 


DETERMINATION  OF  THE  POINT  OF  ZERO  SLOPE 
ON  AN  IMPLIED  STRESS-STRAIN  CURVE  BY  INTERVAL  HALVING 

The  point  of  zero  slope  on  an  implied  stress-strain  curve  or  the 
point  at  which  the  slope  of  the  stress-strain  curve  has  some  prescribed 
nonzero  value  must  be  determined  before  the  linear  extrapolations  of  the 
stress-strain  data  described  in  Section  2.3.3  can  be  implemented.  The 
expression  for  the  slope  of  an  Implied  stress-strain  curve  in  Eq.  (2.20) 
is  a function  of  the  stress  and  the  strain.  However,  for  a given  strain 
the  corresponding  Implied  stress  is  easily  determined  from  the  implied 
stress-strain  relation  In  Eq.  (2.19)  so  the  slope  is  a function  of  the 
strain  only.  The  strain  e*  at  which  the  slope  of  the  Implied  stress- 
strain  curve  is  zero  or  a prescribed  nonzero  value  is  the  root  of  a 
complicated  nonlinear  equation  obtained  by  combining  Eqs.  (2.19)  and 
(2.20).  This  root  Is  easily  found  by  the  numerical  technique  of  inter- 
val halving.  Interval  halving,  In  general,  Involves  halving  an  Interval 
In  which  the  root  of  a nonlinear  expression  lies  until  the  root  is  en- 
closed by  a very  small  interval.  This  technique  is  described  in  this 
appendix  for  finding  the  point  of  zero  slope,  but  the  general  procedure 
Is  the  same  for  finding  a point  of  prescribed  nonzero  slope. 

Before  the  interval  halving  technique  can  be  applied,  an  interval 
in  which  the  strain  e*  lies  must  be  determined.  This  interval  is  found 
by  choosing  an  initial  strain  at  which  the  slope  of  the  implied  stress- 
strain  curve  is  positive  and  then  increasing  this  strain  by  prescribed 
Increments  until  a strain  level  is  reached  at  which  the  slope  of  the  Im- 
plied stress-strain  curve  is  negative,  For  example,  at  the  last  stress- 
strain  data  point  in  Fig.  A-1,  the  slope  of  the  implied  stress-strain 
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FIGURE  A-l  STRAIN  INCREMENTS  FOR  FINDING  POINT 

OF  ZERO  SLOPE  ON  A IMPLIED  STRESS-STRAIN  CURVE 

curve  Is  positive.  The  point  of  zero  slope  Is  located  to  the  right  of 
this  point.  Thus,  the  strain  eQ  corresponding  to  the  last  data  point 
In  Fig.  A-l  Is  used  as  the  Initial  left  bound  of  the  interval  In  which 
e*  lies.  Then,  the  strain  eq  Is  Increased  by  a prescribed  increment  Ae, 
and  the  slope  Is  evaluated  at  this  new  strain  e-j . The  slope  of  the  im- 
plied stress-strain  curve  In  Fig.  A-l  is  positive  at  , so  eQ  Is  re- 
placed as  the  left  bound  of  the  interval  by  the  new  strain  . Next, 
e-|  is  Increased  by  Ae,  and  the  slope  is  evaluated  at  e2>  Again,  the 
slope  Is  positive  so  is  replaced  as  the  left  bound  of  the  Interval 
by  e2.  The  incrementing  process  Is  continued  until  finally,  at  the 
strain  The  slope  of  the  implied  stress-strain  curve  is  negative. 
Thus,  the  strain  e*  Is  located  in  the  interval  e2  < e < E3-  The  size  of 
the  prescribed  increment  Ae  used  in  finding  the  interval  surrounding  e* 
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is  arbitrary  but  should  be  chosen  wisely.  The  Increment  should  not  be 
so  small  that  the  strain  has  to  be  increased  a large  number  of  times  be- 
fore a point  of  negative  slope  is  found.  At  the  same  time,  the  increment 
Ae  should  not  be  so  large  that  the  interval  enclosing  the  point  of  zero 
slope  Is  exceedingly  wide.  However,  due  to  the  rapid  convergence  of  the 
method  after  interval  halving  is  initiated,, reasonably  large  increments 
are  preferable  over  very  small  Increments. 

The  interval  eg  < c < e3  surrounding  c*  In  Fig.  A-l  is  large  and 
must  be  decreased  until  the  difference  between  the  left  and  right  bounds 
of  the  Interval  is  very  small;  i.e.,  until  e*  is  known  to  a specified 
level  of  accuracy.  The  interval  around  the  point  of  zero  slope  is  de- 
creased in  an  orderly  manner  by  use  of  interval  halving.  The  Interval 
halving  technique  is  easily  described  with  the  aid  of  the  slope-strain 
curve  in  Fig.  A-2.  The  initial  Interval  surrounding  c*,  the  strain  at 
zero  slope,  is  bounded  in  Fig.  A-2  by  ep  and  tn,  the  strains  at  which 


FIGURE  A-2  INTERVAL  HALVING  OF  A SLOPE-STRAIN  CURVE 
TO  FIND  POINT  OF  ZERO  SLOPE 
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the  slope  Is  positive  and  negative,  respectively.  The  strains  cp  and 
in  Fig.  A-2  correspond  to  the  strains  c2  and  , respectively,  in  Fig. 

A-l . The  search  for  c*  Is  initiated  by  evaluating  the  slope  at  the  mid- 
point of  the  interval  bounded  by  c and  Ln>  This  initial  strain  is  la- 
beled c i in  Fig.  A-2.  If  the  slope  is  negative  at  e-|  as  in  Fig.  A-2  c* 
lies/between  cp  and  e:^.  If  the  slope  is  positive  at  , e*  lies  between 
l ^ and  en.  Thus,  Ej  lies  in  an  interval  half  as  large  as  the  original 
Interval.  The  slope  is  then  evaluated  at  the  midpoint  of  this  new  inter- 
val, i.e.,  at  eg  *n  Fig*  A-2.  The  slope  at  e2  is  positive  so  c*  lies  be- 
tween e i and  Eg.  Thus,  the  Interval  enclosing  e*  is  halved  again.  This 
interval  halving  process  is  continued  until  the  interval  around  c*  be- 
comes sufficiently  small;  i.e.,  until  c*  is  known  to  a desired  level  of 
accuracy. 

In  conclusion,  the  point  of  zero  slope  on  an  implied  stress-strain 
curve  is  quickly  and  easily  determined  by  the  numerical  technique  of  in- 
terval halving.  The  process  consists  of  two  steps.  First,  an  interval 
of  strain  in  which  the  strain  corresponding  to  zero  slope  is  located  must 
be  determined.  Then,  this  Interval  is  halved  until  it  becomes  very  small 
about  the  strain  corresponding  to  zero  slope.  A point  of  nonzero  slope 
on  the  implied  stress-strain  curve  can  be  found  by  a procedure  analogous 
to  the  one  for  finding  the  point  of  zero  slope.  The  only  difference  be- 
tween the  two  procedures  is  that  for  finding  a point  of  nonzero  pre- 
scribed slope,  the  strain  at  which  the  slope-strain  curve  crosses  the 
l ine  ~j~  = ( J7")  * t 0 is  desired  instead  of  the  stra  in  at  which  the  slope- 
strain  curve  crosses  the  slope  axis  ~ = 0. 
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